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NOMENCLATURE 


An asterisk before the symbols indicates an input variable. All input 
subscripted and non-subs cripted variables are listed separately. 


English 


Symbol 

Definition 

Units Used 


Major and minor semiaxis, respectively, 
of grain exterior used in the ovality 
analysis. 

in. 


Propellant burning rate coefficient above 
and below the transition pressure, re- 
spectively. 

In/sec-psi’ 

*a j. ,a 
c*p Yp 

Pressure sensitivity of C* and y» 
spectlvely. 

— 

*®c*T 

Temperature sensitivity of C* 

/*F 

a ,b 
g g 

Major and minor semiaxis of grain interior, 
respectively, used in ovality analysis. 

in. 

^n’\p ’^s 

Total burning surface associated with 
nozzle end, port and slot surfaces, 
respectively. 

in2 

*\nT 

Burning surfaces (as a function of y) at 
nozzle end for tabular input. 

in2 

*^)pT 

Burning surfaces (as a function of y) of 
port sides for tabular input. 

in^ 

*\sT 

2 . 

Burning surfaces (as a function of y) of slots in 

for tabular input. 

*\hT’*\nT 

Controlling port areas at head and aft ends , 
respectively (as a function of y) for 
purely tabular inputs. 

in^ 

C 

V 

Coefficient of variation; l.e., the ratio 
of the standard deviation to the mean. 

- — 
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English 

Symbol 


C*,*C * 
’ n 


D*,*Dj* 


t ref 


xh* yh 


*e ,*e 
xn* yn 


E ,*E , 

n’ ref 




Definition Units Used 

Characteristic exhaust velocity at ft/sec 

Tgj. and and at standard conditions, 
respectively. 
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respectively, of the nozzle throat. 

Diameter of the nozzle exit. in. 


Length average initial diameter of con- in. 
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Length average outside diameter of circular in. 

perforated grain, excluding lengths in 

closures. 
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nozzle throat erosion equation. 

Initial diameter of the thrust termination in. 
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interior with respect to the center of the 
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Thrust Ibf. 

Gravitational constant (ft/sec^)lbm/lbf 

Integer designating type of grain perforation. 

The half width of the first and second set in. 

of points, respectively, of a wagon-wheel grain. 

Estimated burnout altitude. ft. 

Axial length of nozzle closure. In. 
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English 

Symbol 

eo 


* 1 , 




op 


*I 


po 


^sp*^T 


*I I 
xi* xf 


K 




A 


f 


Definition Units Used 

Integer designating option of ovality 

and malalignment calculations. 

Number of steps in the integrations of 

perimeters in ovality analysis subroutine. 

Integer designating type of program input. 

Integer designating options of plotting 

results and obtaining special outputs. 

Specific and total impulse, respectively. lbf-sec/lbm,lbf-sec 

The initial and final seed numbers, re- 

spectively, for the random number generator. 

Criterion value in Pearson’s system for 

determining frequency curves. 

The length of the two parallel sides of the in. 

first and second set of points, respectively, 
of a wagon-wheel grain. 

Distance between center line of motor and in. 

fillet center of standard star grain. 


A 


s 


Length of sides of truncated star point in. 

excluding fillets. 


*A 


TP 


Initial length of termination passages between in. 
centers of gravity of perimeters of bases. 


*L 


Total initial perforated grain length including in. 
gaps at slots. 


Lg^,*Lg^^ Instantaneous and initial total axial lengths, in. 

respectively, of circular perforated grain 
(not including gaps) . 

*L^j^ Initial slant length of grain surface at the in. 

nozzle end. 

*Lg^j^ Initial total length of star-shaped perforated in. 

grain. 
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English 

Symbol Definition Units Used 



n 




*n 

n 


*n 

P 

*n 

s 


Estimated length of grain at the aft end at In. 

start of first talloff having an additional 
taper not represented by or 0q. 

Number of observations of a statistically 

distributed variable or burning rate exponent. 

Burning rate exponents above and below 
the transition pressure, respectively. 

Number of burning flat end surfaces of a 

star grain located at the extreme nozzle 
end of the chamber. 

Number of star points. 

Number of burning flat end surfaces of a 

star grain not located at the nozzle end 
of the chamber. 


*N(j) Integer designating whether or not a specific 

output plot is desired. 

P Pressure. Ibf/in^ 


max 


Maximum head end chamber pressure calculated psia 
by the program. 


*P 


ref 


*P 


txan 



r 


Reference average nozzle stagnation pressure psia 
used In the nozzle throat erosion equation. 

Transition pressure at which the burning rate psia 
coefficient and exponent change. 

Integer designating grain arrangement. 

Burning rate. In/sec 


Initial inside radius of the propellant web 
of a standard star, truncated star, and 
wagon-wheel grain, respectively. 




*R 

c 


*R. 


OAl 


Distance from center of curvature of a 
spherical end of circular perforated grain 
to the burning surface associated with 0q. 

Outside radius of a star grain. 

The propellant oxidizer to aluminum weight 
ratio. 


in. 


in. 


in. 
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NOMENaATURE (Continued) 


English 

Symbol 

Definition 

Units Used 

*R 

P 

Initial radius of truncated star points. 

in. 

*S 

Number of burning flat slot sides of a 
circular perforated grain not Including 
the nozzle end. 

— 

*S 

op 

Integer designating type of star grain. 

— 

t 

Time 

sec. 

"b 

Calculated total operating time of motor. 

sec. 

*^1 

Estimated time at burnout. 

sec. 


Ignition delay at 60 “F grain temperature. 

sec. 

*t 

max q 

Estimated time at which maximum dynamic 
pressure on the vehicle occurs . 

sec. 

Si 

Earlier time at which tailoff begins in a 
motor pair. 

sec. 

*T 

gr 

Bulk temperature of the propellant grain. 

»F 

*V 

clT 

Initial volume of chamber gases associated 
with tabular input . 

in^ 

W . ,W ,,W 
pi’ p2’ p 

Total masses of propellant burned based 
on mass discharge rate, volume calculations 
and the arithmetic average of and Wp 2 * 

Ibm. 

X 

Value of general statistical variable. 

Units vary 

X »y 

c*''c 

Coordinates of the grain exterior used in 
the ovality analysis. 

in. 

X »y 

g g 

Coordinates of the grain interior used in 
the ovality analysis. 

In. 

*x 

out 

Distance burned at which propellant breaks up 

. in. 

*Sa 

Difference in web thicknesses at ends of L_ . 

la 

in. 

**Tz 

Difference between the initial circular per- 
forated grain diameter and the nozzle end of 
Lcci and the nominal value of D^ less Zq and 

less twice x„ . 

la 

in. 

y 

Distance propellant has burned from initial 
surface. 

in. 
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English 

Symbol Definition Units Used 



Initial and instantaneous differences , re- 
spectively, between web thicknesses due 
to grain bore taper at head and nozzle ends 
of controlling grain length, excluding any 
initial length associated with or 0^. 

in. 

*z 

c 

Initial difference between web thicknesses 
due to grain exterior taper at the head and 
aft ends of the controlling grain length. 

in. 

Greek 

Symbol 



*a , ,*a 
ah an 

The angular orientation of the ovality of 
the grain Interior with respect to the grain 
exterior at the head and nozzle end of the 
grain, respectively. 

degrees 

*«1**“2 

The angle between the slant sides of a wagon- 
wheel grain point and the center line of the 
point for the first and second set of points, 
respectively. 

degrees 

*a , 
eb 

Erosive burning coefficient in the in^*®-ft*^ 

Robillard-Lenoir rule. 

• ®/sec^ • ®lbf 

*a 

n 

Nozzle exit half angle. 

degrees 

*3 

Erosive burning pressure coefficient in 
the Robillard-Lenoir rule. 

— 

3i,32 

Ratio of the square of the third moment 
about the mean to the cube of the second 
moment and the ratio of the fourth moment 
to the square of the second moment, re- 
spectively, for a statistically distributed 
variable. 


AP^/Ay 

Rate of change of chamber pressure with 
respect to distance burned. 

Ibf/in3 

*(AP/Ay) 

Depressurization rate at which propellant is 
extinguished for computation control purposes. 

lb f /in 3 

*AR , *AR 
ch* cn 

One half the difference between the maximum 
and minimum diameter of the grain exterior 
at the head and nozzle end reference planes , 
respectively. 

in 
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Greek 

Symbol 


*AR . ,*AR 
gh* gn 


Ay,*Ayj^ 



*8 

cn 




e 


fw 



*e 

n 

*e 

p 

e 

8 




De finition 


One half the difference between the 
maximum and minimum diameter of the grain 
bore at the head and nozzle end reference 
planes } respectively. 

Incremental distance burned and initial 
value of same, respectively. 

Thrust loss coefficient. 

Approximate acute angle bonded circular 
perforated grain makes with motor center 
line at the head and nozzle closures, 
respectively. Also referred to as the 
grain contraction angles. 

Angular location of fillet center of 
standard star from line of symmetry. 

Angular location of fillet centers with 
respect to radial center line of wagon- 
wheel grain points. 

Angle burning surface element of circular 
grain makes with longitudinal axis of motor 
at the nozzle end of the chamber. 

Nozzle cant angle. 

Angle of standard star point. 

Half angle of star grain sector. 

Angular location of slot side of truncated 
star grain. 

Acute angle between axis of thrust 
termination passage and motor axis. 


X 




Volumetric loading density; i.e., initial 
volume occupied by propellant divided by 
empty case volume. 

Second, third and fourth moments, respec- 
tively, of a statistically distributed 
variable about its mean. 


Units Used 
in. 

in. 

degrees 

degrees 

radians 

degrees 

degrees 

degrees 

radians 

radians 

degrees 
Units vary 
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Greek 

Symbol 


P K 



a 


cri,C2 


T 

S 


*T 

Teff 



ww 


Subscripts 

a 

c 

f 

g,G 

h 

1 

ig 

max 

min 

n 

o 

q 

t 


Definition 


Units Used 


Temperature sensitivity coefficient of 
pressure at constant K = A^/A*. 

Solid propellant density. 


/“F 


slugs/in^ 


The standard deviation of a statistically Units vary 

distributed variable; i.e., the square root 
of the second moment about its mean value. 


The square root of the second moment of a Units vary 

statistically distributed variable about zero 
(a^) and a^ly/T~y respectively. 

Thickness of propellant web at slot bottom in. 

of truncated star grain. 

Estimated "effective" web thickness of in. 

termination port. 

Web thickness of main propellant grain. in. 

Web thickness of standard star (same as in. 

except for some combination grains) . 

Web thickness of wagon-wheel grain. in. 


Value during web action time. 

Case, grain exterior or chamber value. 
Final . 

Grain interior. 

Head end of grain. 

Initial. 

Ignition. 

Maximum value. 

Minimum value. 

Nozzle or nozzle end of grain. 
Stagnation. 

Dynamic pressure. 

Value during talloff. 
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V Value at v/eb time. 

100k Value when thrust has decayed to 100,000 lbs. 

Superscripts 

* Choked throat value. 

- Arithmetic average value over Ay or statistical mean. 

• Time rate of change. 


xvil 



I. INTRODUCTION AND SUMMARY 

This report presents the results of research performed at Auburn 
Unviersity during the period May 31 to November 30, 1974, under modi- 
fication No. 6 to the Cooperative Agreement, dated October 6, 1969, 
between NASA Marshall Space Flight Center and Auburn University. The 
principal objective of the research was to develop a technique for 
statistically investigating the thrust imbalance of a pair of solid 
rocket motors (SRMs) firing in parallel. 

The study of thrust balance and imbalance is of particular interest 
with regard to application to the NASA Space Shuttle because two very 
large solid motors fire in parallel on the Shuttle. Although a similar 
arrangement was utilized in the Titan program with somewhat smaller motors, 
because of the differences in the configuration of the Titan and Shuttle, 
it is more imperative that the thrust on the solid motor booster stage 
be uniform to assure proper guidance and control for the overall vehicle. 

Past analyses of SRM reproducibility have been concerned mainly with 
characteristics of an entire population of single motors rather than pairs. 
Usually a non-time varying parameter such as average thrust and total 
impulse has been of interest, the distributions of variables affecting 
the parameter are assumed normal in the statistical sense, and cross- 
correlation effects are neglected. Reference 1 typifies such an analysis. 

For the present investigation the Monte Carlo technique (Reference 2) 
was selected. Sets of the significant variables are selected on a prob- 
ability basis using a random sampling technique and the imbalance calcu- 
lated for a large number of motor pairs using a mathematical model of the 
internal ballistics. This method is not limited to normal distributions 
of the input variables. Raw data, histograms, or cumulative distribution 
functions may be used for any of forty or more controlling variables to 
specify the statistical input. 

Errors arising from neglecting the cross-correlation of variables 
are minimized in this study by selecting for the most part completely 
independent variables. However, the analysis may be readily extended to 
account for the cross-correlations where they are shown to exist and 
when they are calculable. The imbalance is evaluated in terms of six time- 
varying parameters as well as eleven single valued ones which themselves 
are subject to statistical analysis. 

Application of the Monte Carlo method, which requires evaluation of 
a large number of motor pairs, is made practical and economical through 
the utilization of the simplified computer program for the internal bal- 
listics presented in References 3, 4, and 5. Familiarity of the reader 
with References 3 and 4 or alternatively v;ith 5, which is a consolidated 


- 1 - 
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report on 3 and 4, Is assumed. Although the simplified program makes use 
of numerous approximations, the effects of the vast majority of the 
variables affecting Imbalance are represented. Also, the essence of the 
phenomena controlling the critical tailoff portion of the thrust trace 
are embodied in the program which has been upgraded during the current 
effort to refine the analysis and to incorporate additional input 
variables which may now be statistically distributed. It Is fundamental 
to the approach that, in spite of bias in the performance calculated 
resulting from the approximations in the ballistic analysis, the bias 
will reflect equally in each motor of a pair. Thus the difference in 
performance calculations for a pair of SRMs will be of a higher degree 
of accuracy than the individual motor calculations assuming that the 
significant parameters causing the difference are incorporated into the 
model and their effects precisely determined. 

Owing to the complexity of the problem it was not possible to include 
all variables or give precise representation of the Influences of all 
those that are included.- Nevertheless, it is felt that the model developed 
yields a good first estimate of the difference in performance between a 
pair of SRMs. Early in the study, the sensitivities of performance to 
the variation in the various input parameters was evaluated by changing 
each variable one at a time a small amount between a pair of SRMs and 
conqjuting the performance of each motor. Partial results are given in 
the Appendix in the form of pairs of thrust-time traces obtained for each 
variable. These traces were obtained directly from the computer output using 
the CalComp plotter. Although only very minor performance differences 
were obtained in many cases, it was decided that for the Monte Carlo in- 
vestigation all but a few of the variables would be statistically determined 
in order to minimize cumxilative errors that might result from neglect of 
a large number of seemingly unimportant variables. 

Illustrative of the comprehensiveness of the analysis is the capa- 
bility which has been added to calculate as a computer program option 
the approximate effects of case and mandrel ovality, eccentricity, and 
malalignment. On the other hand, the effects of grain temperature gradient 
(as opposed to bulk temperature) variation between rocket motors are not 
within the capability of the program to evaluate. However, the ovality 
analysis would permit simulation of the approximate results of one special 
situation for temperature gradient from which insight into the effects of 
this important variable may be gleaned. 

The Monte Carlo computer program retains the capability of the pro- 
grams of References 3, 4, and 5 to treat segmented configurations with 
both star and circular perforated grains present in various arrangements 
in the same SRM as well as monolithic grains with either a circular per- 
foration or one of the three most common types of star grains. However, 
the program is more accurate if most of the grain is circular perforated 
as in the Space Shuttle. Also, whereas tabular input may still be used 
to specify portions of the burning surface, the statistical input vari- 
ables influence only the burning geometry that is computed from the program 
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equatlons. Basic assumptions of the program are that the propellant does 
not break up and is not extinguished except by being completely consumed 
by burning normal to the propellant surface. 

Ignition transients are not calculated in the analysis. However, 
the variations in the initial equilibrium chamber pressures as calculated 
by assuming the grains are completely ignited have been incorporated 
along with statistically determined values of ignition delay. It is 
believed that this treatment embodies the significant effects of ignition 
upon the remainder of the trace. 

Perhaps the greatest limitation of the analysis is its reliance on 
the availability of reliable statistical data to specify accurately the 
distributions of the numerous variables. For example, there simply have 
not been enough SRMs built in the size class of those in the Space Shuttle 
stage to provide direct statistical samples. For small motors, where large 
populations exist from which data may be directly applied or scaled to 
another SEM, documentation of manufacturing variations is often incomplete 
or unavailable. In the test cases used in this report, drawing tolerance 
limits were taken in a number of cases to represent the ranges correspond- 
ing to six standard deviations in assiomed normally distributed dimensions. 
Hopefully, this research will provide a stimulus for acquisition of 
statistical data which will reduce reliance on assumptions of this nature. 
As suggested earlier and as is demonstrated in the sample case , the com- 
puter program is virtually unlimited in its ability to treat the various 
types of statistical distributions. 

The program which is written in FORTRAI'I IV requires approximately 
2,000 computer cards not Including the data cards necessary to specify the 
input variables and their statistical distributions. The compilation time 
is approximately 1 minute and 10 seconds on the IBM 370/155 computer. Per- 
formance coinJutation time for 1 pair of SRMs using recommended increment 
sizes in the several integration processes involved is approximately 35 
seconds using the FORTRAN H compiler. 

Section II of this report discusses the changes made to the ballis- 
tic performance analysis program and also describes how the statistically 
distributed variables are treated and sets of random input variables 
selected. The basic mathematical details of the ovality and malalignment 
analysis are also given. A discussion of the program input and output 
is presented in Section III. Although most of the input variables are 
defined in References 3, 4, and 5, a complete discussion is given here 
for a ready reference and guide to the user in specification of his 
problem and the interpretation of the output. Section IV gives the infor- 
mation required to operate the program including data card format and the 
program listing which contains concise printed definitions of all input 
variables along v/ith identification of the statistically determined vari- 
ables to serve as a checklist. In Section V a sample study is presented 
to demonstrate the setup procedures, format and computational capabilities 
of the computer program. 



II. ANALYSIS 


This section of the report describes the changes and additions that 
have been made to the ballistic analyses of References 3, 4, and 5, and 
presents the basic elements of the Monte Carlo analysis. Rationale to 
the selection of input variables is interspersed throughout the discus- 
sion. The section is divided into four parts: the first two parts give 
the changes to the "main” program and the "area" subroutine, the third 
part gives the ovality and malalignment analysis and the final part is 
the statistical presentation. It is assumed the reader is familiar with 
References 3, 4, and 5, and thus in discussing program revisions only 
information sufficient to identify the changes and the bases for them is 
given. Complete documentation of the program logic and computational 
changes is given in the program listing in Section IV. Although a number 
of completely new subroutines have been incorporated into the program, 
the methods used are straightforward applications of basic statistics and 
analytical geometry. Therefore, again concise descriptions of the analyses 
and their capabilities and limitations are given in lieu of detailed math- 
ematical procedures which can be readily identified from the program 
listing. 

Internal Ballistics (Main Program) 

Inert parts and Ignition . The computation of the weights of inert 
parts have been removed from the program since the program is now 
strictly a performance analysis rather than a design analysis. Also, the 
option of calculating the ignition transient by the methods of References 
4 and 5 has been removed because the complications of the ignition tran- 
sients make a detailed Monte Carlo impractical from a computer requirements 
standpoint. However, an Ignition delay time has been incorporated as a 
statistical variable subject to both random variation and systematic 
variation due to the distribution of bulk temperature of the grain. Also, 
the initial equilibrium pressure is subject to variations arising from a 
number of sources. It is believed that this treatment of ignition embodies 
the significant effects of the Ignition transient on the remainder of the 
thrust time trace. 

Throat area . The throat relationship has been modified to account 
for Reynolds* number cind Stanton's number effects of pressure and diameter 
on heat transfer (Ref. 6) as follows: 


E = E - (P /P J0-®(D* 
n n ref on' ref ref' 

D* « D * + 2.0 ZE At 
1 n 

■ •trD*^/4 


- 4 - 


( 1 ) 


A* 


( 2 ) 

( 3 ) 


-5- 


Mach number . To account for throat erosion effects, the Mach number 
at the nozzle exit from which the thrust coefficient is ultimately derived 
is now recalculated at each time step rather than just initially. 

Characteristic velocity and ratio of specific heats . These variables , 

C* and Y» respectively, are calculated based on the oxidizer to aluminum 
ratio of the propellant which is introduced as a new independent 

distributed variable. In Wis way the cross-correlation between C* and y 
has been incorporated. Data for determining separate relationships for 
C* and Y as a function of Rq^^ should be obtained from thermochemical 
analyses. The present progr^ uses typical functional relationships obtained 
from a linear regression analysis of such data. Similar relationships must 
be established to fit the propellant system under consideration. More 
rigorously, additional composition components might be Included. However, 
because of the relatively strong influence of the aluminum on the propel- 
lant chemistry the two-component oxidizer to aluminum system should convey 
the largest effects of composition on C* and y* 


Systematic variations in C* and y have also been incorporated to give 
the effects of chamber pressure and propellant bulk temperature on C* and 
of chamber pressure on y. The relationships are as follows; 


=*ref ■ ) 

(A) 
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The detailed attention given to these variables was motivated not only 
by the cross-correlation but also by the high sensitivity of performance 
to both C* and y as revealed by the sensitivity study discussed in Section 
III. The effect of y is largely due to its influence on exit Mach number 
and hence thrust coefficient. In the subsonic flow regime, Mach number is 
relatively insensitive to variation in y and it is found convenient to neglect 
the variations in y resulting from pressure variations in calculating the 
Mach number at the end of the grain. 

The Introduction of Rq^^j, as an independent variable raises the question 
of possible cross-correlation with other independent propellant variables, 
especially the burning rate coefficient a and density p. When such cor- 
relations can be identified they may of course be treated in a manner 
similar to that given for C* and y. Although the matter has not been investi- 
gated fully as of this writing, it appears that factors such as oxidizer 
particle size distributions, amount of burning rate catalyst and other com- 
position variables will play a more important role than the oxidizer to 
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aluminum ratio in fixing the distributions of burning rate coefficient and 
density; cross-correlation between and p are assumed negligible as 

have, with somewhat less justification, the correlation between a and p, 
themselves. The latter warrants further investigation as it is intuitive- 
ly clear that at least a weak cross-correlation exists between a and p be- 
cause of the effects of oxidizer particle size distribution on these 
parameters . 

Burning rate coefficient . Bulk temperature has been introduced as an 
explicit variable subject to statistical variation. It modifies the burning 
rate coefficient: 


® ” ®60* ® * P [(l-n)(iTp)j^(Tg^-60)] (7) 

Also, the program now makes use of two sets of burning rate coefficients 
and exponents in recognition of the possibility of separate values applying 
above and below a specified transition pressure. In conjunction with this, 
the level above which burning of the propellant is permitted was changed 
from 30 psla to 5 psla as the thrust of some SRMs may still be significant 
at low pressure. Choked flow is still assumed at the lower pressure levels. 

Time increments . Several changes were made in the calculation of time 
increments for the purpose of obtaining more accurate representation of the 
thrust-time characteristics of the Individual SRM. First, the time incre- 
ment is now calculated from the burning rate and Incremental distance burned 
using the burning rate that applies during the Increment under consideration 
rather than the previous increment. This is done except for the purpose of 
computing changes in port areas which must be known before the rate is cal- 
culated. The error in the latter case is minor and a time consuming Iterative 
solution is avoided. Also, time values and corresponding output values are 
now obtained at the precise y positions that the lengths of grain associated 
with X-pa and z begin to burnout. These times are identified in the program 
printout by "initial tailoff begins" and "final tailoff begins," respectively. 
It was found desirable in general to use a y Increment (dy) of approximately 
0.001 of the web thickness. This gives good precision in the calculations 
and also sufficient values to obtain good graphical portrayals. However, the 
smaller Increment size introduces a difficulty in satisfying the mass flow 
continuity relationship during rapid pressure changes; the effect of such 
changes tends to be overestimated in the simplified ballistic model when small 
increment sizes are used. In References 3, A, and 5, this difficulty was 
avoided during tailoff by calculating the nozzle end stagnation pressure 
directly from the pressure gradient. In the new program, the same method is 
also applied to the period involving the burning of the length of grain 
associated with X^a as similar difficulties were encountered in the region 
when small increment sizes were used. 
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Internal Ballistics (Area Subroutine) 

During the present study, an error was found in the first expression 
for i^>nc Eqs. A4c of References 4 and 5. Specifically the first value 
of Rj in the equation should be replaced by 

I(D^/2)2 _ R^2]1/2 + I(D^/2)2 _ + y)2]l/2 

and R 3 should be replaced by Ry throughout the remainder of Eqs. A4c. 

The program listing contains the necessary changes to the computer pro- 
gram. It will be noted that additional changes applicable to Eq. A4c 
have been made to further refine the closure geometry effects for a 
circular perforated grain. Also, application of 4b and c has been modi- 
fied in the computer program to avoid the possibility of negative square 
roots in terms such as [Dq^ - (Di+ 2 y)^]' 2 . 

Ovality and Malalignment (Oval Subroutine ) 

Ovality and lack of concentricity of the grain perforation with 
respect to the motor case can clearly influence the ballistic performance 
of SRMs during the critical talloff period. As a first approach to the 
problem it is assumed that the burning surface geometry embodying the 
effects of ovality and malalignment may be defined by specifying three 
radial reference planes - one near the head end of the rocket, one at the 
aft end of Lpa ^o* reference planes must 

Intersect the cylindrical portion of the rocket motor case to eliminate 
the effect of end closures on the geometric properties to be calculated. 

The implied assumption here is that the central portion of the rocket 
dominates the influences of ovality and malalignment. 

The geometry of the referencf; planes is illustrated by Figure II-l. 

To fix ideas, hotli the. exte.rior and inte.rior (bore) surfaces of the grain 
are assumed lo (»e distor:te<i into oval sha]jed surfaces from nominally cir- 
cular ones. However, as discussed later, v/ith some restrictions and loss 
in accuracy, the analysis is also applicable to grains with the star-shaped 
perforations treated in the program. 

The exterior and interior grain perimeters are next, assumed to be 
elliptical and the interior grain surface to remain elliptical as burning 
progresses. The latter assumption is not rigorous for burning normal to 
the surface but the error introduced is insignificant for the small degrees 
of ovality to be encountered in practice. Thus, the ellipse defining the 
burning perimeter is given by Equation ( 8 ) where y is the distance burned 

X ^/(a + y)^ + y_^/(b + y)^ = 1 

So 00 


( 8 ) 
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and 3g and bg are semi-axes of the ellipse defining the initial burning 
grain perimeter. The exterior grain perimeter is expressed in the co- 
ordinates Xg and Yg of the interior perimeter by the equation 

(x cos a - y sin a + e )^ + (a ^/b (x sin a + y cos a+e)^=a^ (9) 

g g X'' ' c ' c g ■'g y c 


where a^. and b^, are the semi-axes of the ellipse defining the grain ex- 
terior. The burning perimeter at each reference plane is determined as 
follows. Equations (8) and (9) are rewritten in polar coordinates r and 
0. The perimeter is determined through numerical integration of the 
equation . 



( 10 ) 


0 


At each value of 0 in the integration the radial coordinate of the grain 
bore and exterior is computed from (8) and (9) , respectively. If the 
exterior value exceeds the bore value at the 0 position, r is given the 
value calculated for the bore; otherwise, a zero value is assigned. 


Burning perimeters for each reference plane thus vary from one SRM to 
another as a result of eccentricity of the grain bore and exterior, which 
is specified by the independent variables ex and ey, orientation of the 
ovality of the bore with respect to the exterior, and the out-of-'roundnesses 
of both the bore and exterior. It is possible that some correlation may 
exist between these variables in particular cases. However, if the case 
is segmented, it would appear that the head end and aft end geometric 
features are Independent and they are treated as such in the remainder of 
this report. On the other hand, in general, there should be a close cor- 
respondence between the exterior grain geometries and the locations of 
the bore centers at the aft ends of Ljg and z^, and it has been assumed 
that these are identical. The geometry of the interior grain is determined 
by the distributed variables Di, Zj^, and x-pa; l.e., these determine the 
initial semi-axes and the exterior geometry is determined by D© and z^ where 
Zc has been introduced to account for axial variations in the outside dia- 
meter of the grain. It is clear that the independent specifications of ex 
and ey at both ends of the rocket provides for a statistically distributed 
malalignment of the grain bore and exterior. 


From the burning perimeter values obtained as outlined above, correction 
factors are next calculated and applied to the standard calculations of 
perimeter to account for the ovality and malalignment. For this purpose 
when solution of the ovality equations indicates that bumthroughs have 
occurred at adjacent reference planes, the burning perimeters are assumed 
to vary linearly between the planes. When bumthrough has occurred at only 
one of two adjacent reference planes, the portion of the correction factor 
applicable to each end is weighted in proportion to the corresponding 
length that has or has not experienced bumthrough to determine the overall 
correction factors (computer symbols; SEN, SEH, CHINAV, CHIN and CHIH) . 
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It is recognized that the approach represents considerable idealiza- 
tion of the general behavior to be expected, particularly with respect to 
the assumed qualities of the axial distributions of the parameters between 
reference planes. Nevertheless, the model seems to capture the essence 
of the performance effects associated with ovality and malalignment. Even 
with highly nonlinearly distributed ovality and eccentricities between 
reference planes the effects on performance should be roughly the same as 
determined here. Also, factors such as the precise shape of the ovals 
would appear to play a secondary role in influencing performance variations. 

Care has been taken in the above discussion to differentiate between 
the grain exterior and the case Interior. As far as the equations go the 
pertinent item is the grain exterior. Its shape, however, may be Influenced 
by the case and it is a choice of the user as to whether or not the statis- 
tical variations of ovality of the case alone should determine the variation 
in the grain exterior or if variation in Insulation thickness should be 
statistically combined with the case variation to arrive at the qualities of 
the grain exterior. 

As mentioned earlier, the analysis may also be applied to star grain 
configurations. In this case, the star points are in effect disregarded. 

The rationale for this is that by far the most important effects of malalign- 
ment and ovality occur just before and during tailoff and are dominated by 
the behavior of the remaining propellant web which may be approximated for 
this purpose by a circular perforated grain of the same web thickness. The 
capability to treat star grains has thus been incorporated into the computer 
program. However, it is applicable only when the entire grain is a star 
grain. When both circular perforated and star grains are present together, 
it is assumed that the circular perforated grains determine the performance 
characteristics in as far as the effects of ovality and malalignment are 
concerned. 

An additional application for the ovality and malalignment is note- 
worthy. This has to do with the temperature gradient within the propellant 
grain which is an important variable that has not been taken directly into 
account in the analysis because no convenient way has yet been foxind to do 
so. Some insight into the gross effects of such gradients can be gleaned 
by statistically incorporating additional variations in the out-of-roundness 
of the grain to represent the general effects of having a biaxial burning 
rate which could be different in magnitude and orientation at the two main 
reference planes. 

Statistical Analysis 

This section describes the statistical analysis used in selecting values 
for the distributed variables. The basic computational methods as related 
to the logic of the statistical procedure used in the computer program are 
presented. 

The computer routine is designed to accept several different types of 
data, perform the specified operations required to obtain a frequency dis- 
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tribution for each variable, and select a value based on Its statistical 
frequency curve. The routine is divided into two subroutines. One sub- 
routine, "setup," is called once and generates the required frequency curves 
for the statistical variables given. The other, "input," is called as re- 
quired and provides specific values for each variable given. The latter 
subroutine utilizes the Monte Carlo technique for selection of values which 
will be discussed later. Figure II-2 is a general flow chart of the setup 
and input subroutines. 

Frequency curves . The primary task of the setup subroutine is to 
obtain a frequency curve for each statistical variable from the data given. 
For each variable, the ultimate product of this subroutine, a cumulative 
distribution function (CDF), is obtained from its frequency curve. A 
CDF is a step function which jumps at regular intervals and is constant 
between jump points. At each jump point the magnitude of the jump is the 
probability that the variable will be within that interval; thus, as the name 
implies, the probabilities are accumulated over the range of the statistical 
variable. Examples of the types of CDFs produced by the setup subroutine are 
shown in Figure II-3. This CDF can be obtained in the present program from 
several types of input data, ranging from raw data points to specifying the 
actual CDF directly (See Section IV) . 

Basically, there are two classes of requests allowed for input statis- 
tical variables. The first class contains all variables which require little 
if any statistical analysis. This includes such requests as to hold certain 
variables constant or to obtain the CDF directly from a given histogram, as 
illustrated in Fig. II-3A. The other class requires statistical analysis to 
obtain the frequency curve. The analysis consists of obtaining the first 
four moments of the statistical variable and generating an equation which 
approximates the actual frequency curve for the data. The method selected 
for obtaining the frequency curve from the first four moments is known as 
Pearson's system (Ref. 7). The variables and sequence of calculations were 
chosen to parallel those of Pearson's system which is discussed in sufficient 
detail in Ref. 7 to permit direct adaptation to a computer program. 

Basically, Pearson's system consists of a family of curves with a cri- 
terion vsilue used to determine which equation of the family best describes 
the data. This criterion value, K, is evaluated as follows: 

K = 3^(32 + 3)2/4(432 - 333^)(232~ 33^^ - 6 ) (11) 


where 


Here ^ 2 * 1 ^ 3 » 3*^® the second, third and fourth moments of the variable 

about its mean, respectively. 
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The curve is determined to be one of the three main types as follows; 


K < 0 

0 < K < 1 

K < 1 

TYPE I 

TYPE IV 

TYPE VI 


Transitional types exist, but are not incorporated in the present program. 

The statistical analysis section of the subroutine initially obtains 
the first four moments. It then evaluates K and determines the appropriate 
type of curve to use. From these moments and the total number of data 
points for each statistical variable, the appropriate parameters for the 
proper equation are evaluated. Thus the equation representing the fre- 
quency curve is established. This equation is evaluated over the specified 
range and using Simpson's rule to Integrate the curve, the CBF is obtained 
for each statistical variable, as shown in Fig. II-6B. The CDF for each 
variable is then stored in an array. 

Monte Carlo . When called, the Monte Carlo or input subroutine generates 
values for all statistical variables. These values are determined based on 
the frequency curve for each variable. The Monte Carlo subroutine performs 
this function by obtaining a random number and using the CDF array produced 
earlier. 

The CDF consists of values from 0. to 1. with corresponding variable (x) 
values over the appropriate range. A random number (RANDU) with a value 
between 0. and 1. is obtained from the random number generator discussed 
below. This value establishes the appropriate CDF value, and the correspond- 
ing X value is then selected from the array to fix the proper value of the 
variable. Each time the Monte Carlo routine is called, it determines new 
values for all the statistical variables; that is, it selects a new set of 
variables for each SRM whose performance is to be analyzed. This is accom- 
plished as follows. 

In the computer program, the CDF is stored in an array of 100 elements 
(a different number of elements may be specified by program changes). 

A second array (x) also contains 100 values corresponding to the possible 
values which the statistical variable may assume. For any x array value, the 
corresponding CDF array value is the accumulated area under the frequency 
curve to that point, expressed as a percentage of the total area: 

x(N) x(lOO) 

f(x)dx/ ^f(x)dx (16) 

x(l) x(l) 

where f(x) is the frequency distribution. The values of the CDF array, 
thus, will progressively Increase from an initial value of zero to a final 
value of one at a rate depending on the given frequency curve. To randomly 


CDF(N) 


I 
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select a value for a given statistical variable, the random number generator 
(RANDU) is invoked to obtain a random number between zero and one. This 
random number is compared against CDF array values until the smallest CDF 
value that is larger than the random number is found. The value of the 
corresponding element of the x array is assigned as the appropriate random 
value for the given statistical variable. The random numbers, obtained from 
a rectangular distribution from zero to one, are thus transformed by the CDF 
in such a way that as a series of random numbers is selected for a single 
variable, the corresponding variable values occur at a frequency correspond- 
ing to their probability distributions. 

Random number generator . Based on its size, five statements, this sub- 
routine would not appear to be worthy of a lengthy discussion. The subject, 
however, is quite critical to the statistical portion of this program. With- 
out some understanding of the random number generator (RNG) the results may 
prove both surprising and less than satisfactory. 

The numbers generated by this routine are not actually random numbers, 
but rather are referred to as pseudo-random numbers. This fact may cause 
a less than comfortable feeling. It turns out, however, that the character- 
istics of pseudo-random numbers are adequate for the present purposes. Real 
random numbers do not cycle. The pseudo-random numbers have a cycle, but 
the period is reasonably large. For IBM 370 hardware configurations, approx- 
imately 2^*^ selections are made. It is unlikely that this program will be 
used extensively enough for this to become a problem. The major consideration 
when using the RNG is that it requires a number to initiate the random gen- 
eration. After the initial "seed" number, the routine will automatically 
reset this number for the next iteration as is done when a new set of variables 
is to be selected for an SRM. For any given seed number, the random number 
produced is always the same and the new seed nunber generated will also be 
equivalent. Thus, if one hundred random numbers are generated from a given 
initial seed number, the same, exact set of numbers will be generated given 
that same initial seed number. Often this is a useful characteristic; how- 
ever, it is essential to luiderstand that to obtain two different sets of 
random numbers, two different initial seed values must be given. 


III. DISCUSSION OF INPUT AND OUTPUT 

In this section, each Input parameter is defined in the order in 
which it is encountered in the program. The English or Greek symbol is 
given first followed by the computer symbol in parentheses. The English 
or Greek symbols are provided mainly for convenience in consulting the basic 
analyses of References 3, 4 and 5 as only a few of these are used in the 
present report. Where appropriate, additional discussion of the variable 
and recommended or typical nominal numerical values are given. Sketches 
of geometric characteristics are presented for clarification. Although 
many of the variables are the same as used in References 3, 4 and 5, the 
discussion relative to these is repeated here for ready reference. In 
addition, the present outputs of the program are defined. It is the aim 
of this section to provide a guide to the user in the specification of his 
problem and the Interpretation of the outputs. 

Concise printed definitions of all input variables are also given 
with the computer program listing (See Section IV) . The definitions 
appear in groups throughout the "main" program and the "area" subroutines. 

In general, each group is divided in the computer program and in this 
section into two subgroups: the first group containing the variables de- 

scribing the SRM which are always fixed for the Monte Carlo Investigation 
and the second group containing those that are subject to being statistic- 
ally distributed. 

The original basis for classification of the two types of variables 
was to be the analysis of the sensitivities of the performance of the SRM 
to the variable under consideration. Of course, each variable considered 
would also have to be an independent one or at least relatively so. Each 
variable whose variation could reasonably be expected to influence the 
performance calculations for the SRM was examined one at a time. A few 
variables whose effects were obviously of very minor significance were 
omitted at this point and classified without further evaluation as non- 
statlstical (fixed value) variables. Among these were the temperature and 
pressure sensitivity coefficients of the propellant properties. The 
thrust loss coefficients were also classified as fixed, as there is no 
practical way of distinguishing in experimental data analysis between the 
variations in thrust due to the variation in the statistical variables 
and that due to other factors. A baseline design was then selected based 
on the nominal values for the sample case SRM of Section V and the 
performance of a pair of SRMs computed. One motor in the pair had the 
baseline values of the variable and the other had the baseline values 
except for the variable under consideration which was changed a small amount. 
The amount of the change was somewhat arbitrary but represents a rough 
estimate of the maximum range of variation for a wide variety of SRMs 
scaled to the size of the baseline SRM; i.e., the range reflects an esti- 
mate of the maximum variation to be expected when no special attention is 
given to control of the variables during manufacture. In the case of 
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the examination of the effects of grain ovality and malalignment, several 
variable changes were examined together in order to reduce computational 
time . 

Partial results are given in the Appendix in the form of the pairs 
of thrust-time traces obtained for each variable. These plots should be 
of Interest to the SRM developer as they indicate possible results of 
deviations from manufacturing specifications and tolerances. Their 
utility as far as the present study goes is somewhat limited as it was 
decided that for the Monte Carlo investigation all but a very few of the 
variables would be statistically determined in order to minimize cumulative 
errors that might result from neglect of a large number of seemingly 
unimportant variables. It will be noted that there are no plots in the 
^pendlx for a few of the distributed variables listed in the computer 
program and in this section of the report. This is because either the 
variable has no influence for the design or option selected (e.g., 0j.n 
and 0j.jj for Cgp = 1) ; the effects of the variable is an obvious one (6) ; 
data is not presently available for a meaningful investigation (a and 0) ; 
or the final selection of variables differed from that considered during 
the earlier sensitivity, investigation (e.g,, was substituted for y 

and C* and %g and were both Introduced as independent variables 
eliminating Ts) . 

The present section provides information only on nominal values and 
the units of the variables. Instructions for preparing the input for the 
description of the distribution of statistical (distributed) variables is 
given in Section IV based on the analysis of Section II. It should be 
noted here that the statistical variables may be held constant if desired 
by use of the proper code on the data cards for the variables. With minor 
program modifications it is also possible to accommodate additional dis- 
tributed variables, the number of such variables being limited only by 
the core storage capacity of the computer. The listing of variables 
follows. 

Seed Number 

I^^ (IXl) The seed number for the random number generator. An odd 

eight-digit integer should be used. This number initiates 
the generation of random numbers forming the basis for 
selection of the various sets of input variables for each 
SRM. The seed number must be changed to change the sets 
of variables generated when repeating an analysis of the 
same motor pairs (see Section II for further discussion) . 

User's Options - Fixed Values 

I^^(IEO) 0 For no consideration of grain ovality or malalignment. 

1 For consideration of above. Calculations of ovality 

and malalignment effects approximately triples computer 
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time requirements. These calculations should be 
by-passed by use of the option provided if there is a 
basis for assuming in a particular case that the 
effects are negligible. 


0 

For no plots and no statistical or difference 
analysis of results. 


1 

For plots and tabular output of motor pair differences. 


2 

For tabular output of motor pair differences without 
plots . 


3 

For plots and statistical analysis of results only. 
Codes 1 and 2 will also yield statistical analysis 
of the motor pair results. 


N(j) (NUMPLT(J) An Integer designating whether or not a specific output 
plot is desired: 

0 If a specific plot is desired. 


1 If a specific plot is not desired. 

The order of specification of NUMPLT(J) is as follows: 

1 F versus T for the motor pair. 

2 F versus T for the motor pair during talloff. 

3 Difference in F between the pair of motors versus T. 

A Difference in ITOT between the pair of motors versus T. 

5 Time Integral of absolute value of difference in F 

of motors versus T. 


Propellant Characteristics - Distributed Values 

p (RHO) Density of the solid propellant (slugs/in^) 

P 

aj^(Al) ,32(A2) Propellant burning rate coefficients in the equation 
r = aP*^ above and below the transition pressure, re- 
spectively (in/sec-psia’^) . 


n^(Nl) ,n2(N2) Burning rate exponents corresponding to a^ and a 2 » 
respectively (1) . 

a , (ALPHA) Erosive burning coefficient in the Robillard-Lenoir burning 

^ rate rule (Equation III-ll, Ref. 3) (in^* ®-f t*^ • ®/sec^ • ® Ibf®*®). 
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B (BETA) 
Koaa(ROal) 


Basic Motor 
L (L) 

(TAU) 


Basic Motor 
Dg (DE) 
(DTI) 
(THETA) 
(ALFAN) 


Erosive burning pressure coefficient in the Robillard- 
Lenolr rule (1). 

The oxidizer to aluminum ratio (1). Variations, in this 
quantity determine variation in the the nno chemically 
determined characteristic exhaust velocity CSTABN and the 
ratio of specific heats GAMN of the chamber gases. The 
latter variations are determined by a regression analysis 
which gives CSTARN and GAMN as functions of ROAL. The 
resultant functional relationship must be used in the 
program where indicated on the program listing (Section IV) . 
CSTARN and GAMN are determined at 1000 psi chamber pressure 
and 60“F propellant bulk temperature. If no aluminum is 
present in the propellant CSTARN and GAMN should be set at 
this point to their nominal values at the reference con- 
ditions. 

Dimensions - Fixed Values 


Total initial perforated grain length including gaps at 
slots (in.). This is used only in the erosive burning rate 
equation. An estimate will suffice. 

Estimated web thickness of the controlling propellant 
length (in.). (See ZO below for definition of control- 
ling length.) The actual web thickness is calculated from DI 
and DO (or RC and RSI, RIWS or RIWl'7)and includes statistical 
variations. If the grain is tapered, the length average 
value should be specified excluding lengths having addition- 
al taper not represented by ZO and segments located any- 
where having a step decrease in thickness. Such step de- 
crease must be handled by the additive tabular input option 
if they introduce two significantly different web thick- 
nesses for the same grain type; e.g., two circular perforated 
grains. If a circular perforated grain is used, it is 
assumed that it will have the approximate average web thick- 
ness of the controlling propellant length. 

Dimensions - Distributed Values 


Diameter of the nozzle exit (in.). 

Initial diameter of the nozzle throat (in.). 

Cant angle of the nozzle with respect to the motor (degrees) . 
Exit half angle of the nozzle (degrees) . 
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Isj,^ (LTAP) 


*Ta 

Zq (ZO) 


z (ZC) 
c 


Estimated axial length of grain at the aft end at start 
of first tailoff having an additional taper not repre- 
sented by ZO or THETAG (in.). See Figure III-l. This 
variable permits the designer to specify an additional 
taper at the noz 2 le end of a circular perforated or star 
grain to produce regresslvlty shortly before tailoff. 

Difference in web thickness at ends of LTAP (in.). 

See Figure III-l. 

Initial difference between web thicknesses due to grain 
bore taper at the head and aft ends of controlling grain 
length, not including any initial length associated with 
LTAP or THETAG (in.). The controlling length of the grain 
is the axial distance between the head end of the grain 
and the position of expected maximum Mach number in the 
port. In general this length terminates whenever there 
is an abrupt decrease in web thickness near the aft end 
of the chamber. 

Initial difference between web thicknesses due to grain 
exterior taper at the head and aft ends of the controlling 
grain length (in.). 


Ovality and Malalignment - Distributed Values 

(Not required if lEO = 1). The following variables characterize the 
ovality and lack of concentricity of the grain interior and exterior at 
two reference radial cross-sections - one near the head end of the grain 
and one near the nozzle end. Based on the mathematical analysis of the 
burning perimeters at these two planes, correction factors are calculated 
and applied to the burning surface calculations to account for approximate 
effects of ovality and malalignment. The effects of interior and exterior 
grain taper are taken into account through the parameters ZO, ZC, and XT. 
See Section II for details. 


AR (RONDCN) , 
cn 

AR , (RDNDCH) 
cn 


One half the difference between the maximum and minimum 
diameter of the grain exterior at the nozzle and head end 
reference planes, respectively (in.). 


AR (RONDGN) , 
gn 

AR , (RONDGH) 
gh 


One half the difference between the maximum and minimum 
diameter of the grain interior at the nozzle and head end 
reference planes, respectively (in.). 


e (EXN) 
xn 

e (EYN) 
yn 


The eccentricity of the center of the grain interior with 
respect to the center of the grain exterior in the x^, and 
directions, respectively (See Fig. II-l) at the nozzle 
end reference plane (in.). 




Figure III-l. Basic motor dimensions. i-Gci* ^Gni» used only for circular perforated grains 

The controlling grain length may be or some lesser value depending on the position 

of maximum port Mach number, is the length average web thickness over the grain 
length excluding the lengths associated with Lj.^, 0 q and the head end dome. 
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e^(EXH) , 

eyh(EYH) 


The eccentricity of the center of the grain interior with 
respect to the center of the grain exterior in the 
and directions, respectively, at the head end reference 
plane (in.). 


o^^(ALPHAN) , The angular orientation of the ovality of the grain 

interior with respect to the grain exterior at the nozzle 
a^(ALPEAH) and head end reference planes, respectively (degrees). 


Basic Performance Constants - Fixed Values 


Ay^^ (DELTAY) Initial desired bum Increment (in.). This Increment will 
be used by the program for initial 5% of the web thickness 
burned, for the period shortly preceding and following 
talloff and at such other times as the rate of pressure 
change is large. Larger increments will automatically be 
used at other times. An initial increment size of approx- 
imately O.OOIt^, is recommended for the Monte Carlo program. 
If the Increment size used is twice the recommended value, 
the maximum thrust imbalance calculated may be decreased 
by as much as 5 percent, representing a loss in program 
accuracy. 


Ij (II) Number of steps in the integrations of perimeters in the 

OVAL subroutine. Approximately 25 steps appears to offer 
a good compromise between accuracy and computation time, 
the latter being more strongly Influenced than the former 
by the step size. 

x^^^ (XOUT) Distance burned at which the propellant breaks up (in.). 

This permits the option of specifying termination of 
burning resulting from possible structural breakup of 
propellant. If this option is not desired, XOUT should 
be set to some large value; e.g., 1000 in. 


(AP/Ay)out 

(DPOUT) 


Cj. (ZETAF) 


Rate of change of pressure with respect to distance burned 
at which the propellant is extinguished (Ib/in^). This 
permits the option of specifying termination of burning 
when it is determined that an abrupt tailoff will not 
permit the computer to handle the rapid change in surface 
area. If a gradual tailoff is expected, DPOUT may be set 
to some large value; e.g., 10,000psia/in. For large 
motors (120 in. dia. and up) where the tailoff is expected 
to be abrupt, a value of 500 psla/in. is recommended. In 
general, larger values may be used for smaller motors. 

Thrust loss coefficient (1). In the absence of data to 
the contrary, a value of 0.98 is recommended. 
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t 

max 


q(TMAXQ) 


Estimated time at which maximum dynamic pressure on the 
vehicle occurs (sec.)* This permits the program to com- 
pute an estimate of the thrust differential of the motor 
pair at the time of maximum dynamic pressure. 


t^i(TB) 


Estimated burning time (sec.). This and HB below permit 
the program to calculate delivered specific impulse and 
thrust based on an assumed trajectory which was de- 
termined from analysis of typical large SRM applications. 


(HB) Estimated burnout altitude (ft.). To obtain sea level 

performance characteristics, HB should be set equal to 
zero. 

P^^^(PREF) Reference average nozzle stagnation pressure used in the 

nozzle throat erosion equation (psia) . See also ERREF 
below. 


^^^(DTREF) Reference nozzle throat diameter used in the nozzle throat 
erosion equation (in.). See also ERREF below. 

(it )j^(PIPK) Temperature sensitivity coefficient of pressure at con- 
" stant ratio K of burning surface to throat area (/®F) . 


a ^^(CSTART) Temperature sensitivity of CSTAR at constant K (/®F). (See 
^ Eq. 4, Section II). A typical value is 0.000038. 


a ^ (CSTAEP) Pressure sensitivity of CSTAR (l)(See Eq. 5, Section II). 
c p 

Ptran^^^^^^ Transition pressure at which the burning rate coefficient 
and exponent change (psia) . 


a (GAMP) 
YP 


Pressure sensitivity of the ratio of specific heats (1) 
(See Eq. 6, Section II). A typical value is 0.00527. 


Basic Performance Constants - Distributed Values 


E^^^ (ERREF) 


T (TGR) 

T. (TIGR) 
xgr 


Reference nozzle throat erosion rate (in/sec) (See Eq. 1 
Section II.) A set of typical values for a carbon tape- 
phenolic impregnated throat is 0.00763 in/sec, 560 psia 
and 57.285 in. for ERREF, PREF and DTREF, respectively. 

Bulk temperature of the propellant grain ("F). 

Ingition delay (sec.) at a TGR of 60°F. The time required 
to reach ninety-five percent of initial equilibrium pressure 
should be used. 


The Program and Basic Grain Configuration and Arrangement - Fixed Values 
I^p (INPUT) 1 For only tabular input. 
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G „ (GKAIN) 
op 


S (STAR) 
op 


NOTE: 

(NT) 

Q (ORDER) 
op 


2 For only equation input. 

3 For a combination of 1 and 2. 

1 For an entirely circular perforated grain. 

2 For star grain only. 

3 For a combination of 1 and 2. 

0 For an entirely circular perforated grain. 

1 For standard star (See Fig. III-2) . 

0 

2 For truncated star (See Fig. III-3) . 

3 For wagon-wheel (See Fig. III-4) . 

The different types of star grains may not be combined in a 
single configuration. 

Number of thrust termination passageways in grain (1) 

NT is zero if there are no thrust termination passage- 
ways . 

1 If a star grain is at head end and a circular 
perforated grain at aft controlling end. 

2 If a circular perforated grain is at both ends. 

A star grain segment may still be present. 

3 If a circular perforated grain is at head end and 
a star at the aft controlling end. 

4 If a star grain is at both ends. 

If grain = 1, value of order must be 2. 

If grain = 2, value of order must be 4. 

It is Important to realize that ORDER establishes the 
controlling port area equations to be used. Thus if the 
nozzle end segment is not indeed the controlling one 
(the one that establishes the maximum Mach number in the 
port) , ORDER should be specified to designate the actual 
controlling segment as the nozzle end segment. GRAIN, 

STAR, NT and ORDER are not used for INPUT =1, but values 
must be assigned for continuity of computer operations. 






(final) bu; 









Figure III-4. Wagon-wheel grain cross-section. Calculated variables are circled 
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C (COP) For extreme ends of a circular perforated grain only: 

op 

0 If both ends are conical or flat. 

1 If head end is conical or flat and aft end is spherical. 

2 If both ends are spherical. 

3 If head end is spherical and aft end is conical or flat. 
Tabular Burning Surface and Port Areas - Fixed Values 


\pi 

\sT 

Vt 

Vt 




Burning area in the port (in^) . 

Burning area in the slots (in^) . 

Burning area at the nozzle end (in^) . 

Port area at the nozzle end of the controlling grain 
length (in^) . 

Port area at the head end of the grain (in^) . APHK and 
APNK are not required when INPUT is 2 or 3; the equation 
Inputs must be used to provide the information in these 
cases. Values of all A's must be specified to completely 
describe the burning surface and port areas versus distance 
burned. The program computes intermediate values by 
linear interpolation. The number of values required is 
arbitrary and limited only by the storage capacity of 
the computer, but values must be specified for y = 0. 

Also, as an example of the procedure for specifying 
terminal values which must be followed, burning surfaces 
must be specified as zero at burnout of the tabular 
surfaces and at a y exceeding the highest anticipated 
calculated performance value. Separate input cards must 
be prepared for each value of y and arranged as described 
in Section IV of this report. The use of these tabular 
values in conjunction with the equation inputs (INPUT = 3) 
increases the flexibility of the program considerably. 
Frequently it is easy to estimate burning surface effects 
which end constraints on the equation inputs neglect. In 
this case a table of input values can be readily prepared 
from estimates of the effects. Also, the program can be 
run first without the tabular values and outputs used as 
an aid in obtaining the estimate. For example, the burning 
perimeters of a star grain can be determined in this way 
and the values used to estimate the effects of a head end 
closure on the star grain. 

Initial volume of chamber gases associated with tabular 
input (in^). 



Geometry for 
(XTZO) 

S (S) 

Geometry for 
Dq (DO) 

\ (DI) 

0- (THETAG) 

Loci (LGCI) 
L^l (L®I) 
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Circular Perforated Grain - Fixed Values 


Difference between the initial circular perforated grain 
diameter at the nozzle end of LGCI and the nominal value 
of DI less ZO and less twice XT (in.) (See Fig. III-l) . 

Number of burning flat ends of a circular perforated 
grain not including an extreme aft grain end (1) . 

Circular Perforated Grain - Distributed Values 


Length average outside diameter of circular perforated 
grain, excluding lengths extending into the closure (in.). 

Length average inside initial diameter of circular per- 
forated grain (in.). Only the controlling length exclud- 
ing LTAP should be considered in the averaging. 

Angle burning surface element of circular perforated 
grain located at the extreme nozzle end of chamber makes 
with the motor axis (degrees). See Figure III-l. THETAG 
must be set to zero if a star grain is located at the 
nozzle end (GRAIN = 3, ORDER = 3) or if aft end burning 
surfaces are represented by tabular values. THETAG is 90“ 
if the circular perforated grain represented by equations 
has a flat burning surface located at the extreme nozzle 
end of the chamber. If THETAG is less than or equal to 5“, 
a value must be assigned (zero is satisfactory) , but the 
effect of THETAG on burning surface area is not computed 
(See also LGNI and LTAP) . THETAG is zero if the end sur- 
face is flat and inhibited. 

! 

Initial total axial length of circular perforated grain 
represented by equation inputs not including gaps (in.). 
LGCI excludes lengths associated with THETAG. 

Initial slant length of a burning conical circular per- 
forated grain at the nozzle end (in.). LGNI is set 
equal to zero if THETAG is less than or equal to 5“. In 
this case the length otherwise associated with LGNI should 
be added to LGCI. If the error in burning surface area 
thus introduced is deemed significant, a correction may 
be introduced by making use of tabular inputs in combi- 
nation with the equation inputs. Basic effects of small 
THETAG on tailoff may be accounted for by specification 
of LTAP. If a nozzle end burning surface is flat (THETAG 
= 90“) LGNI equals one half the difference between inside 
and outside local grain diameters. 
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6 (THETCN) The contraction angle of a circular perforated grain 

bonded to the nozzle closure (degrees). See Figure III-l. 
Use an estimated value which yields approximately the 
correct volume of propellant burned. If a star shaped 
grain is located at the extreme nozzle end of the chamber 
or if tabular values are used to represent downstream 
burning surfaces, THETCN is 90“ . THETCH is also 90“ if 
the extreme aft end of the grain is inhibited, but only a 
flat-ended, inhibited grain (THETAG = 0) which does not 
extend into the nozzle closure may be accurately represented. 
THETCN is zero for a burning flat end (THETAG = 90“) which 
does not extend into the closure. THETCN must be assigned 
a value even if COP is 1 or 2, but it will not affect the 
numerical results. 

6 ^^ (THETCH) The contraction angle of a circular perforated grain bonded 
to the head end (degrees). See Figure III-l. Use an 
estimated value which yields approximately the correct 
volume of propellant burned. THETCH is 90“ if the extreme 
forward end of the circular perforated grain (bonded or not) 
represented by equations is flat. A head end flat burning 
surface is treated by proper specification of S. THETCH 
must be assigned a value even if COP is 2 or 3, but it will 
not affect the numerical results. 

Basic Geometry for Star Grains - Fixed Values (The wagon-wheel is considered 
a type of star grain for the purpose of this program.) 

Number of burning flat end surfaces of a star grain not 
located at extreme nozzle end of the chamber (1) . 

The number of star points (1). 

Number of burning flat end surfaces (0 or 1) of a star 
grain located at the extreme nozzle end of the chamber (1) . 


Hg (NS) 

Up (NP) 

n (NN) 
n 


Basic Geometry for Star Grains - Distributed Values 


LGsi(LGSI) 


\ (RC) 
f (FILL) 


Initial total axial length of star-shaped perforated 
grain represented by equations (in.). No provision 
comparable to the use of LGNI for circular perforated 
grains is made here to treat effects resulting from 
THETAG greater than 5“. Adjustments may be made, however, 
by use of tabular input values in conjunction with the 
equation Inputs. Also, effects of taper, including addi- 
tional small taper at the nozzle end, on tailoff may be 
treated by use of the variables ZO, XT and LTAP. 

The star grain outside radius (in.). 

The fillet radius at star valleys (in.). See Figure 
III-2, 3, A. 
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Special Geometry for Wagon-Wheel Grain - Distributed Values 


The length average initial radius of the inside of the pro- 
pellant web (in.). 

The lengths of the pairs of parallel sides of the first and 
second sets of grain points, respectively (in.). See Figure 
III-4) . 

Th® angles between the slant sides and the center lines of the 
CATPHAi ATPHA 9 ^ the first and second sets of grain points, respec- 

' ' ^ tlvely (degrees). The angles should not exceed 90 degrees. 

h (HW) The half-width of the star points (in.). HW must not exceed 

” TAUWW (HW < TAUWW). 


Riww(RIWW) 

Zl,i2 (L1.L2) 


Special Geometry for Truncated Star Grain - Distributed Values 


The length average initial radius of the truncation (in.). 
See Figure III-3. 

The length average initial radius of the Inside of the 
propellant web at the bottom of the slots (in.). 

Special Geometry for Standard Star Grain - Distributed Values 

0 (THETAF) Angular location of the fillet center of standard star from 

^ the line of symmetry (degrees) . 

The apex angle of the star point (degrees) . 

The initial radius of the inside of the propellant web of 
standard star grain (in.). See Figure III-2. If the 
grain is tapered, the length average value should be used. 

Geometry of Thrust Termination Passageways - Distributed Values 


0 (THETAP) 
P 

R (RIWS) 
Iws 


Rp <RP) 
Rls (RIS) 


(LTP) 



(DTP) 


0^p (THETTP) 


Teff (TAUEFF) 


Initial length of the termination passageway between the 
centers of gravity of perimeters of the bases (in.). See 
Figure III-l. 

Initial diameter of the termination passage (in.). 

The acute angle between the axis of the passage and the motor 
axis (degrees) . 

Estimated effective web thickness at the termination port (in.) 
The user must judge the distance burned at which the effect 
of the termination passage on modification of the burning 
surface geometry ceases to be significant. In general, this 
should be between two-thirds and full web thickness. The 
equation used to account for the burning surface is based on 
a passageway in a circular perforated grain terminated at the 
case by a flat inclined plane. Thus only a rough estimate 
of the effect of the termination passage is provided. 
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Speclal Equation Inputs -(Fixed Relationships required only at the option 
of the user. May be used when INPUT * 1, 2, or 3.) 




Additive burning surface input as function of y for port 
burning surface (in^). 


»bs (BBS) 


Additive burning surface input as function of y for slot 
burning surface (in^). 


“b„ (BBN) 


Additive burning surface input as function of y for 
nozzle end burning surface (in^). In order to make use 
of the option of specifying the B's, a minor program 
modification is required. The B's are all set equal to 
zero in the present program. If this option is to be used, 
the program statements assigning values to the B's are 
easily replaced with the desired equation inputs. 


Program Outputs 

The variables whose values are printed by the present program are 
defined below. Additional variables may also be printed with minor program 
modifications. In addition to the variables listed below, the present pro- 
gram prints out values of all input variables including those selected for 
each SRM from statistical distributions. The input characteristics of the 
statistical distributions are also printed. 

Time dependent data - single motors 


t (T) 


Operating time (sec.). This is calculated from the time 
of initiation of ignition. 


y (Y) 

P (PONOZ) 
on 

P, (PHEAD) 
h 

A. + A, + A. 
Dp Ds on 

(SUMAB) 

F (F) 

I^ (ITOT) 


Average distance burned (in.). 

Stagnation pressure at the nozzle end of the chamber 
(Ib/in^). 

Pressure at the head end of the chamber (Ib/ln^). 

The total burning surface of the propellant (in^) . 


The delivered thrust based on the assumed trajectory (Ibf) . 
Losses are included. 

Total delivered specific impulse (Ibf-sec) . Losses are 
included. 
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Time Independent data - single motors 


(WPl) 


Propellant weight calculated from mass discharge rates 
(Ibm) . 


W 2 (WP2) Propellant weight calculated from the products of burning 

" surfaces and Incremental distances burned (Ibm) . 

W (WP) Arithmetic average of WPl and WP2 (Ibm) . A check on the 

^ calculation accuracy Is provided by comparison of this 

with WPl and WP2. 


Maximum head end chamber pressure calculated by the pro- 
gram (Ib/in^). 

Initial seed number for the current configuration. New 
seed numbers are automatically selected for each variable 
based upon this seed number. The initial seed number may 
be used as an input seed number to reproduce the calcula- 
tions for the configuration. See Section II. 

Final seed number for the current configuration. This final 
seed number may be used as an input seed number to continue 
the calculations without random number cycling before 2^*^ 
random numbers have been selected. 

Time dependent data - motor pairs . The following data are available 
in tabular and/or graphical form, subject to the option of the user (See 
IPO) . To compute the imbalances between two motors the data for the motor 
which has the fewer computational y-steps is subtracted from the data of 
the other motor to determine the imbalances. No difficulty with regard to 
the interpretation of the results arises from this since in most instances 
it is the absolute value of the differences which is important. For those 
variables whose absolute values are not taken, changes in sign of the dif- 
ferences indicate points where the crossings between the two motors' traces 
occur. 


P, (PHMAX) 
h max 


hi <™> 


(IX) 


AF (FDIFF) 


The difference in thrust between the two motors (Ibf) . 


AI^ (IDIFF) 


The difference in total impulse between the two motors 
(Ibf-sec) . 


Iai^I(iadiff) 


The absolute value of the difference in total impulse 
between the two motors (Ibf-sec). 



( 17 ) 
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Tlme independent data - motor pairs 


(FMAX1,TFMX1, The maximum and minimum algebraic thrust Imbalance (Ibf) 

FMIN1,TFMN1) and the time (secs) at which they occur, respectively, 

during web action time. 


(FMAX2 , TFMX2 , The maximum and minimum algebraic thrust imbalance (Ibf) 

FMIN2,TFMN2) and the times (secs) at which they occur, respectively, 
during talloff. 


(TDFT01,TDFT02, The time at which talloff begins for the first motor of 
DTW) a pair to begin talloff, for the final motor to begin 

talloff and the absolute value of the difference between 
the two times, respectively, (secs). 


(FW1,FW2,DFW) The thrust at the beginning of tailoff for the first 

motor of a pair to begin tailoff, for the final motor to 
begin tailoff and the absolute value of the difference 
between the two thrusts, respectively, (Ibf). 


(DFT01,DFT02) The absolute value of the thrust imbalance when the first 
motor of a pair begins tailoff and the final motor begins 
tailoff, respectively. 


(DFMQjTMAXQ) The absolute value of the thrust Imbalance (Ibf) which 
exists when the maximum dynamic pressure occurs on the 
vehicle and the estimated time (secs) at which this event 
occurs, respectively. 

(AFMAXjTFMAX, The absolute value of the maximum thrust imbalance (Ibf) 
AFMAXT , TFMAXT) which exists and the time (secs) at which it occurs during 
web action time and talloff, respectively. 


(FDIFIG,TD1FIG) The absolute value of the thrust Imbalance (Ibf) which occurs 
during the Initial portion of operation (t < 0.02 t^,) and 
the time (secs) at which it occurs, respectively. 


(DIT.ADIT) 


The total Impulse imbalance and the absolute value of the 
total Impulse imbalance (Ibf-secs) accumulated during tailoff, 
respectively: 



(18) 



( 19 ) 


where tj.j^(TDFT01) is the earlier time at which tailoff 
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(DFIOOK.TIOOK) 


begins In the two motors and t^, is the time at which 
operation of both motors ends. 

The absolute value of the thrust imbalance (Ibf) which 
exists when the last motor reaches 100,000 lb. thrust 
during tailoff and the time (secs) at which it occurs, 
respectively. 



IV. THE COMPUTER PROGRAM 

This section contains the instructions for the preparation and arrange- 
ment of the data cards. Also, a complete listing of the program statements 
is given. The program was written for use on an IBM 370/155 computer and 
requires approximately 168K storage locations on that machine. The program 
also is designed to be used with a CALCOMP 663 drum plotter. The plotter 
requires one external storage device (magnetic tape or disk) . In addition 
to the one storage device required for the plotter, three other 
external storage units are required. Unit 1 is used to store the output 
data, pertinent to the imbalance calculations , for the first motor in each 
pair of motors. Unit 2 is used to store the nonstatlstlcal data which remain 
constant for all of the motors. Unit 4 is used to store the values of the 
statistical variables for use with each motor. Only minor program modifica- 
tions are required to eliminate the plotting capability of the program. 

Also, Unit 2 can be eliminated by using repeated sets of data cards for the 
nonstatlstlcal variables. Hence, it is relatively simple to modify the pro- 
gram to require only 2 external storage units. Elimination of the other 
two external storage units would require significant program modification. 

Input Data 


The discussion below gives the general purpose, order and FORTRAN 
coding information for the input data. 

Card 1 Total nimiber of individual motors to be analyzed (42X,I2) 

Col. 1-42 NUMBER OF CONFIGURATIONS TO BE TESTED = 

43-44 Number of rocket motors to be analyzed 
Card 2 Initial seed number (IIP) 

Col. 1-10 Initial 8-10 digit seed number 

It is necessary to describe one type of statistical analysis for each 
statistical input variable. The method for doing this is described below 
using Cards 3 through 9. Note that only one type of statistical analysis 
may be requested for each variable. Hence, only the card or cards necessary 
for that particular type of statistical analysis are input for each variable. 
For example, to obtain a Type II analysis only Card 5 and Cards 5A would be 
would be used. In addition, it is necessary that the data cards for the 
variables to be used in a given configuration be placed in the order in 
which they are input into the computer program. In some cases certain 
variables are not required for an analysis. In such cases, the cards for 
those variables should be omitted. As many Cards 5 through Cards 9A as 
required may be used. 
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Card 3 Variable name (2A4)(one card for each variable) 


Col. 1-8 Name of statistical variable. 

NOTE; One Card 3 Immediately precedes the Card 4 thru Card 9B 
used for each variable. Also, END should be used as the 
last variable name before using Card 9B below. 

Card 4 Input for Type I statistical analysis (12, 2X, 7E10.0) 


Col. 1-2 


5-14 

15-24 


Code = 10 


Code = 11 


Raw data given; obtain CDF directly from 
histogram. 

Raw data given; obtain CDF from Pearson's 
equation of the frequency curve. 


XI = Number of raw data points given. 


X2 = Mean value of first interval of histogram. 


25-34 X3 = Histogram interval width. 

35-44 X4 = Number of intervals in histogram. 

45-74 Blank 

Card 4A Subsequent Type I data cards (10E8.0 ) 


.Col. 1-8 
9-16 

. t 

72-80 


Raw data points equivalent to the number 
specified in XI. Ten data points per card for 
as many cards as required (e.g., 46 data points 
would require 5 data cards with the last card 
having the final four fields blank) . 


Card 5 Data input for Type II statistical analysis (12, 2X, 7E10.0) 


Col. 

1-2 i 

^ Code = 20 

I 

Histogram given; obtain CDF directly from 
histogram. 



Code = 21 

Histogram given; obtain CDF from Pearson's 
equation of the frequency curve. 


5-14 XI = Number of Intervals in histogram. 

15-24 X2 = Mean value of first interval of histogram. 

25-34 X3 = Interval width. 


35-74 


Blank 
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Col. 


Col. 


Col. 


Col. 


Card 5A Subsequent Type II data cards (10E8.0) 


1-8 

9-16 


The same number of data 
points as specified in XI, 
for as many data cards as 
necessary 


72-80 

Card 6 Input for Type III statistical analysis (12, 2X, 7E10.0) 

= 31 Four moments given; obtain CDF from Pearson's 

equation of the frequency curve. 

First moment about zero. 


1-2 

Co< 

5-14 

XI 

15-24 

X2 

25-34 

X3 

35-44 

X4 

45-54 

X5 

55-64 

X6 

65-74 

X7 


Fourth moment about mean. 


NOTE: No data cards required. 

Card 7 Input for Type IV statistical analysis 


1-2 

Code = 40 CDF given; read in 

the given 

5-14 

XI = Number of Intervals in CDF. 


15-24 

< 

^ - Mean value of first interval 

of CDF. 

25-34 

X3 = Interval width. 


35-74 

Blank 


Card 7A 

Subsequent Type IV data cards (10E8.0) 


1-8 CDF values corresponding to the cumulative 

9-16 frequency up through each Interval. Data 

I should be provided for as many Intervals as 

indicated by the vadue given for XI. 


72-80 
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Card 8 Input for Type V statistical malysis (Use appropriate 
card below) 

Card 8A Normal distribution to obtain CDF. 


Col. 1-2 


Code = 51 


5-14 XI = Mean of normal distribution. 

15-24 X2 = Standard deviation. 

25-34 X3 = Beginning X value of CDF (optional) . 

35-44 X4 = Ending X value of CDF (optional) . 

45-74 Blank 

NOTE: If either X3 or X4 is omitted, a three-sigma limit is assumed; 

thus, if both values are left blank, a six-sigma limit will be 
generated by the program. If a zero value is desired for X3 
or X4, ±.0000001 should be used instead. 

Card 8B Rectangular distribution to obtain CDF (12, 2X, 7E10.0) 


Col. 1-2 

Code =52 

5-14 

XI = Beginning X value. 

15-24 

X2 = Ending X value 

25-74 

Blank 

Card 8C 

J-Distribution to obtain CDF 

Col. 1-2 

Code = 53 

5-14 

XI = Mean (beginning X value) 

15-24 

X2 = Standard deviation. 
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Col. 


Col. 


Col. 


Col. 


Col. 


Card 8C (Cont'd) 

25-34 X3 = Ending X value (optional) 

35-74 Blank 

NOTE: The J-dlstrlbutlon Is defined herein as the right half of a 

normal frequency curve. The Xl value specified should be the 
mean as If the full normal curve were being specified. The 
X3 value is optional; if not specified, a three sigma limit 
will be assumed. If zero is desired for the X3 value, 
+.0000001 should be used instead. 

Card 9 Input for Type VI statistical analysis (use appropriate 
card below) 


Card 9 A Use a constant for this value (12, 2X, 7E10.0) 


1-2 

Code = 60 

Use a constant value for this variable. 

5-14 

XI = Desired 

constant value. 

15-74 

Blank 


Card 9B 

Indicates end 

of data (12) 


1-2 Code = 90 

Card 10 Initialization of variables (22F3.1) 

1-66 Zero's or blank card 

Card 11 Ovality and output options (5X, II, 5X, II, 9X, 51) 



7-11 IPO 




0 


12 




No ovality analysis 
Ovality analysis 

No plots or statistical analysis 
Plots, statistical analysis and tabular output 
Tabular output and statistical analysis 
Plots and statistical analysis 



-40- 


Col. 


Col. 


Col. 


Card 11 (Cont'd) 


13-17 NUMPLT(J) 

Plot thrust time trace 

18 

1 

Do not plot thrust time trace 
Plot tailoff thrust time trace 
Do not plot tailoff thrust time trace 
Plot thrust imbalance 
Do not plot thrust i^alance 
Plot Impulse imbalance 
Do not plot impulse imbalance 
Plot absolute Impulse imbalance 
Do not plot absolute Impulse imbalance 
Card 12 Nonstatistlcal motor dimensions (3X, F10.2, 5X, F10.3) 


19 


20 


21 


22 


C 

C 

c 

c 

c 


1-3 L = 

4-13 Value of L 

12-18 TAU = 

19-28 Value of TAU 

Card 13 Nonstatistlcal performance constants (requires 4 data cards) 


Card 13A (8X, FlO.3, 4X. 14. 6X. F10.2. 7X, F10.2. 7X. F10.4 ) 

1-8 . DELTAY = 

8-18 Value of DELTAY 

19-22 II = 

23-26 Value of II 

27-32 XOUT = 


33-42 


Value of XOUT 
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Col. 


Col. 


Col. 


Card 13A (Cont'd) . . ; . ; 

43-49 DPOUT = 

50-59 Value of DPOUT 

60-66 ZETAF = 

67- 76 Value of ZETAF 

Card 13B (4X. FlO.l. 4X, FlO.l, 6X. F10.2, 7X. F10.3. 6X. F10.5) 
1-4 TB = 

5-14 Value of TB 
15-18 HB = 

19-28 Value of HB 

29-34 PREF = 

35-44 Value of PREF 

45-51 DTREF = 

52-61 Value of DTREF 

62-67 PIPK = 

68- 78 Value of PIPK 

Card 13C (8X, FlO.7, 7X, F10.2. 8X, F10.7. 6X. F10.7) 


1-8 

CSTART = 


9-18 

Value of 

CSTART 

19-25 

PTRAN = 


26-35 

Value of 

PTRAN 

36-43 

CSTAEP = 


44-53 

Value of 

CSTARP 

54-59 

GAMP = 



60-69 


Value of GAMP 



Col. 


Col. 


Col. 
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Card 13D (7X, F10.3) 

1-7 TMAXQ = 

8-17 Value of TMAXQ 

Card 14 Description of type of grain configuration (9X, 12. 9X, 12 . 

8X, 12, 6X, F4.0, 9X, 12. 7X. 12) 

1-9 INPUT = 

10-11 Value of INPUT (1, 2 or 3) 

12-20 GRAIN = 

21-22 Value of GRAIN (1, 2 or 3) 

23-30 STAR = 

31-32 Value of STAR (0, 1, 2 or 3) 

33-38 NT = 

39-42 Value of NT 

43-51 ORDER = 

52-53 . Value of ORDER (1, 2, 3 or 4) 

54-60 COP = 

61-62 Value of COP (0, 1, 2 or 3) 

Card 15 Tabular values for geometry at y = 0.0 (requires 2 data cards ) 
(Not required if INPUT = 2) 

Card 15A (6X, F6.2, IPX, F11.2. IPX. F11.2. 8X, Fll.2) 


1-6 

YT = 

7-12 

0.0 

13-22 

ABPK = 

23-33 

Value of ABPK 

34-43 

ABSK = 


44-55 


Value of ABSK 
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Col. 


Col. 


Col. 


Col. 


Card 15A (Cont'd) 

55-62 ABNK = 

63-71 Value of ABNK 

Card 15B (22X, F11.2. 9X. F11.2. 8X. F11.2) 

1-22 APHK = 

23-33 Value of APHK 
34-42 APNK = 

43-53 Value of APNK 

54-61 VCIT = 

62-72 Value of VCIT 

Card 16 Tabular inputs for y greater than 0.0 (requires 2 data cards 
for each y value) (Not required for INPUT = 2) 

Card 16A (6X. F6.2. IPX. F11.2. IPX. F11.2. 8X. F11.2) 


1-6 

YT = 


7-12 

Value of 

YT 

13-22 

ABPK = 


23-33 

Value of 

ABPK 

34-43 

ABSK = 


44-54 

Value of 

ABSK 

55-62 

ABNK = 


63-73 

Value of 

ABNK 

Card 16B 

(22X, F11.2. S 

1-22 

APHK - 


23-33 

Value of 

APHK 


34-42 APNK = 


43-53 Value of APNK 
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Col. 


Card 17 Non-statistlcal c.p. grain geometry (Not required for GRAIN 
= 4)(6X. F10.3, 3X. FIO.O) 

1-6 XTZO = 

7-16 Value of XTZO 

17-19 S = 

20-29 Value of S 

Card 18 Non-statistlcal star grain geometry (Not required for 
GRAIN = 2)(4X. FIO.O, 4X, FIO.O. 4X. FIO.O) 


1-4 

NS => 



5-14 

Value 

of 

NS 

15-18 

NP = 



19-28 

Value 

of 

NP 

29-32 

NN = 



33-42 

Value 

of 

NN 


Finally, Figure IV-1 is a schematic representation of the data deck 
construction, and Table IV-1 presents an example set of data. This is the 
same data as used in the sample case presented in Section V of this report. 
Note that these are all the data cards which are required for this example 
for any number of configurations. 


Program Listing 

Figure IV-2 shows a block diagram of the overall program and Table 
IV-2 presents the complete program listing. 

As previously mentioned, the program has been designed to produce 
graphical presentations of the computational results. Program statements 
that must be removed in order to delete the plotter compilation require- 
ments are indicated by check marks (►) in Table IV-2. Removal of these 
statements is necessary if the user's computer is not equipped for CALCOMP 
plotting. However, if other plotters are available, generally only the 
plotting subroutines need be replaced. 
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Figure IV-1. Schematic of data deck. 





TABLE IV-1. EXAMPLE DATA SHEETS 





TABLE IV-1. EXAMPLE DATA SHEETS (CONT'D) 
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TABLE IV-1. EXAMPLE DATA SHEETS (CONT'D) 











TABLE IV-1. EXAMPLE DATA SHEETS (CONT'D) 





SETUP INPUT OVAL AREAS PAIR SIGBAR OUTPUT PLOTTER 















TABLE IV-2 


C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 




♦ MCNTF CARLO PERFORMANCE ANALYSIS OF SRM PAIRS * 

♦ PREPAREO AT AUHURN UNIVERSITY ♦ 

♦ UNDER MUD. NO. 12 TO CCOPERATIVE AGRFEMFNT WITH ♦ 

♦ NASA MARSHALL SPACE FLIGHT CENTER * 

♦ * 

♦ BY ♦ 

♦ R. H. SFORZINI, W. A. FOSTER, JR. AND J. S. JOHNSON, JR. ♦ 

♦ AEROSPACE ENGINEERING DEPARTMENT ♦ 

♦ NOVEMBER 197A ♦ 


INTEGER GRAIN 
REAL IDIFF 

REAL MGENtMDI S.MNOZ ,MNl , JROCK,N,L,MEl ,ME, ISP, I TOT , MU, MASS , I SPVAC 

REAL N1 ,N2,NSEG,K1 ,K2 ,KCH,KEN,NS,LCC,LTAP 

REAL M2,MCRAR, ISP2, I T VAC , KA , KH , LAMBDA 

COMMO.N/CONSTl / ZW, AE , AT,THE TA, ALFAN 

COMMO.N/CUNST2/CAPGAM,HE ,BUr ,ZE TAF , TB ,HB , GAM 

C 0 M M 0 N / C 0 N S T B / S , N S , G R A I N 

COMMON/CONSrA/CELDI , DC , D I , ZC , X T , ZO 

C0MMCN/C0NST5/KPLT 

COMMQN/VAR I Al/T,nELY,DELTAT,PCNOZ,PllEAD,RNOZ ,RHEAD, SUMAB,PHMAX 

COMMON/ VAR I A2/AB PORT, ABS LOT, ABNOZ,APHE AD, APNOZ,DADY ,ABP2,ABN2,ABS2 

COMMON/VARIA3/ITnr, I T VAC , JROCK , I SP , I SP V AC , f'D I S , MNCZ , SG, SUMMT 

COMMON/VARI AA/RNT,RHI ,SUM2,Rl ,R2,R3,RHAVE ,RNAVE,RBAR,YH,KCUNT 

COMMCN/VARI AS/ABMAIN, ABTO,SUMCY,VCI , VC 

COMMON/VAR I A6/YDI , TE 

COMMON/ VAR I A 7/ Y, THRLST 

COMMGN/PLGTT/n'n.'iDt'M, I P T , I OP 

COMMCN/PLCT2/NUMPL I 

CGMMGN/0/ALA/ChIH,CHlN, SEN,SEH,CFINH, AZ,BZ 
C0MMCN/0VALB/CHINN,CH INAV,S£NN 

COMMGN/OVALC/KUNOCN,RONDCH,RUMjGN,RONnGH,EXN,EYN,eXH,EYH, 

2ALPHAN, ALPHAH 
COMMON/QVALM/Z ,ZQ 
COMMON/OVALM2/KK I , I I 
COMMON/SEEC/ IX 

COMMCN/PA IRl/Twl , T W2 , D T W , F « I , Fw2 , DF W 1 , DF W2 , DEW , TMAX C , DFMQ , 
2FDIFF,TDIFF,NX 

COMMCN/PAIR2/FMAXI , TFMXl ,FMIN1 , TFMNl, 

2 FMAX2 , TFMX2 ,FHI N2 , FFMiJZ 

COMMC.N/PA IR3/AFMAX, TFMAX ,AFMAXT , IFMAXT 
COMMON/UUT 1/FDIF IG, TD IF IG,DI r ,ADI T 
C0MMCN/UUT2/CF100K, TIOOK 
COMMON/DA TA2/I 

COMMON/TOFF/DFTOl ,DFT02, TDFTOl ,TDFT02 
DIMENSIOi FCIFF( AGO ) ,TDIFF(ACO) 

DIMENSION NUMPLT(5) 
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TABLE IV-2 (CONT'D) 


DATA PI ,G/3. IA159»32. 1725/ 

READ(5,500) NRUMS 

C ♦ READ IN THE NUMBER OF CONFIGURATIONS TO BE TESTED ♦ 

500 FORMAT I A2X, 12 ) 

NPAIRS=NRUNS/2 

I0P = 0 

NPL0T=0 

KPL0T=0 

TW1=0.0 

FW1=0.0 

WRITEI6, 11112) 

11112 F0RMAT(20X, 'DATA FOR STATISTICAL ANALYSIS PROGRAM*) 

CALL SETUP 
DO 901 I=1»NRUNS 
REWIND 2 
IX1=IX 
CALL INPUT 
WRITE(6,A02) I 

602 FORMAT! 1H1,A2X, ‘CONFIGURATION NUMBER *,I2) 

IF(I-l) 5CCC, 5000, 5001 

5000 READ (5, 499) SUMDY , ANS , ZW , Y , T , DELTAT , RNOZ , RHEAD, SUMAB, PHMAX, SUM2, I T 
10T,RHT,RNT ,R1 ,R2 ,R3,RHAVE,RNAVE,RBAR, I TVACtSUMMT 

WRITE (2, 49 9) SUMD Y , A NS , Z W , Y , T , DEL T AT , RNOZ , RHEAD , SUM AB , PHMAX , SUM2 , I 
1T0T,RHT,RNT,R1,R?,«3,RHAVE,RNAVC,RBAR, ITVAC.SUMMT 
GO TO 5002 

5001 RE AD (2, 49 9) SUMDY , ANS , 7 W , Y , T , DEL TAT , RNOZ , RHEAD , SUMAB , PHMAX, SUM2 , IT 
lOT, KH T, RN T,R1,R 2, R 3, RHA VE, KNA VE ,RBAR, I TVAC, SUMMT 

5002 CONTINUE 

C ♦ SET INITIAL VALUES OF SELECTED VARIABLES EQUAL TO ZERO ♦ 

C ♦ ♦♦♦NOTE’^'** THESE VALUES MUST BE ZEROED AT THE BEGINNING OF ♦ 

C ♦ EACH CCNFIGURATION RUN * 

Q ti:)::t‘***4ii‘«‘***‘******************<‘*********^********4‘****$*******4l****** 

499 F0RMATI22F3.1 ) 

IF(I-l) 5003,5003,5004 

5003 READ(5,491) I EC , 1 PO , ( NUMPLT I J P ) , JP= 1 , 5 ) 

WRITE (2, 491) IEO,IPO, ( NUMPLT UP ) , JP= 1 , 5 ) 

GO TO 5005 

5004 READ(2, 491) I EC , I PO , ( NUMPLT (J P ) , JP= I , 5 ) 

5005 CONTINUE 

491 F0RMAT(5X,I1,5X,I1,9X,5II) 

Q t4^4‘***<‘****4‘*********4‘*:^*****1‘’^******************************* ******* 


C * READ IN THE USER'S OPTIONS ♦ 
C ♦ ♦ 
C ♦ VALUES FOR lEO ARE ♦ 
C ♦ 0 FOR NO OVALITY ♦ 
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TABLE IV-2 (CONI' D) 


C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


♦ I FOR OVALITY ANALYSIS 

♦ VALUES FOR IPG ARE 

♦ 0 FOR NO PLOTS AND NO STATISTICAL ANALYSIS 

♦ I FOR PLOTS AN!) TABULAR OUTPUT 

♦ 2 FOR TABULAR OUTPUT ONLY 

♦ 3 FOR PLOTS ONLY 

CONTINUE 

♦ VALUES FOR NUMPLT(J) ARE (NOT REQUIRED FOR IPO=0,3) 

♦ 0 IF SPECIFIC PLOT IS DESIRED 

♦ I IF SPECIFIC PLOT IS NOT DESIRED 

♦ ORDER OF SPECIFICATION OF NUPPLT(J) IS 

♦ I THRUST VS TIPE (ENTIRE TRACE) 

♦ 2 THRUST VS TIME (TAILCFF PORTION ONLY) 

» 3 THRUST imbalance VS TIME 

« 4 TOTAL IMPULSE IMBALANCE VS TIME 

♦ 5 ABSOLUTE TOTAL IMPULSE IMBALANCE VS TIME 


♦ 

♦ 

♦ 

* 

* 

♦ 

* 

« 

♦ 

* 

♦ 

♦ 




WRITE (6,492) I£0, IPQ, ( NUMPLT ( JP ) , JP= 1 , 5 ) 

492 FORMAT!//, 2CX, ’OPTIONS*, /,13X,‘ IEO= • , 1 1 , / , 1 3X , • I P0= *,11, 
2/,l3X, *NUMPLT(J)= *,512) 
mil F0RMATIE16.9) 

READ(4, 1 U ll ) RHO,Al ,Nl, A2,N2,ALPHA,BETA,R0AL 


C 

C 

C 

C 

C 

C 

C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 




♦ READ IN BASIC PROPELLANT CHARACTERISTICS 

♦ ♦ 

♦ THE FCiLLGWING VARIABLES ARE OBTAINED FROM THE STATISTICAL ♦ 

♦ ANALYSIS PROGRAM ♦ 

♦ ♦ 

♦ RHO IS THE DENSITY OF THE PROPELLANT IN LBM/IN«'*3 » 

♦ A1 IS THE BURNING RATE COEFFICIENT BELOW THE TRANSITION ♦ 

♦ PRESSURE * 

♦ N1 IS THE BURNING RATE EXPONENT BELOW THE TRANSITION PRESSURE * 

♦ A2 IS THE BURNING RATE COEFFICIENT ABOVE THE TRANSITION ♦ 

♦ PRESSURE ♦ 

♦ N2 IS THE BURNING RATE EXPCNENT ABOVE THE TRANSITION PRESSURE « 

♦ ALPHA AND BETA ARE THE CONSTANTS IN THE EROSIVE BURNING ♦ 

♦ RELATION OF ROBILLARD AND LENIOR ♦ 

♦ ROAL IS THE OXIDIZER TO ALUMINUM RATIO ♦ 

♦ ♦ 

♦ DEFINE CSTARN AND GAMN ♦ 

♦ * 

♦ CSTARN is the nominal THERMOChEM ICAL CHARACTERISTIC EXHAUST ♦ 

♦ VELOCITY IN FT/SEC AT IGOO PSI AND 60 OEG F * 

♦ GAMN IS THE NOMINAL RATIO OF SPECIFIC HEATS FOR THE ♦ 
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TABLE IV-2 (CONT'D) 


C ♦ PROPELLANT GASES ♦ 

Q ************************** ****** ^**«****^ 

c * ♦ 

CSTARN=-I7.84 75<'ROAL + 52 39.7 
GAMN=ROAL*5.67357E-3+ 1.11707 

C ♦ ♦ 

Q ************************************** ******^f********* ******** ******* 

WRITE (6, 60 3) RHO, A 1 , N I , A2 , N2 , AL PHA , BE T A , ROAL ,CSTARN ,GAMN 
603 FORMAT! // ,20X ♦* PROPELLANT CHARAC TER I S T ICS * , / , 1 3X, • RHO= »,F8.6,/,1 
23X,'A1= • ,F7.5,/, 13X, *N1= • , F*i. 3 , /, 13X , • A2= * , F 7. 5 , / , 1 3X , • N2= *, 
3F5.3,/»13X, *ALPHA= » ♦ F A . I , / » 1 3X , • BE T A= • . F5 . I , / * 1 3X , • ROAL= *»F7.A 
4*/, 1 3X, 'CSTARN= • , 1 PE 1 1 . A ♦ / ♦ 1 3X ♦ 'GAMN= '♦IPEll.A) 

IF! IPO)4002, 4002,3999 

► 3999 IF!I.EQ.U CALL GS I 2£ ! 1 2C0. 0 , 1 1 .0, 1 12 1 1 
IFl!-!)*^!) 4000,4000,4001 

4000 REWIND I 
KPLT=l 

GO TO 4002 

4001 KPLT=2 

4002 CONTINUE 
RHO=RHO/G 

IFII-1) 5006,5006,5007 

5006 READ!5,502) L , TAU 
WRITE!2,502) L,TAU 
GO TO 5008 

5007 READ!2,502) L,TAU 
5C08 CONTINUE 

IF!IEO) 6CCC,6C00,600l 

6000 READI4, I 11 1 I ) CE , DT I , THE T A, AL FAN, L T AP, X T , ZO , ZC 
GO TU 6002 

6001 REA0I4,mil) CE,DT I , THETA, ALFAN,LTAP,XT, ZO,ZC, 
2RONDCN,R()NDCH,RONOGN,RONOGH,EXN,EYN,EXH,EYH,ALPHAN, ALPHAH 

6C02 CONTINUE 


c ♦ READ IN BASIC MOTOR DIMENSIONS ♦ 

C ♦ ♦ 

C ♦ L IS THE TOTAL LENGTH OF THE GRAIN IN INCHES ♦ 

C ♦ TAU IS THE ESTIMATED AVERAGE WEB THICKNESS OF THE CONTROLLING ♦ 

C ♦ GRAIN LENGTH IN INCHES * 

C * ♦ 

Q ********************************************************************* 

C * THE FOLLOWING VARIABLES ARE OBTAINED FROM THE STATISTICAL ♦ 

C ♦ ANALYSIS PROGRAM * 

Q ********************************************************************* 

C ♦ ♦ 

C ♦ DE IS THE DIAMETER OF THE NOZZLE EXIT IN INCHES ♦ 

C ♦ DTI IS THE INITIAL DIAMETER OF THE NOZZLE THROAT IN INCHES * 

C ♦ THETA IS THE CANT ANGLE OF THE NOZZLE WITH RESPECT TO THE ♦ 
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TABLE IV-2 (CONI' D) 


C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


♦ MOTOR AXIS IN DEGREES * 

♦ ALFAN IS THE EXIT HALF ANGLE OF THE NOZZLE IN DEGREES ♦ 

♦ LTAP IS THE LENGTH OF THE GRAIN AT THE NOZZLE END HAVING ♦ 

♦ ADDITIONAL TAPER NOT REPRESENTED BY ZO IN IN’CHES * 

♦ XT IS THE DIFFERENCE IN WEB THICKNESS ASSOCIATED WITH LTAP ♦ 

♦ ZO IS THE INITIAL DIFFERENCE BETWEEN WFH THICKNESSES IN * 

« INCHES DUE TO GRAIN BORE TAPER AT THE HEAD AND AFT ENOS ♦ 

♦ OF THE CONTROLLING GRAIN LENGTH ♦ 

» ZC IS THE INITIAL DIFFERENCE BETWEEN WEB THICKNESSES IN * 

♦ INCHES DUE TU GRAIN EXTERIOR TAPER AT THE HEAD AND AFT ♦ 

♦ ENDS OF THE CONTROLLING GRAIN LENGTH ♦ 

CONTINUE 

♦ RUNDCN AND RCNDCH ARE ONE HALF THE DIFFERENCE IN INCHES ♦ 

♦ BETWEEN THE MAXIMUM AND MINIMUM DIAMETER OF THE GRAIN * 

♦ EXTERIOR AT THE NUZZLE AND HEAD END REFERENCE PLANES * 

♦ RESPECTIVELY ♦ 

♦ RONDGN AND RONDGH ARE ONE HALF THE DIFFERENCE IN INCHES ♦ 

0 BETWEEN THE MAXIMUM AND MINIMUM DIAMETER OF THE GRAIN ♦ 

♦ INTERIOR AT THE NOZZLE AMD HEAD END REFERENCE PLANES ♦ 

♦ RESPECTIVELY ♦ 

« EXN.EYN.EXH AND EYH ARE THE ECCENTRICITIES IN INCHES OF THE ♦ 

♦ CENTER OF THE GRAIN INTERIOR Vil TH RESPECT TU THE GRAIN ♦ 

♦ EXTERIOR AT THE NUZZLE AND HEAD END REFERENCE PLANES ♦ 

♦ RESPECTIVELY ♦ 

♦ ALPHAN AND ALPHAH ARE THE ANGULAR ORIENTATIONS IN DEGREES ♦ 

♦ OF THE ovality OF THE GRAIJ INTERIOR WITH RESPECT TO ♦ 

♦ THE GRAIN EXTERIOR AT THE NOZZLE AND HEAD END REFERENCE * 

♦ PLANES RESPECTIVELY ♦ 


50Z FORMAT13X,FIO,?,5X,F10.3) 

IFIIEU) ftC03,6C03,600A 

6003 WRI TE (6,60A0) L , T AU , OE , D T I , THETA , AL FAN , L T A P , X T , ZO , ZC 
GO TO 6005 

600A WRITE(6,604) L , T AU , l)E , D T I , THE TA , ALF AN » L T AP » X T » ZO , ZC » 

2R0NDCN,RaNCCH,R0NDGN,RUNDGH,EXN,EYN,EXH,£YH, ALPHAN, ALPHAH 
6005 CONTINUE 

604 F0RMATI//,20X, 'BASIC MOTOR D I MENS I ONS • , / , 1 3X , • L= » , F 8 . 2 , / , I 3X , 
I'TAU^ • ,F6.3,/, 13X, 'DE= 

21PE 11.4,/ , 13X, 'DTI = • , IPEl 1.4, / , 13X, • THE TA= ' , 1 PE I I . 4 , / , I 3 X , • ALFAi^^ 
3 • ,IPE11.4,/,13X,'LTAP= ' , I PE 1 1 . 4 , / , 1 3X , ' X T = • , 1 P E I I . 4 , / , I 3X, • Z0 = 
4* , lPEll.4 ,/, 1 3X, ' ZC= *, 

51P611.4,/, 13X, •RfJNOCN= • , IP L I I . 4 , / , 13X , 'RONDCH= • , I PE 1 1 . 4 , / , 1 3X , 
6*R0NDGN= • , lPEll.4,/, 13X, 'RONDGH= • , IPE I I . 4 , / , 1 3X , • EXN= ',1PE11.4, 
7/,13X,'EYN= • , IPEU.4,/,13X, 'eXH= • , IPE 1 1 .4 , / , I 3X , • £YH= ',IPEU.4, 
8/, 13X, • ALPHAN= • , I PE 1 1 . 4 , / , I 3 X , • ALPHAH= ',1PE11.4) 

6040 FORMAT (// ,20X, • BAS IC MOTOR U I ME NS IONS ' , / , I 3X , • L= • , F8 . 2 , / , 1 3X , 
l'TAU= • ,F6.3,/,13X, 'CE= ', 

21PE I 1.4, /, 13X, ‘DTI = • , IPE I 1.4, /, 13X, • THETA= ',1PE11.4,/,13X, • ALFAN = 


- 55 - 



TABLE IV-2. (CONI' D) 


C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


3 • , IPtl l.A,/, 13X, *LTAP= •,IP£11.4,/,13X,*XT= • , I PE 1 1 . 4 , / , 1 3X, » ZO= 
4* ,IPE11.4,/,13X, *ZC= ',IPE11.4) 

THETA = THFTA/‘i7.2 95 78 
ALFA.N = ALFAN'/5 7.29 57 3 
ALPHAN = ALPhA;M/b7.29 5 78 
ALPHAH=ALPHAH/ 57.29578 
IF(I-n 5CC9,5C09,5010 

5009 REA0(5,50 3) CELT AY, I I t XOUT , DPGUT , ZETAF , T8,HB, PREF ,OTREF ,P IPK, 
2CSTART,PTRAN,CSTARP,GAMP,TMAXQ 

WRITE (2,503) CELT AY, I I , XfJUT , OPGUT , ZET AF ,TB, HB, PREF , OTREF , P I PK , 
2CSTAKT,PTRA.N,CSTARP,GAMP,TMAXQ 
GO TO 5011 

5010 RE AD (2, 50 3) CELT AY, 1 1 , XQUl , DPGUT , ZETAF , TB,HB , PREF ,DTREF , P I PK , 
2CSTART,PTRAN,CSTARP ,GAMP,TMAXQ 

5011 CONTINUE 

READI4,lllll) ERREF,TGR,TIGR 


* READ IN BASIC PERFORMANCE CONSTANTS ♦ 

♦ ♦ 

♦ DELTAY IS THE DESIRED BURN INCREMENT DURING TAILOFF IN INCHES ♦ 

♦ II IS THE NUMBER OF INTEGRATION STEPS USED IN OVAL « 

=0" XOUT IS THE DISTANCE BURNED IN INCHES AT WHICH THE PROPELLANT ♦ 

♦ BREAKS UP ♦ 

♦ OPOUT IS THE DEPRESSURIZATION RATE IN LB/IN**3 AT WHICH THE ♦ 

♦ PROPELLANT IS EXTINGUISHED ♦ 

♦ ZETAF IS THE THRUST LOSS COEFFICIENT ♦ 

♦ TMAXQ IS THE ESTIMATED TIME AT WHICH THE MAXIMUM DYNAMIC ♦ 

♦ PRESSURE OCCURS ON THE VEHICLE ♦ 

♦ TB IS THE ESTIMATED BURN TIME IN SECONDS ♦ 

♦ HB IS THE ESTIMATED BURNOUT ALTITUDE IN FEET ♦ 

♦ PREF IS THE REFERENCE NOZZLE STAGNAIICN PRESSURE ♦ 

♦ DTREF IS THE REFERENCE THROAT DIAMETER ♦ 

♦ PIPK IS THE TEMPERATURE SENSITIVITY COEFFICIENT OF PRESSURE ♦ 

♦ PER DEGREE F AT CONSTANT K ♦ 

♦ CSTART IS THE TEMPERATURE SENSITIVITY PER DEGREE F OF CSTAR ♦ 

♦ AT CONSTANT PRESSURE ♦ 

♦ CSTARP IS THE PRESSURE SENSITIVITY OF CSTAR ♦ 

♦ PTRAN IS THE PRESSURE IN PSIA ABOVE WHICH THE BURNING RATE ♦ 

♦ EXPONENT CHANGES * 

♦ GAMP IS THE PRESSURE SENSITIVITY OF GAM ♦ 

*^** **^ *!* ********* ’tf*’^**** ^ 4 ^ ************ 4 ********* 

CONTINUE 

♦ THE FOLLOWING VARIABLES ARE OBTAINED FROM THE STATISTICAL ♦ 

♦ ANALYSIS PROGRAM ♦ 

♦ * 

* ERREF IS THE REFERENCE THROAT EROSION RATE IN IN/SEC ♦ 
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TABLE IV-2 (CONI' D) 


C * TGR IS THE BULK TE^!PERA^URE OF THE GRAIN IN DEGREES F ♦ 

C ♦ TIGR IS THE IGNITION DELAY IN SECONDS AT 60 DEGREES F ♦ 

503 FORr'ATIBXtFlO.atAX. U,6X,F10.2,7X,F10.2.7X,F10.A,/,^«X,F10.1,AX, 
2F10.1,6X,F10.2,7X,F l0.3,6X,F10.5t/,8X,Fl0.7t7X,Fl0.2t8X,Fl0. 7, 
36X,F10.7»/,7X,F10.3 ) 

WRI TEI6,606) DELTAY, I I,X0UT,DP0UT,ZETAF,TR,H8,ERREF,PREF,DTREF 

2, TGR,PIPK,CSTART,PTRAN.CSTARP,TIGR,GARP,TMAXQ 

606 F0RMAT(//,20X, 'BASIC PERFORMANCE CONSTANTS I 3X , *DELTAY= *,F5.3, 
l/,13X, • I 1= • , lA, 

1/, 13X, 'XOUT= • »F7.2,/, 13X, »DPOUT= • , F9 . 2 , / , 1 3X , • ZET AF= ',F6.A,/,13 
2X,'TC= • ,F5. I,/, 13X, 'HB= • ♦ F 7 . 0 , / , 1 3X , • ERRE F= • 

3, F8.5,/tl3X,»PREF= * , F8 . 2 , / , I 3X , ' DTREF = * , F 7 . 3 , / , 1 3X , ' TGR= ',F7.3, 
A/, 13X, *PIPK= • ,F7.5,/, 13X, •CSTART= • , F 1 0 . 7 , / , 1 3 X , • P TR AN= *tF8.2 
5,/, 13X, 'CSTARP= • ,F10.7,/, 13X, »TIG= • , F 7. A , / , 1 3X , • GAMP= »,F10.7, 
6/,13X, 'TMAXC= ',F7.3) 

TIG=TIGR^EXPIPIPK*(60.0-TGR) ) 

T=TIG 

A=Al 

N=Nl 

C STARR=C ST AKN*EXP(C START* ( TGR-60. ) ) 

GAM=GAMN 

Q=A*EXP(PIPK*( l.-M)*(TGR-60.I ) 

KKI=0 

CHIH=1.0 

CHIN=1.0 

CHIN :N=l.O 

CHINAY=l.O 

SEN=0,0 

SENN=0.0 

SEH=0.0 

CHINH=l.O 

NDUM=0 

IPT = 0 

MN1=.85 

M£l=7.0 

Z = Z(J+ZC 

ZO=ZO 

S=0.0 

NS=0.0 

K0UNT=0 

KEWAT=0 

ABMAIN=0.0 

ABTO=0.0 

TW2=0.0 

UTW=0.0 

FW2=0.0 

DFW=0.0 
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TABLE IV-2. (CONI' D) 


D6LY=DELrAY 

TOP=GAM+l. 

BOT=GA«-l. 

2AP = TOP/(2.«'BOT) 

CAPGAN=SORT(GAM)*(2./TOP»**ZAP 

AE=PI*OE«DE/^. 

1 IF(XT.LE.O.O) TE=0.0 

TCALL= ( TAU-XT-ABS( Z/2. I)/1.05 
IFnEO.EQ.l.AND.Y.GT.TCALL ) CALL OVAL 
IFCXT.LE.0.0) GO TO AO 
TL=lY-TAU+XT+Z/2. )«LTAP/XT 
IF(TL.LE.O.O) TL=0.0 
IF(TL.GE.LTAP) TL=LTAP 
TE=LTAP-LTAP*CHINAV 
IF(IEO.EO.O) TE=TL 
AO IF(r-TIG) A1,A1,A2 
Al OT=DTI 

CSTAR=CSTARR 
GO TO A3 

A2 RADER = ERREF<'( ( PONQZ/PREF )**0.8 ) ♦( COTREF/DT ) ♦♦0.2) 
DT=DT+(2.0^RAnER^DELTAT ) 

A3 AT=PI«OT^OT/A. 

CALL AREAS 
IFU.LE.O.O) VC = VCI 
IF( ABS(ZW).GT.O.O) GO TO 20 
IF( SUMAB.LE.O.O) GO TO 31 
X=( ABPORI+AGSLOT )/SUMAB 

90 MN0Z = AT*X/APNQZ»(2.*( l.-*-BQT/2.*PNl*PNl ) /TOP) ♦♦ZAP 
IF( AOS(MNGZ-MNU .LE. 0.002) GO TO 2 
KNl=MNOZ 
GO TO 90 

2 VNOZ=GAK^CSTAR^PNOZ^SQRT ( ( { 2 . /TCP ) ♦♦ ( TOP/BOT ) )/{ 1 .+B0T/2 . ♦MNOZ^MNO 
IZ) ) 

PRAT=(1.+0GT/2.^MNOZ*MNOZ)**(-GAM/BOT ) 

JROCK=AT/APNOZ 

8 32 SUMYA=OELY^( ACP2 + ABN2+ABS2 ) 

IF(Y.GQ.O.O) SUMYA=0.0 

VC = VC + SUi'^YA 

IF( Y.GT.0.0) GO TO ll 

P0N0Z=(0*RH0^CSTAR^SUMAB/AT)*^( l./( l.-N) )♦( l.+(CAPGAM^JR0CK)^A2/2. 
l)^^(N/(l.-N) ) 

PON = POiWOZ 

CSTAR=CSTARK^( PONOZ/IOOO. I^^CSTARP 

MDI S=AT^PCNCZ/CSTAR 

P2=P0N0Z 

P0NGZ2=P0N0Z 

PN07. = PRAr*PCNUZ 

PA = 2.^MDIS^VN0Z/( APHEAD + APNOZ )*PN0Z 
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TABLE IV-2 (CONT'D) 


IFIGRAIN.EG.:^) P^=MDIS*VNOZ/APNOZ + PNOZ 
*> PN0Z = PHAT*PC:>1CZ 

PHEAD = 2.<=NDI S*VMOZ/ { APHEAD*APNGZ )+PNOZ 

IF(GRA1,''J.£G.3) PH t A n = MD I S^VNOZ / APNOZ 4-PJOZ 

IF(PHeAD.LT.PrRAN).N = iJl 

I F ( PHEAD. L T. P i R AN ) A = A 1 

IF ( PIlEAD.Gt .PTRAN )N = N2 

IF(ini'EAD.GE.PrRAnA = A2 

HHEAD = Q<=PH£AD--S'N 

Z I T=MDI S*X/APNOZ 

RNl=RHEAO 

PHEAD2=PHEAC 

IF( punoz.lt.ptranin^.ni 
IF( PONOZ.lt. P TRAN) A = A 1 
IF(P0)N0Z.GE.PTRAN)N = N2 
IF{ PGNOZ.GE.pt RAN) A=A2 

3 RN0Z=RN1-I (RNl-0*PNOZ<'«=N-ALPHA<=ZIT<'*.8/ ( L <■«= . 2 #E XP { H ETA >PRN 1*RH0/ Z I T 
1) ) )/( l. + ALPHA«2I T*«.8«=BETA<=RH0/ZIT / I L . 2 XP I BE T A^'RN 1 *RHO/ Z I T ) ) ) ) 

IF (ABS(RNl-RNOZ).LE. 0.002) GO TO ^ 

RNl=HNOZ 
GO TO 3 

4 AVEl= (KHEAC+RNGZ ) /2. 

IF(Y.GT.O.O) GO TO 7 

RN2=RNQZ 

RH2=RHEA0 

PONJ=PGNOZ 

0PCDY=0.0 

AVE2=AVE I 

7 RNAVE = (RrJCZ + RN2) /2. 

RHAVC=(RHEAn+RH2 )/2. 

MG6N = RhfJ/2,4=( (R.NOZ + RHEAD 1 => ( A8P0RT+ABSL0T ) +2 . «0*P0N0Z*<'N<‘ABN0Z ) 
DRDY= ( AVE 1-AVE2 ) /DELY 
RBAR=(AVE1+ WC2)/2. 

GMAX = 1.0G02<=MDIS 
GMIN = 0. 99984^^015 
IF(Y.GT.O.C) GO TO 12 
GMAX=l,001*PniS 
GMIN = C.9994'RD1 S 

IF( RGEiN.GE. GRIN. AND. MGEN.LE.GKAX) GO TO 6 
MOI S=MGEM 

PGNQZ=PD IS*CSTAR/AT 
GO TO 5 
6 PONJ=PONOZ 

1 7 GAM = GAMN4'(PGN0Z/1000. )«4!GAMP 
TOP-GAM+1. 

R0T=GAM- I . 

ZAP = TOP/{ 2.<'B0T) 

CAPGAM = SORT (GAR)*(2./T0P)*<‘ZAP 
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TABLE IV-2 (CONI' D) 


ME = SQRT(2./BOT»« TOP/2. «( AE*MEl/An««( 1. /ZAP )-l. » I 
IF(ABS(ME-MEl).LE. 0.002) &0 TO 9 
ME1=ME 
GO TO 17 

9 IF(Y.Le.O.O) CALL OUTPUT 
IF(Y.LE.O.O) GO TO 10 
DELTAr=2.*DELY/(RHAVE^RNAVE) 

Z=Z+OELTAT*IRNAYE-RHAVE) 

ZO=ZQ+UELTAT*(RNAVE-RHAVE) 

T=T+DELTAT 
CALL OUTPUT 

10 IF( Y.LE..05*TAU) GO TO 16 
SINK1=VC/(CAPGAM*CSTAR)*»2*RBAR*0PCDY/12. 

MASS=.01*MDIS 

ANS^=Y+10.0*DELTAY 
IF(KOUi'JT.GT.O) GO TO 16 

IF(ABS(SINK1 ).LE. MASS. AND. ANIS4.LE.ANS-XT) GO TO 18 
GO TO 16 

18 OELY-10.*OELTAY 
GO TO 55 
16 DELY=OELTAY 
55 YL£D*Y 
Y=Y+DELY 

IFIY.GE. (TAU-XT-Z/2. l.ANO.KEHAT.EQ.O) DELY*TAU-XT-Z/2.-YLE0 
IFIY.GE. <TAU^XT-Z/2. ) .AND.KEWAT.EQ.O) V»TAU-XT-Z/2. 
IF<Y.GE.<TAU-XT-Z/2.).AN0.KEWAT.EQ.0) KEKATal 
ANS=TAU-ABS(Z/2) 

IFCY.GE.ANS.ANO.KOUNT.EQ.O) 0ELV=A.9S-YLED 
IFCY.GE. ANS.AiNID.KOUNT.EQ.OI Y=Ai^S 
DELTAT=2.*0ELY/«RHAVE+RNAVE) 

SUM2=SUMAB 
RN2=RN0Z 
RH2=RHEAD 
AVE2=AVEI 
GO TO I 

11 CSTAR=CSTARR*(PONOZ/1000.I#*CSTARP 
MOl S=AT*PONQZ/CSTAR 

GO TO 5 

12 DPCDY=(1./(1.-N> )♦( (PHEA02*PO’^OZ2)/( IABP2^AB,^2^ABS2)*2. )*DADY) 

IFC ABS(DPCi)Y) .GE.OPOUT.OK.Y.GE.XOUT) GO TO 25 

SINK1=VC/(CAPGAM*CSTAR)**2*RBAR*0PCDY/12.^IPHEAD2*P0N0Z2)/2.*CRNAV 
lE+RHAVE)/2.*( ABP2^ABN2+A8S2I/C12.*ICSTAR*CAPGAH)**2) 
STUFF=MGEN-SINK1 
MDIS=STUFF 
PONOZ=MOIS*CSTAR/AT 

IFI2.*Y+DI+0ELDI.GE.D0) PONOZ=PDNJ^DPCDY*DELY 
IF( STUFF. GE.GMIN. AND. STUFF. LE.GMAX) GO TO lA 
GO TO 5 


^ 0 - 



TABLE IV-2 (CONI' D) 


lA Pl=PONOZ 
POi'JJ=PO'J(]Z 
PC'JfJZ2= ( P 1 + P2 ) /2. 

P2 = PO;^OZ 
P3=PHEAn 

PHEAD2=(P3+P4)/2. 

PA=PHEAO 

MDI S = AT*PGNfJZ/CSTA« 

IF(KEWAr.E0.1 ) GO TO 2221 
GO TO 2222 

2221 WRl TE (6,2223) 

2223 FOR>'AT I / / , 37X, + 

237 X, • lOiriAL TAIL OFF P.FCIIS 3 7X , 

KEWAT=KEWAT+ 1 

2222 CONTINUE 
IF(Y.LT.ANS) GC TO 17 
Zw = Z 

YW=Y 

SUMBA=SUPAB 
P1=PGNGZ 
RH2=RUEAU 
RN2=RNUZ 
RAVE = AV£ 1 
ABMAIN=SURA8 
ABT0=0.0 
WRI Tt (6, 31 ) 

51 FORf^AT(//,37X, '«*«♦ FINAL 

2TAIL GFF BEGINS 3 7X ,'>«=.'*♦*««««« 4^ ) 

IF( (-1 )«=*=! .LT.C) TWl = T 

IF( (-1 I ♦<=! .LT.O) F„ i = ThKUST 

IF( I-l)*«=I.GT.O) TW2=T 

IFt (-1 M'4=I .GT.C) FW2=THRLSr 

IF ( TW2.NE.0. ) OTW^AHS ( TW2-TW1 ) 

IFI TW2.NE.O. ) CFW = ARS ( FW2-FW1 ) 

20 ANS2 = TAU + ABS( Zw/2. ) 

KOUNT = KOUNT+ 1 
IFiKGUNT.tC. 1 ) GO TO 17 
DELY'rt = DFL TAY 
0Y2=DELYW 
IF(ZK) 32,32,33 

32 IF(Y.LT.ANS2.AND.A8S(ZW).GT.DY2) GO TO 211 
SUMAB=ARPAIN 
GO TO 31 

211 SUMOY = SUyGY + nELYV, 

SUMAB = ( 1 . +SUMDY/ZW-nELYWZ (2.»ZW) ) «ABTO- ( SUMCY/ZW-DELYW/ ( 2 .«ZW ) ) ♦AB 
IMA I N 
GO TO 31 
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TABLE IV-2 (CONT'D) 


33 IF(Y.LT. ANS2.AND.2W.GT.0Y2) GQ TO 21 
SUMAB=ABTO 
GO TO 31 

21 SUMDY=SUNDY+DeLYW 

SUMAB=( l.-SUMDY/ZW+l)eLYW/(2.*ZW) )*ABMAIN+)SUWDY/ZW>DELYW/ (2.»ZW) )* 
lABTO 

31 IF( SUMA8.LE.0.0) PCN0Z=P0N0Z/2. 

IF(SUMAB.LE.O.O) GO TO 25 
CSTAR=CSTA3R>9‘(PONOZ/1000. )**CSTARP 
MDIS=AT*PCNCZ/CS1 AR 
ABAVE=( SUMAB+SUMBA) /2. 

SUMYA=DELY*ABAVE 

VC=VC+SUMYA 

DADY=( SU.VAB-SUPBA)/0£LY 
PBAR=(Pl*P0N0Z)/2. 

SUMBA=SUPAH 

22 DPCDY=PBAR/{ 1.-N)«1./A6AVE*DA0Y 
IF(P0N0Z.LE.5.0) GO TO 25 
KNOZ=Q»PCNOZ»*N 

RHEAD=RNQZ 

R0AR=(RHEAD+RAVE )/2. 

MGEN=RH04'(RN0Z + RHEA0) /2.*SUMAB 

GMA.X=l.0002*MDIS 

GMIN=0.g998#MClS 

SINKl=VC/(CAPGAM«CSTAR)<«*2*R8AR*DPCOY/l2.+PBAR#ABAVE/U2.*(CAPGAM 

♦♦CSTAR)**2)*RBAR 

STUFF=MGtN-SINKl 

IF( STUFF. GE.GMIN. AND. STUFF. LE.GPAX) GO TO 23 
MDIS=STUFF 

PONOZ=PCHJ+CPCCY*OELY 

IF( PGNOZ.LE.0.0) PaNQZ=0.0 

PBAR=(Pl+PCNOZ)/2. 

GO TO 22 

23 RHAVE=(RH2+RH6A0)/2. 

RNAVE=(RN2+KNOZ) /2. 

RH2=RHEAD 

RN2=RN0Z 

PHEA0=P0N0Z 

RAVE=RHEAD 

Pl=PONOZ 

PONJ=PQNOZ 

KDIS=AT«PCNOZ/CSTAR 

IF(ABS(DPCDY).GE.DPOUT) GO TO 25 

IF( Y.GE.XOUT) GO TO 25 

GO TO 17 
25 SUMAB=0.0 
RHEAD=0.0 
RNOZ^RHEAD 
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TABLE IV-2 (CONI' D) 


PHEAD=PONOZ 

MDIS=AT*PGNOZ/CSTAR 

WRirE(6,3l8) 

318 FORMAT{//,37X» ,/,37Xf •♦♦♦♦ BEG I 

IN HALF SECOND TRACE ♦♦♦♦•,/, 37X 
I* ) 

DELTA T=2.0*DELY/(RHAVE+RNAVE) 

T = T*^DELTAT 
CALL OUTPUT 

IF! PCNOZ.LE.0.0) GO TO 100 

TIME=T 

DELTAT=.5 

TIM=TIME+S. 

PHT=PHEAD 

PONT=PUNOZ 

SG=0.0 

29 T=T*DELTAT 

CSTAR=CSTARR*( PONOZ/ 1000. )*»CST ARP 
PHEAD=PHT/EXP(CAPGAF*<'2»AT*CSTAR/VC<‘( T-TIME)«12. ) 

PONOZ=PHEAD 
MDI S=PONOZ*AT/CSTAR 
/ Y=Y+.5*RHEAD 
CALL OUTPUT 

28 IF( T.LT.TIM.ANC.PHEAD.GE.5.0) GO TO 29 
100 WP1=G«SUNMT 

WP2=RH0*( VC-VCI)»G 
WP=(WPl+WP2)/2. 

ISP=ITOT/WP 
ISPVAC=I TVAC/rtP 
WRITE 16, 1022) 

1022 FORMAT!//, 20X, • INDIVIDUAL MOTOR DATA') 

WRI TE 16, 102) WPl ,WP2,WP,PHMAX, IXl, IX 
102 FORMAT! 13X, 'WP 1= • , IPE 1 1 .4 , / , 1 3X, • WP2= * , IPE 1 1 . A , / , 1 3X , * WP= *,0000 
llPEll.A,/, liX, *PHMAX= • , IPEll.A,/, 13X, • 1X1= • , I 1 0 , / , 1 3X , • I X= 

2110) 

NDUM=l 

IF!IPO.NE.O) CALL OUTPUT 
IFIIPO.EQ.O) GO TO 901 
IF! !-l)>»<'I.LT.O) GO TO 901 
CALL PAIR 

CALL SIGHAR! AFMAX,S1. S2, SAFMAX,BAFMAX, I , NPA I RS , SG 1 , SG2) 

CALL SIGBAR! TFMAX,SA,SR,STFMAX,BTFMAX, I , NPA I RS , SG3, SGA ) 

CALL SIGBAR! AF MAX T, S3, S4,SAFMXT, OAF MX T, I , NP A I RS , SG5 , SG6 ) 

CALL SIGBAR! TFM.AXT, SC,SO,STFMXT,BTFMXT, I , NP A I RS , SG7 , SG8 ) 

CALL SIGRARIDFTOl ,ST,SU,SDFT01 ,BDFr01, I , NPA I RS , SG19 , SG20 ) 

CALL SIGBAR!TDFT01 ,ShA,SIA,STpFTl,BTDFTl, I,NPAIRS,SG9,SGI0) 

CALL SIGBAR!DFT02,SV,SW,SDFT02,HDFT02, I ,NPAIRS,SG21 ,SG22) 

CALL SIGBAR!TDFT02,SJ,SK,STDFT2,BTDFT2, I , NP A I RS , SGI 1 , SG 12 ) 
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TABLE IV-2. (CONI' D) 


i 

I 


CALL SIGBAR(DTW,SP»SQtSDTW,OOTW,I,NPAIRS,SGl3»SGlA) 

CALL SIGBAR{FWl,SLA,SMA, SFWl ,HFWl, I , NP A I RS , SG I t SGI 6 ) 

CALL SIGBAR(FW2, SN, SO,SFW2,bFW2,I ,i\|PAIRS#SGl7,SGl8) 

CALL SIGBAR (UFw, Zl, Z2 , SDFw , 60FW, I , .NP A [ R S* SG 37 , SG 38 ) 

CALL SIGBAR(DFH0,SX,SY,SDFM0,BDFM0, I,NPAIRS,SG23,SG2A) 

CALL SIGBAR(FDIFIG, SE,SF,SFDFIG,BFDFIG, I,NPAIRStSG25,SG26) 

CALL SIGBAR(TDIFIG,SEA,SAF,STCFIG,BTDFIG, I,NPAIRS,SG27,SG28) 

CALL SIGBAR(DI r,S5,S6,Sni r.HDI r, t ,NPAIRS,SG29,SG30I 
CALL SIGBAR( ADIT,S7,S8, SADIT, 8AD1T, I , NP A I RS , SG 3 1 1 SG32 ) 

CALL SIGBAR(DFlOOK,Ol,D2,SFlOOK ,RFIOOK, I , NP A I RS t SG3 3 , SG3A ) 

CALL SIGBAR(T100K,D3,D4,ST100K,Bri00K, I , NPAI RS ♦ SG35 , SG 36 ) 

901 CONTINUE 

IF(IPO.EO.C) STOP 
WRITEI6.H87) 

887 FORMAT!//, 20X,'STANCAR0 DEVIATIONS AND MEANS FOR MOTOR PAIR DATA', 
2/, lAX, • VAR. • ,6X, * STD. DEV. *,5X,' MEAN •) 

WRI TE (6, H88 ) SAF M AX , BAFM AX , S TFMAX , BTFMA X, SAFMXT ,DAFMXT , 
2STFMxr,BTFMXT, 

2S0FTQI ,B0FT01 , STOF T 1 , BTDFT 1 , SDF T02 , BDF T02 , STDFT2 , BTDFT2 , 
2SDTW,BDTW, SFWl ,BFW1 ,SFW2,BFW2 ,SDFW.ROFvJ,SDFMG,BOFMQ, 
2SFDFIG,BFDFrG,STDFIG,BTOFIG,SDIT,BDIT,SADIT,BADIT,SF100K,BF100K, 
2Sri00K,BTlC0K 

888 FORMATI 1 JX, 'AFMAX • , iiX , IPE 1 1 . 4 , 5X , I PE 1 1 .4 , / , 

21 3X, 'TFMAX ',3X,lPE11.4,5X,lPEll.4,/, 
213X,'AFMAXT',5X,lPEll.4,5X,lP£ll.4,/, 

213X, ' TFMAXT' ,3X,1PE11.4,5X, lPEll.4,/, 

213X,'DFT01 ',5X,1PE11.4,5X,1PE11.4,/, 

213X,'TDFI01' ,SX, IPE11.4,5X, lPEll.4,/, 

21 3X, 'DF TCJ2 ',5X,1PE11.4,5X,1PE11.4,/, 

213X, • TDF102',3X, 1PE11.4,5X,1PE11.4,/, 

213X,'DTW ',5X,1PE11.4,5X, 1PE11.4,/, 

213X,'FW1 ' ,5Xy 1PE11.4,5X, lPEll.4,/, 

213X,'FW2 ' ,5X, 1PE11.4,5X, lPEll.4,/, 

213X,'DFW ',5X,1PE11.4,5X,1PE11.4,/, 

213X,'DFMU ' ,5X, 1PE11.4,5X, 1PE11.4,/, 
213X,*FDIFIG',5X,1PE11.4,5X,1PE11.4,/, 
213X,'T0IFIG',5X,1PE11.4,5X,1PE11.4,/, 

213X,'DIT ',5X,1PE11.4,5X,1PE11.4,/, 

213X, 'ADI T • , 5X, IPE 1 1.4, 3X, I PEI 1.4,/, 

213X, 'DFIOOK ',3X,1PE11.4,5X,IPEI1.4,/, 

213X,'T10UK ' ,5X, IPE11.4,3X, lPEll.4) 

WRI TE (6, f389) SG 1 , SG2 , SG5 , SG6 

889 FORMAT!//, 20X, 'ALTERNATE DISPERSION VALUES FOR THRUST IMBALANCE DA 

2TA' ,/,14X, 'VAR. • ,6X, ' SIGMA I ',t>X,' SIGMA 2 ',/, 

313X,'AFMAX ’ , 3X , 1 PE 1 1 . 4 , 5X , 1 P E 1 1 .4 , / , 1 3X , ' AF MA XT ' , 5X , 1 PE 1 1 . 4, 

45X, 1PE11.4) 

► IFIIOP.NE.O) CALL PLOT !0. 0,0. 0,999) 

STOP 

END 
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TABLE IV-2 (CONT’D) 


SUBROUTINE AREAS 

C ♦ SUBROUTINE AREAS CALCULATES BURNING AREAS AND PORT AREAS FOR * 

C * CIRCULAR PERFORATED IC.P.) GRAINS AND STAR GRAINS OR FOR A ♦ 

C ♦ COMBINATION OF C.P. AND STAR GRAINS ♦ 

INTEGER STAR, GRAIN, ORDER, COP 

REAL MGE.N,MDIS,MNOZ ,MN1 , JROCK , N , L , ME I , ME , I S P t I TOT , MU , MASS , I SPVAC 

REAL LGCI ,LGNI , NS , NN , NP , LGS I , NT , LTP , LGC , LS , LF 

REAL M2,MDBAR, ISP2, I T VAC , L 1 , L2 , LF W , LF WSOD 

COMMON/CONSTl/ZW, AE,AT, T HE T A, AL FAN 

C 0 M M 0 N / C 0 N S T 3 / S , N S , G R A I N 

COMMON/CONS TA/UELD I , 00 , D I , ZC , X T , ZO 

COMMON/VARI A1/T,OELY,OELTAT,PONOZ,PHFAD,RNUZ,RHEAD, SUMAB,PHMAX 
COMMON/VARI A2/AHP0RT , ABSLOT, ABNOZ,APHEAD, APNOZ, DADY,ABP2, A6N2,ABS2 
C0MM0N/VARIA3/ITUT, I T VAC , JROC K , I SP , I SP V AC , MD 1 S , MNOZ , SG, SUMMT 
COMMON/VARI AA/RNT,RHT,SUM2,R1 ,R2,R3,RHAVE,RNAVE,RBAR,YB,K0UNT 
COMMON/VARI A5/ ABMA IN, A8 TO, SUMDY, VC 1 , VC 
COMMON/VARI A6/YDI , TE 
COMMON/VARI A 7/ Y, THRUST 

C0MM0N/0VALA/CHIH,CHIN, SEN,SEH,CHINH, AZ,BZ 
COMMON/DATAZ/ IDATA 
DATA PI/3.14159/ 

21111 F0RMATIE16.9) 

A8PC=0.0 

ABNC=0.0 

ABSC=0.0 

ABPS=0.0 

ABNS=0.0 

ABSS=0.0 

DABT=0.0 

SG=0.0 

VCIT=0.0 

ANUM=PI/4. 

PID2=PI/2. 

RNT=RNT+RNOZ*OELTAT 
RHT=RHT+RHEAD«DELTAT 
IF(Y.LE.O.O) AGS=0.0 
K=0 

IF(ABSIZW).GT.O.O) K=1 
YB=Y 

IFIK.EQ.l) Y=YB-SUMDY/2. 

2 IF(K.£0.2) Y=YB+ABS ( ZW) /2.-SUMDY/2. 

IFIIOATA-l) 5000,5000,5001 

5000 IFIY.LE.0.0) READ(5,500) I NPUT , GRA I N , S T AR , NT , ORDER , COP 
IF(Y.LE.O.O) WRITE(2,500) I NPUT , GRA I N, ST AR , NT , ORDER , COP 
GO TO 5002 

5001 IF(Y.LE.O.O) READ(2,500) INPUT , GRA IN, STAR , NT , ORDER, COP 
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TABLE IV-2 (CONI' D) 


C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


5002 CONTINUE 


* 

* 

* 

♦ 

« 

« 

« 

« 

♦ 

♦ 

♦ 

♦ 

♦ 

* 

« 

4c 

« 

« 

« 

« 

« 

* 


♦ 

♦ 

♦ 

♦ 

♦ 

4e 

* 


READ THE TYPE OF INPUT 
CONFIGURATION AND 
VALUES FOR INPUT ARE 

1 FOR ONLY TABULAR INPUT 

2 FOR ONLY EQUATION INPUTS 
INTO THE SUBROUTINE) 

3 FOR A COMBINATION OF I AND 
FOR GRAIN ARE 

1 FOR STRAIGHT C.P. GRAIN 
STRAIGHT STAR GRAIN 
COMBINATION OF C.P. 

ARE (WAGON WHEEL IS 
THIS PROGRAM) 

STRAIGHT C.P. GRAIN 
STANDARD STAR 
TRUNCATED STAR 
WAGON WHEEL 


FOR THE PROGRAM AND THE BASIC GRAIN 
ARRANGEMENT 


(EQUATIONS ARE BUILT 


VALUES 


VALUES 

STAR 


VALUES 


2 FOR 

3 FOR 
FOR STAR 
GRAIN IN 

0 FOR 

1 FOR 

2 FOR 

3 FOR 
FOR NT ARE 

0 IF THERE 
X WHERE X 


AND STAR GRAINS 
CONSIDERED A TYPE 


OF 


ARE NO TERMINATION PORTS 
IS THE NUMBER OF TERMINATION 


VALUES OF ORDER ESTABLISH HOW 
GRAIN IS ARRANGED 

1 IF DESIGN 

2 IF DESIGN 

3 IF DESIGN 
A IF DESIGN 

♦ ♦♦NOTE^v^' IF 
♦♦♦NOTE*** IF 


A COMBINATION C.P, 


PORTS 
AND STAR 


IS 

IS 

IS 

IS 


star 
C.P. 
C.P. 

STAR 
GRAIn = 1 , 
GRAI N=2, 


AT 

AT 

AT 

AT 


HEAD 
HEAD 
HEAD 
HEAD 
VALUE OF 
VALUE OF 


END 

END 

END 

END 


AND 
AND 
AND 
AND 
ORDER 
ORDER 


CONTINUE 

VALUES 


C.P. 

C.P. 

STAR 

STAR 

MUST 

MUST 

ONLY) 


AT 

AT 

AT 

AT 

BE 

BE 


NOZZLE 

NOZZLE 

NOZZLE- 

NOZZLE 

2 


AFT END IS 


END IS 


♦ 

♦ 

♦ 

* 

* 

♦ 

* 

♦ 

♦ 

♦ 

« 

* 

* 

* 

♦ 

♦ 

« 

♦ 

♦ 

♦ 

♦ 

* 

* 


FOR CO? ARE (APPLICABLE TO C.P. GRAINS 

0 IF BOTH ENDS ARE CONICAL OR FLAT 

1 IF HEAD END IS CONICAL OR FLAT AND 
HEMISPHERICAL 

2 IF BOTH ENDS ARE HEMISPHERICAL 

3 IF HEAD END IS HEMISPHERICAL AND AFT 
CONICAL OR FLAT 

500 FORMAT (9X, I2,9X,I2,3X,I2,6X,FA.0t9X,I2,7X,I2) 
IF(Y.LE.n.O) WRITE(6,607) 

607 FORMAT!//, 20X, 'GRAIN CONFIGURATION') 

IF(Y.LE.O.O) WRITE(6,600) I NPUT, GRA IN, STAR ,NT ,ORDER ,COP 


600 F0RMAT(13X, ' INPUT= ' , I 2 , / , 1 3 X , • GRA I N= ',12,/ 
l,'NT= • ,F^.O,/ ,13X, 'ORDER= ' , I 2 , / , I 3X , ' COP = 
IF! INPUT. EQ. 2) GO TO 12 
IF(Y.LE.O.O) GO TO 6 
IF ( YT.LE.Y. AND.K.LT .2 ) GO TO 8 


13X, ' STAR= 

f 12,// ) 


*,I2,/,13X 
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TABLE IV-2 (CONI' D) 


9 DENOM=YT-YT2 

SLOPEl=f ABPK-ABPK2)/DEN0M 
SLOPE2={ ABSK-ABSK2) /DENOM 
SLOP£3= ( A8rJK-ABNK2 ) /DENOM 
SL0PE4=( APHK-APHK2 )/DENOM 
SLOPE5= ( APNK-APNK2 ) /UENQM 
Bl=ABPK-SLOPEl*YT 
B2 = ABSK-SLOPE2<=YT 
B3=A6NK-SLOPE3*YT 
BA=APHK-SLOPEA*YT 
B5=AP:MK-SL0PE5«YT 
ABPT=SLOPEl*Y+Hl 
ABST=SLQPE2«Y+B2 
ABNT=SLOPE3*Y+83 
APHT=SL0PEA*Y+B4 
APNT=SLOPE54«Y + R5 
IFIINPUT.E0.3) GO TO 3 
GO TO 52 

6 IF(IDATA-l) 5003,5003,5004 

5003 READ(5,5071 YT , APPK , ABSK, ABNK, APHK, APNK, VCl T 
WRITE (2, 507) YT , ABPK , ABSK, ABMK , APHK , APNK , VC I T 
GO TO 5005 

5004 REA0(2,507) YT , ABPK , ABSK, A8NK, APHK, APNK, VC I T 

5005 CONTINUE 


C * READ IN TABULAR VALUES FOR Y=0.0 (NOT REQUIRED IF INPUT=2) ♦ 
C ♦ ♦ 
C ♦ ABPK IS THE BURNING AREA IN THE PORT IN IN**2 ♦ 
C ♦ ABSK IS THE BURNING AREA IN THE SLOTS IN 1N<«42 ♦ 
C * ABNK IS THE BURNING AREA IN THE NOZZLE END IN IN4*2 ♦ 
C ♦ APHK IS THE PORT AREA AT THE HEAD END IN IN^^Z ♦ 
C ♦ APNK IS THE PORT AREA AT THE NOZZLE END IN IN**2 ♦ 
C ♦ VCIT IS THE INITIAL VOLUME OF CHAMBER GASES ASSOCIATED WITH ♦ 
C ♦ TABULAR INPUT IN IN*«3 ♦ 


507 F0RMAT(6X,F6.2,l0X,Fll.2,l0X,Fll,2,8X,Fll.2,/,22X,Fll.2,9X,FU.2, 
1 6X,F11.2) 

WRITE(6,610) 

610 FORMAT! 13X, • TABULAR VALUES FOR YT EQUAL ZERO READ IN') 

KRI TE(6,583) ABPK , ABSK, ABNK, APHK, APNK 

583 FORMAT! 1 3X, *ABPK=* , IPEl 1 .4 , 5X , * ABSK= • , I PE 1 1 .4, 5X , • ABNK=* , IPE 1 1.4, 

1 5X, 'APHK=' , lP£ll.4,5X, •APNK=»,lPEll.4) 

WRITE(6,584) VCIT 

584 FORMAT! 13X, •VCIT=« , IPEll. 4,//) 

ABPT=ABPK 

A8ST=ABSK 

ABNT=ABNK 

APHT=APHK 
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TABLE IV-2 (CONT'D) 


APfU = APNK 
YT2=YT 

IF( IiMPUT.EQ.3) GO TU 3 
vci=vcir 
GO TO 52 
8 YT2=Yr 
A8PK2=AQPK 
AB,NK2=ABMK 
ADSK2=A8SK 
APHK2=APHK 
APNK2=APNK 

TF(lDATA-l) 5006,5006,5007 

5006 READ(5,505) Y f , APPK , ABSK , ADrjK , APHK, APNK 
WRI TE (2 , j 05 I YT, AOPK, AbSK,AQ:'^K,APHK, APNK 
GO TO 5008 

5007 READ(2,5U5) YT , A HPK , ABSK , ABNK , APHK , APNK 

5008 CONTINUE 

C ♦ REAU IN TABULAR VALUES FOR Y=Y (NOT REQUIRED FOR INPUT=2) ♦ 

505 F0RMAT(6X,F6.2, 10X,F11.2,10X,Fll.2,8X,Fll.2,/,22X,Fll.2,9X,Fll.2) 
WRITE(6,6ll) YT 

611 FORMAT! ///, 13X,» TABULAR VALUES FOR YT= •,F7.3,' READ IN* I 
WRIT£(6,583) ABPK, ABSK, ABNK, APHK,APNK 
GO TO 9 
12 ABPT=0.0 
ABNT=0.0 
ABST=0.0 

3 IF( GRAIN.NE.2) go TO 4 
ABPC=0.0 

ABNC=0.0 
ABSC=0.0 
GO TO 7 

4 IFIIDATA-l) 5009, 5009,,5010 

5009 IF(Y.LE.O.O) READ(5,50l) XTZO,S 
IF(Y.LE,0.0) WRITE(2,501) XTZO.S 
GO TO 5011 

5010 IF(Y.LE.0,0) READ(2,50l) XTZO,S 

5011 CONTINUE 

IF( Y.LE.O.O ) RE AO (4,2 1 1 1 1 ) DO , D I , THE TAG, LGC I , L GN I , THE TCN, THETCH 

Q *>*****#**** 


c * READ IN BASIC GEOMETRY FOR C.P. GRAIN (NOT REQUIRED FOR ♦ 
C ♦ STRAIGHT STAR GRAIN) * 
C ♦ XTZO IS THE DIFFERENCE BETWEEN THE I N 1 1 1 AL - I NT ERN AL GRAIN ♦ 
C » DIAMETER AT THE NOZZLE END OF LGCI AND DI IN INCHES ♦ 
C ♦ LESS TWICE XT AND LESS ZO ♦ 
C * S IS THE NUMBER OF FLAT BURNING SLOT SIDES (NOT INCLUDING ♦ 
C « THE NOZZLE END) ♦ 
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TABLE IV-2 (CONT’D) 


C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


♦ ♦ 

♦ THE FOLLOWING VARIAOLES ARE OBTAINED FROM THE STATISTICAL ♦ 

♦ analysis program ♦ 

« « 

♦ DO IS THE AVERAGE OUTSIDE INITIAL GRAIN DIAMETER IN INCHES ♦ 

♦ DI IS THE AVERAGE INITIAL INTERNAL GRAIN DIAMETER IN INCHES ♦ 

♦ THETAG IS THE ANGLE THE NOZZLE END OF THE GRAIN MAKES WITH ♦ 

♦ THE MOTOR AXIS IN DEGREES ♦ 

♦ LGCI IS THE INITIAL TOTAL LENGTH OF THE CIRCULAR PERFORATION ♦ 

♦ IN INCHES * 

♦ LGNI IS THE INITIAL SLANT LENGTH OF THE BURNING CONICAL ♦ 

♦ GRAIN AT THE NOZZLE END IN INCHES * 

♦ THETCN IS THE CONTRACTION ANGLE OF THE BONDED GRAIN IN DEGREES* 

♦ THETCH IS THE CONTRACTION ANGLE AT THE HEAD END IN DEGREES ♦ 


501 FORMATI6X.FIO. i,3X,F10.0) 

IF(Y.LE.O.O) WRITE(6,60l ) DO, DI , XTZO. S , THETAG, LGCI , LGN I * THE TCN , TH 
lETCH 

601 F0RMATI20X, 'C.P. GRAIN GEOMETRY* ,/, 13X, • D0= • , F7. 3 , / , 1 3X , *DI= *,F7 
1.3,/, 13X, 'XTZO= • ,F7.3,/, 13X, *S= • , F4.0 , / , I 3X , • THETAG= *,F8.5,/,13 
2X,*LGCI= * ,F7.2,/, 1 3X, *LGNI= • , F6. 2 , / , 1 3X , • THETCN= • , F 8. 5 , / , 1 3X , 
3*THETCH= ',FR.5,//) 

IF(Y.LE.O.O) TAU=(D0-0n/2.0 

IF( Y.LE.0.0) DELOI=2.0*XT+ZO+XTZO 

IF (Y.LE. 0.0) THETAGsTHET AG/57. 29578 

IFI Y.LE.0.0) rH£TCN=TMETCN/57.29578 

IF (Y.L£.0.0)TH£TCH=THFTCH/57. 295/8 

00SQ0=0()*0Q 

DISQD=DI*DI 

BNUM=ANUM*DOSQD 

TLL=TE 

IFI0RDER.GE.3) TLL=0.0 
YDI=2.*Y*0I 
YDISQD=YD1*YDI 
ABSC=S*ANUM*(DOSQO-YDISQO) 

IF! ABSC.LE.0.0 ) ABSC = 0.0 
IF( YDI.GE.OOIGO TO 100 
IF( THETAG. GT. 0.08727) GO TO lOl 
IF(COP.EQ.O) GO TO 700 
IF(COP.EO.l) GO TO 701 
IFIC0P.EQ.2) GO TO 702 
CHCK1=DOSOD-YDISOD 
IFICHCKI.lt. 0.0) CHCK1=0.0 

LGC=LGC I - ( SORT I DOSQO-D I SQD ) -SORT I.CHCKl) ) /2.-Y*C0TAN I THETCN ) 

GO TO 710 

702 CHCKI=DOSQD-YD1SOD 
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TABLE IV-2 (CONT'D) 


IF(CHCKl.Lr.O.O) CHCKI=0.0 
LGC=LGCI-(SQRT(DOSOD-DISGn)-SORT(CHCKl ) ) 

GO TO 710 

701 CHCK2=00SUD-( YDI+DELDI )**2 
IF(CHCK2.LT.0.0) CHCK2=0.0 

LGC = LGCI-(SORT(DOSOD-(DI+nELDI ) ) -SORT ( CHCK2 ) )/2. 

2-Y«'C0TA.J( THETCH) 

GO TO 710 

700 LGC = LGCI-Y*(COTAN( THE TCN) ♦COTAN (THE TCH) ) 

710 ABPC = PI*YDH'(LGC-TLL-S»Y) 

ABNC=0.0 
GO TO 732 
lOl CONTINUE 

IF(COP.EO.O.OR.COP.EO.l) GO TO 720 
CHCKl=DOSQD-YniSQO 
IF(CHCKl.LT.O.O) CHCK1=0.0 

LGC=LGCI-( SORT(DOSOn-OISOU)-SORT(CHCKl ) )/2.-TLL 
2-( S+TANI ThETAG/2. ) ) «Y 
ABPC=PI«YDI*LGC 
GO TO 730 

720 LGC=LGCI-Y«COrAN( THETCH)-TLL-(S+TAN(THETAG/2. ) l«Y 
ABPC=PI*YDI*LGC 

730 IFICOP.EO. 1.0R.C0P.EQ.2) GO TO 731 

ABNC = PI *(LG JI-Y^C0TAN( THETAG^THETCN)-Y*TANI THETAG/2. ) V*(DI + 

1 06Lni4Y+LGNI*SlNl THETAG)+Y*SIN(THETCN)/SIN( THETAG+THETCN) ) 

GO TO 732 

731 IF(Y.LE.O.O) GO TO 7311 
GO TO 7312 

73U R7=( (DI+OELUI ) /2 . + L GN I ♦ S I N ( T H E T AG) ) *COS ( T HE T AG ) -S I N ( THE TAG ) ♦ 

1 SORT( (00/2. )**2-( (OI^OELDI ) /2.+LGN I «S I N ( THET AG ) ) ♦♦2 ) 

7312 IF(R7+Y.LT. (DQ/2. )«C(JS(THETAG) ) GO TO lllll 

ABNC=PI» ( LGNl+( I ./SIN( THEFAG) ) ♦ ( (DO/2. ) -L GN I I N ( THET AG )- ( D I 
2-DELDI ) / 2. )- Y>»=CO T AN ( THE FAG)- Y*COTAN( THE TAG/2. ) ) » ( ( D I -DELD I ) /2 . 
3+Y+00/2.) 

GO TO 22222 

lllll RPR = SORT( ( (DO/2. ) *<=2 ) -R7«*2 ) -SQRT( ( { DO/2 . ) ♦♦2 ) - ( R7 + Y ) *#2 ) 

ABNC=PI*( LGNI-RPR-Y^TANI THETAG/2. ) )♦( ( DI+DELDI ) /2. ♦ SQR T ( ( DO/ 

1 2. )<'*2-(R7 + Y)«^2)*SIN(THEFAG)+Y + (R7 + Y)«CnS(THETAG) ) 

22222 CONTINUE 

732 IF( ABPC.LE.0.0) ABPC=0.0 
IF( A3NC.LE.0.0) A8nC=0.0 
GO TO G 

100 ABNC=0.0 
A8PC=0.0 

p APHT = ANUM*(0I+2. <'RHF )**2 
IF( APHT.GE.BNUM) APHT=BNUM 
IFIK.LT.2) APHTl=APHF 
APNF = ANUM<'(DI>DELDI ♦2 .*RnT )«*2 
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TABLE IV-2 (CONI' D) 


IF( APNT.GE.6NUM) APNT=BNUM 

IFlGRAIN.NE.l) GO TO 7 

ABPS=0.0 

ABSS=0.0 

ABNS=0.0 

GO TO 50 

7 IF(IDATA-l) 5012,5012,5013 

5012 IF(Y.LE.O.O) REA0(5,502) NS,NP,NN 
IFlYiLE.0.0) WRITE(2,502 )' NS,NP,NN 
GO TO 501A 

5013 IF(Y.LE.O.O) REA0(2,502) NS,NP,NN 

5014 CONTINUE 

if(y.le.o.o) read(4, 21111) lgsi,rc, fill 


C ♦ READ IN BASIC GEOMETRY FOR STAR GRAIN (NOT REQUIRED FOR ♦ 

C ♦ STRAIGHT C.P. GRAIN) ♦ 

C ♦ NS IS THE NUMBER OF FLAT BURNING SLOT SIDES (NOT INCLUDING ♦ 

C ♦ THE NOZZLE END) ♦ 

C ♦ NP IS THE NUMBER OF STAR POINTS ♦ 

C ♦ NN IS THE NUMBER OF STAR NOZZLE END BURNING SURFACES ♦ 

C ♦ ♦ 

C * the following variables are OBTAINED FROM THE STATISTICAL » 

C ♦ ANALYSIS PROGRAM ♦ 

£ 4i4E4> 4= 4c 4c 4c 4^41^4(4(4: 41* 4c 4 4 4>4t 4: 44>4c4c 4 : 

C ♦ 4c 

C * LGSI IS THE INITIAL TOTAL LENGTH OF THE STAR SHAPED 4 « 

C ♦ PERFORATED GRAIN IN INCHES * 

C * RC IS THE AVERAGE STAR GRAIN OUTSIDE RADIUS IN INCHES ♦ 

C 4C FILL IS THE FILLET RADIUS IN INCHES * 


C 4c44i4[4(4e4(4c«4c4(4:4(4(4(4(4c4(*4(4:4(4(*4(4(4i«*4:4(4(4(4(4(4(4c4(4(4c4[4c4(4(4(4(4c4(4(44(*4i*4>44(*«4>«*44(4(4>*4(4( 

502 FORMAT (4X, FI 0.0, 4X, FI 0.0, 4X,F 10.0) 

IF(Y.LE.O.O) WRITE(6,602) NS , LGS I ,NP , RC , F I L L ,NN 
602 F0RMAT(20X, 'BASIC STAR GEOME TRY* , / , 1 3X , • NS= • , F4 . 0 , / , I 3X, • LGS 1= », 
IF7.2,/,1:3X, *NP= ',F4.0,/,13X, *RC= • , F7. 3 , / , 1 3X, • F I L L= »,F7.3,/,13X 
2,»NN= *,’F4i0,//) 

IF( y.le.;o.6.and.grain.eo.2) DO=2.0*RC 

PIDNP = Pl’/NP 

RCS0D=RC4'RC 

FY=FILL+Y 

fysqd=fy*fy 
IF(STAR.EQ. 1 ) GO TO 20 
IF(STAR.E0.2) GO TO 201 
IF(Y.GT.O.O) GO TO 179 

READ(4,21111 ) RI WW,Ll ,L2,ALPHA1, ALPHA2,HW 

C :^44^4c4c4i*4c4[4[:4*4c4c4(*44c4c4c4*4i44(4(*44:44(4:4 4 44(* 4 4( 4c * 4 4(4: * *4(4(4i4(4(4(4( **4( 4c 4* 444444c 4(4: 4( 

C 4 READ IN GEOMETRY FOR WAGON WHEEL (NOT REQUIRED FOR STANDARD 4 

C * OR TRUNCATED STAR GRAINS) ♦ 
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TABLE IV-2 (CONT'D) 


C ♦ ♦ 

Q «««««« «««««« ««««««« ««««««« ««««««*« 

C * THE FfJLLCWiNG VARIAflLES ARE OBTAINED FROM THE STATISTICAL ♦ 

C ♦ ANALYSIS PROGRAM ♦ 

c ♦ ♦ 

C ♦ RIWW IS THE AVERAGE RADIUS OF THE INSIDE OF THE PROPELLANT ♦ 

C ♦ WEB IN INCHES ♦ 

C ♦ LI AND L2 ARE THE LENGTHS OF THE TWO PARALLEL SIDES OF THE ♦ 

C ♦ TWO SETS OF STAR POINTS IN INCHES ♦ 

C ♦ ALPHAl AND ALPHA2 ARE THE ANGLES BETWEEN THE SLANT SIDES OF * 

C ♦ THE STAR POINTS CORRESPOND! NG TO Ll AND L2» RESP6CTI VELY, « 

C ♦ AND THE CENTER LINES OF THE POINTS IN DEGREES ♦ 

C ♦ HW IS HALF THE WIDTH OF THE STAR POINTS IN INCHES * 


WRI TE (6,^22 ) RIWVULl ,L2t ALPHAl,ALPHA2,HW 
A22 F0RMATI20X, 'WAGON WHEEL GEOMETRY* ,/. I 3X R I WW= • ,F5 .2 , / , 1 3X , 

1 *L1= • ,F5.2»/, 13X, *L2= • , F 5. 2 , / , I 3X, • ALPHA 1= • ,F7.5, /, 13X, 

2 *ALPHA2= ',F7.5,/,l3Xt *HW= ',F5,2.//I 
IFIY.LE.0.0) TAUWW=RC-RI WW 

IF( Y.LC.C.O.AND.GRA IN.EQ.2) DI=D0-2.0*TA0WW 

ALPHA1 = ALPHA 1/57, 20 578 

ALPHA2= AL PH A2/57. 29578 

ALP2=ALPHA2 

XL2=L2 

LFW=RC-TAUWW-FILL 
LFWSOD=LFw<'LFW 
THETFW=ARSINI I HW+F I LL ) /LFW ) 

SLFW=LFW«SIN(THETFw) 

179 KKK=0 
SG=0.0 

ENUM=(KCSQD-LFWSQD-FYSQ0J/(2.*LFW*FY) 

ALPHA2=ALP2 

L2=XL2 

190 YTAN=YOTAN( ALPHA2/2. ) 

COSALP=COS { ALPHA2) 

SINALP=SINIALPHA2 ) 

IFI YTAN.GT.L2) GO TO 182 
IF( FY.GT.SLFW) GO TO 181 

SGW=NP*(L2-2.«YTAN+ ( SLFW-FILL ) / S I NALP-Y*COT AN ( ALPHA2 ) +F Y* 

1 IPID2+THETFW)+(LFW+FY)*|PI0NP-THETFW) ) 

GO TO 183 

181 IFI Y.GT. TAUWW) GO TO 18^r 

SGW=NP«( FY*(PI0NP4^ARS INISLFW/FYI ) ♦ I P I DNP-THETFw ) *LF W ) 

GO TO 183 

18A SGW^NPOFY^ITHETFW+ARSINI SLFW/FY)-ARC0S<ENUM)) 

GO TO 183 
102 YPQ=-SLFW 
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TABLE IV-2 (CONI' D) 


IF(ALPHA2.GE.PID2) GO TO 222 
Q=-F1LL+L2*TAN(ALPHA2 )-Y/COSALP 

XPI =(-Q*TAN( ALPHA2 )-SQRT (-0*Q+FYSQD/COSALP*COSALP ) ) *COSAl P^COSALP 
YPI=XPI*TAN( ALPHA2)+Q 
XPO=( YPn-Q)«CQTAN( ALPHA2) 

GO TO 223 

222 XPI=Y-L2 

YPI=-SQRnFYSOD-XPI«'XPI ) 

XPO=XPI 

223 FYLS = SQRT(SLFW<'SLFW+XPI#XPn 
XPI02 = (XPI-XP0)<=(XPI-XP0) 

YP102=( YPI-YPO)«lYPI-YPO) 

IF(FY.GT.FYLS) GO TO 186 

IF( Y.GE.TAUWW) GO TO 185 

SGW=NP#( SORT( XPI02+YPI02 )+FY*(PID2+THETFW-ARSIN{XPI/FYn+CLFW+FY)* 
I (PIDMP-THETFWn 
GO TO 183 

185 SGW=NP<'( SORT (XPI02 + YPr02}+FY*(PID2-AHSIN(XPI/FY)-ARC0S( ENUM) ) ) 

GO TO 183 

186 IF( Y.GT.TAUWW) GO TO 187 

SGW=NP#( FY*(PIDNP+ARSIN( SLFW/FY) )+ ( P I DNP-THET FW ) ♦LF W ) 

GO TO 183 

187 SGW=NP«FY*( THETFWf ARSIN( SLFW/FY)-ARCOS(ENUM) ) 

183 IF(SGW.LE.O.O) SGW=0.0 

IF(Y.GT.O.O) GO TO 188 

AGS2=.5«( PI*RCSQD-NP*LFW#SLFW*(COS(THETFW)-Sl N( THETF W ) *COTAN I ALPHA 

1 2)-2.*( L2+FILL*rAN( ALPHA2/2. ) I /LFW ) - ( P I -THETF W^NP ) *LF WS0D-2.*NP*F 

2 ILL»(L2 + SLFW/SINALP+LFW#(PI0NP-THETFW) + (PIDNP + PrD2-l./Sr.^<ALP)« 

I FILL/2. ) ) 

AGS=AGS*AGS2 

188 continue 

SG=SG+SGW 

IF(KKK.EQ.l) GO TO 2A 
L2=L1 

ALPHA2=ALPHAI 
KKK = 1 

. GO TO 190 

201 IF(Y.LE.O.O) REA0(4, 211111 RP,RIS 

C ♦♦♦♦♦♦♦♦♦♦■♦♦♦♦♦♦♦♦♦♦♦A ♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦ ♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦ 


C * READ IN GEOMETRY FOR TRUNCATED STAR I NOT REQUIRED FOR ♦ 

C ♦ STANDARD STAR OR WAGON WHEEL) ♦ 

C ♦ ♦ 

C ♦ THE FOLLOWING VARIABLES ARE OBTAINED FROM THE STATISTICAL ♦ 

C ♦ ANALYSIS PROGRAM ♦ 

C ♦ ♦ 

C ♦ RP IS THE INITIAL RADIUS OF THE TRUNCATION IN INCHES ♦ 
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TABLE IV-2 (CONI' D) 


C ♦ RIS IS THE AVERAGE RADIUS OF THE INSIDE OF THE PROPELLANT ♦ 

C ♦ WEB IN INCHES * 

C ♦ OF THE SLOTS IN INCHES ♦ 

IF(Y.LE.O.O) WRITE(6,603I RP.RIS 

603 F0RMAT(20X, 'TRUNCATED STAR GEOMETRY *,/, 13X, »RP= » ♦ F7. 3, /# 13X. *R I S= 
I •,F7.3,//) 

IFIY.LE.0.0) TAUS=RC-RIS 

IFI Y.LE.O.O.AND.GRA IN.EQ.2) OI=DO-2.0*TAUS 

THETAS=PIDNP 

RPY=RP+Y 

LS=RC-TAUS-FILL-RP 

RPL=RP+LS 

THETSl=THETAS-ARSIN(FY/RPY) 

IF( THETSI.LE.0.0) GO TO llO 
IFI Y.LE.TAUS) GO TO 103 

THETAC=ARSIN( (RCS0D-RPL*RPL-FYSQDI/(2.*FY*RPL) ) 

IF(THETAC.GE.O.O) GO TO lOA 
IFCY.LT.RC-RP) GO TO 105 
SG=0.0 
GO TO 14 

103 SGs2.*NP*(RPY*THETSl+LS-(RPY*C0SCTHETAS-THETSl)-RP)^PID2*FYJ 
GO TO 14 

104 SG=2.*NP=9'(RPY«'fHETSULS-(RPY*COS(THETAS-THETSl)-RP) + FY*THETAC) 

GO TO 14 

105 SG=2.<‘NP*(RPY*THETSUSQRT(RCSQD-FYS0D)-S0RT(RPY*RPY-FYS0D) ) 

14 IFI Y.LE.0.0) AGS = PI*IRCS0D-RP*RP)-NP*IPI*FILL<‘FILL/2.<-2.<'LS«FILL) 
GO TO 31 

110 THETAF=THETAS 

THETAP=2. ★THETAS 
TAUWS=TAUS 
GO TO 111 

20 IFIY. LE.O. 01 REA0{4, 21111) THETAF, THETAP» Rl WS 


C ♦ READ IN GEOMETRY FOR STANDARD STAR INOT REQUIRED FOR ♦ 

C ♦ TRUNCATED STAR OR WAGON WHEEL) ♦ 

C ♦ ♦ 

C ★ THE FOLLOWING VARIABLES ARE OBTAINED FROM THE STATISTICAL ★ 

C ★ ANALYSIS PROGRAM ♦ 

C ♦ ♦ 

C ♦ THETAF IS THE ANGLE LOCATION OF THE FILLET CENTER IN DEGREES ★ 

C ♦ THETAP IS the ANGLE OF THE STAR POINT IN DEGREES ★ 

C ★ RIWS IS THE AVERAGE RADIUS OF THE INSIDE OF THE PROPELLANT ♦ 

C ★ WEB IN INCHES ♦ 

Q ★«★★★«★★★★★★★★★★ ★★«★★★★★«★★★★★★«★★★★*»★*★ ★★★★★★★★★ 

IFIY.LE.0.0) WRITE(6,604) THETAF , THE TAP , R I WS 
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TABLE IV-2 (CONT’D) 


60^ F0RMAT(20X, 'STANDARD STAR GEOMETRY* ,/, 13X, • THETAF= • , F7.5,/ , I 3X, • T 
XHETAP= • ,F7.5»/, 13X, 'RIWS= •♦F6.3,//) 

IF(V.LE.O.O) TAUWS=RC-RIWS 

IF( Y.LE.O.O.A 'JD.GRAIN.EQ.2) DI = D0-2.0^TAUWS 
THETAF=THETAF/57.29578 
THETAP=THETAP/ 57. 29578 
THETAS=PI/NP 
THETSl=l.OO 
111 LF=RC-TAUWS-FILL 
CNUM=(Y+FILL)/LF 
DNUM=SIN( THETAF) /SI N( THE TAP/2. ) 

ENUM=(RCS0D-LF*LF-FYS0D)/(2.<'LF«FY) 

FNUM=SIN( THETAF )/COS( THE TAP/ 2. ) 

IF(CNUM.LE.FNUM) GO TO 106 
IF( Y.LE. TAUWSIGO TO 107 

SG=2.*NP*FY*(THETAF+ARSIN(SIN(THETAF)/CNUMI-ARC0S(ENUMI ) 

GO TO 23 

106 IF( Y.LE.TAUWS) SG=2 .♦NP^LF* ( DNUM+CNUM* I PI D2 + THETAS-THETAP/2 . 
l-COTANlTHETAP/2. ) ) + THETAS-THETAF ) 

IF(Y.LE.TAUWS) GO TO 23 

SG=2.*NP#(FY«ARSIN( ENUM-ITHETAS-THETAP/2. J » +LF*ONUM-F Y^COTAN ( THETA 
lP/2. ) ) 

GO TO 23 

107 SG=2.<'NP«LF«(CNUM*( THETAS + ARSIN(SIN(THETAF)/CNUM) I + THETAS-THETAF ) 

23 IF(THETSl.LE.O.O) GO TO 1^ 

IF( Y.LE.0.0) AGS=PI#RC*RC-NP*LF*LF*( SIN (THETAF) ♦(CDS! THETAF )- 
ISINl THETAF )*COTAN( THE TAP/2. ) THETAS-THETAF+2 . ♦F ILL /LF* ( S I N ( THETAF 
2)/SIN(THETAP/2.)+THETAS-THETAF+FILL/(2.*LF)*(PID2+THETAS“THE 
• 3TAP/2.-C0TAN(THETAP/2. n ) ) 

24 CONTINUE 

31 IF( SG.LE.0.0) SG=0.0 

IF(K.E0.0.0R.K.E0.2) SGN=SG 
IF(K.LE.l) SGH=SG 
IF(Y.LE.O.O) SG2=SG 
IFIK.E0.2) GO TO 37 
RAVEDT=Rl+(SG+SG2)/2.*RBAR*DELTAT 
. RNDT=R2+(SG+SG2)/2.#RNAVE»0ELTAT 
RHDT=R3+( SG+SG2)/2.«RHAVE*DELTAT 
R1=RAVE0T 
R2=RNDT 
R3=RHDT 
SG2=SG 
GO TO 38 

37 IFIKOUNT.NE.l) GO TO 39 
SG3=SG 
RA=R1 
R5=R2 
R6=R3 
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TABLE IV-2 (CONI' D) 


39 RAVEDT = .l<»^(SG+SG3)/2.*R0AR*0ELT4T 
RNDT=R!>*( SG+SG3) /2.*RNAVE*DELTAr 
RHDT=R6*(SG+SG3)/2.*RHAVE*0ELTAT 
R4=RAVEDT 
R5=RNDT 
R6=RHDT 
SG3=SG 

38 A8SS=( AGS-RAVEDT )*NS 

IF( ABSS.LE.O.O.OR.SG.LE.O.OI ABSS=0,0 
ABNS= ( AGS-RNDT ) *NM 

IF( ABiMS.LE.O.O.OR.SG.LE.0.0) ABNS=0.0 
IF(ORDER.LE.2) ABPS= ( LGS I-Y* C MS+NN) )*SG 
IF(0RD£R.LE.2) go TO 36 
ABPS=(LGSI-TE-Y*C4S*.M.Mn*SG 
36 PIRCRC=PI«RCSQD 

APHS=PIRCRC-AGS+RHOr 

IF( APHS.GE.PIRCRC.OR.SG.LE.O.OI APHS=PIRCRC 
APNS=PIRCRC-AGS+RNDT 
!FCK.LT.2) APHSl=APHS 
IF( APN'S.GE.PIRCRC) APNS = PIRCRC 
50 IF(MT.EO.O.O) GO TO 371 

IF(Y.LE.O.O) READ(4, 21111) LT P, DTP» THET TP , TAUEFF 


C ♦ READ IN geometry ASSOCIATED WITH TERMINATION PORTS (NOT ♦ 
C * REQUIRED IF NT=0) * 
C ♦ ♦ 

C ♦ THE FOLLOWING VARIABLES ARE OBTAINED FROM THE STATISTICAL ♦ 
C ♦ ANALYSIS PROGRAM ♦ 

C * * 
C ♦ LTP IS THE INITIAL LENGTH OF THE TERMINATION PASSAGES ♦ 
C ♦ IN INCHES ♦ 
C ♦ OTP IS THE INITIAL DIAMETER OF THE TERMINATION PASSAGE ♦ 
C ♦ IN INCHES ♦ 
C * THETTP IS THE ACUTE ANGLE BETWEEN THE AXIS OF THE PASSAGE ♦ 
C ♦ AND THE MOTOR AXIS IN DEGREES ♦ 
C ♦ TAUEFF IS THE ESTIMATED EFFECTIVE WEB THICKNESS AT THE ♦ 
C ♦ TERMINATION PORT IN INCHES ♦ 


Q «««««« 4c4c999«99 9 9999999 9«99999 9««99«9999«9« *999 99 9 94I99999999999999999 9 

IF(Y.LE.O.O) WRITE(6,606) LTP.DTP, THETTP, TAUEFF 
606 F0RMAT(20X, 'TERMINATION PORT GEOMETRY* ,/, I3X, 'LTP* • , F6.2 , / , I 3X , ' D 
ITP= • ,F5.2,/,13X, *THETTP= ' , F7.5 ,/ , I 3X , • T AUEFF= ',F6.3,//) 

THE TTP=THETTP/57. 29578 

DABT=NT93.1A1599( (DTP+2,*Y) ♦(LTP-Y/SINC THETTP) ) - ( 0TP92 .*Y ) *♦2/4.+ 

1 1 Y9DIP/2. )*(DTP/2. )*( l.-l./S INI THETTP ) ) ) 

IFIY.GE. TAUEFF) DABT=0.0 
371 IFIY.GT.O.O) GO TO 52 
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TABLE IV-2 (CONT’D) 


IF(NT.NE.O.O) go to 45 

LTP=0.0 

DTP=0.0 

45 IFIGRAIN.NE.2) GO TO 49 
LGCI=0.0 
LGNI=0.0 
DISQ0=0.0 
D0SQD=4.*RCSQD 
49 IF(GRAIN.EO. I ) LGSI=0.0 

VC I = 1. !♦ ( ANUM*DI SQD« ( LGC I + LGM II ♦ ( ANUM*DOSQD-AGS ) ^LOS I ♦NT*LTP*ANUM* 
1 OTP*OTP)+VCIT 
52 BBP=0.0 
BBS=0.0 
BBN=0.0 

IF(K.NE.O) GO TO 521 
IFIKKL.ECJ.O.AND.KKM.EQ.O ) GO TO 521 
IF( KKL.EO.O ) ABPC = A8PC’«'< BZ+AZ*CHIN) 

IF( KKM.EO.O) ABPC=AfiPC«( AZ*CHIH+BZ ) 

IF(KKL.EQ.0.AND.GRA1N.EQ.2 ) A8PS=ABPS* ( BZ+AZ*CH IM ) 

IF( KKM.EO.O. AMD.GRAI^.EO.2 ) ABPS=ABPS*( AZ*CH IH+BZ ) 

521 ABPORT=ABPT+A8PC+A8PS+OA8T*CBP 
A8SLOT=ABST+AHSC+ABSS+BOS 
ABNOZ=ABn»T+AHNC+ABNS + BBM 
IF(K.GE.2) GO TO 55555 
SUMAR=ABPORT*ABSLOT+ABNOZ 
55555 CONTINUE 

IF(K.EO.O) GO TO 99 
IF(ZW) 322,321,323 

322 IFIK.EQ.n ABPOR T = ABPORT«CHl N 
IF(K.EQ.2.AND.GRAIN.NE.2 ) ABPORT=ABPORT+SEH*LGC 

IF( K.EQ.2, AND.GRA IN.EQ.2 ) ABPORT=ABPORT+SEH« ( LGS I-Y* ( NS+NN I » 

GO TO 33333 

323 IF(K.EO.l) ABPORT=AHPORT*CHIH 
IF(K.EQ.2.AN0.GRAIN.NE.2 ) ABPORT=A8PORT+SEN*LGC 

IF( K.EQ.2. AND.GRA IN. E 0.2) A8P0RT=ABP0RT ♦SEN*( LGS I-Y*( NS+NN) ) 

33333 IF(K.EQ.l) ABMAI N=ABPORT +A8SL0T+ABN0Z 
K=K + l 

IFIK.GT.2) GO TO 69 
GO TO 2 

69 ABT0=ABP0RT4-ABSL0T+ABN0Z 
99 CONTINUE 

if(y.gt.o.o) go to 70 

ABP1=ABP0RT 

ABN1=ABN0Z 

ABS1=ABSL0T 

70 ABP2=(ABPl+ABP0RTI/2. 

ABN2=(ABNl+ABN0Z )/2. 

ABS2=(AHS1 +485101) /2. 


( 
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TABLE IV-2 (CONT'D) 


IF( INPUT. EQ.n GO TO 76 
GO TO ( 71 f 72*73,74) .ORDER 

71 APHEAD=APHS1 
APNOZ=APNT 
SG=SGH 

GO TO 75 

72 APHEAD=APHTl 
APNOZ=APNT 
SG=0.0 

IF(GRAIN.EQ.3) SG= ( SGH+SGN ) /2 . 
GO TO 75 

73 APHEAD=APHTl 
APNOZ=APNS 
SG=SGN 

GO TO 75 

74 APHEAD=APHS1 
APNQZ=APNS 
SG=SGN 

GO TO 75 

76 APHEAD=APHT 
APNQZ = APiMT 

75 Y=YB 

DIFF=SUMA8~SUM2 
DAOY=DIFF/DELY 
A6P1=A8P0RT 
ABN1 = AI3NCJZ 
ABSl=ABSLOT 

IF( ZW.GE.0.0) GO TO 77 
A8Ml=ABMAIN 
ABMAIN=A6T0 
ABTO=ABMl 

77 RETURN 
END 
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TABLE IV-2 (CONI' D) 


SUBROUTINE OUTPUT 


C ♦ SUBROUTINE OUTPUT CALCULATES BASIC PERFORMANCE PARAMETERS * 
C ♦ AND PRINTS THEM OUT . ♦ 
C ♦ (WEIGHT CALCULATIONS ARE PERFORMED IN THE MAIN PROGRAM) ♦ 
C ♦ T IS THE TIME IN SECS ♦ 
C ♦ Y IS THE DISTANCE BURNED IN INCHES ♦ 
C ♦ SUMAB IS THE TOTAL BURNING AREA OF PROPELLANT IN IN**2 ♦ 
C ♦ (IF ANY) ♦ 
C ♦ F IS THE THRUST IN LBS * 


REAL MDIS.ME, IT0T,M2.MDBAR, ITPLOT, ITPLTl, IDIFF,IADIFF 
COMMON/CGNSTl/ZW, AE« AT , THETA, ALFAN 
C0MM0N/CQNST2/CAPGAM,ME,B0T,ZETAF,TB,HB,GAM 
COMMON/CONST5/KPLT 

COMMON/VARI Al/T,OELY,DELTAT ,P0N02,PHEAD,RN0Z,RHEAD, SUMAB, PHMAX 
COMMON/VARI A3/ITOT, I T VAC , JROCK , I SP , I SP VAC , MD I S , MNOZ , SG, SUMMT 
COMMON/VARI A5/ABMA IN, AB TO, SUMDY, VC I , VC 
COMMON/VARI A7/Y, F 

COMMON/PAIRl/TWl ,TW2,DTW,FWl,FW2,DFWl,DFW2,OFW,TMAXQ,DFM0, 
2FDIFF,TOIFF,NX 
C0MM0N/PL0TT/IP0,N0UM,NP,10P 
COMMON/ PLOT2/NUMPLT 
COMMON/OUTl/FDIF IG, TDIFIG,OIT,ADIT 
COMM0N/OUT2/DFl00K,T100K 
COMMON/OATA2/ [COUNT 

DIMENSION T0FPLT( AOO) ,T0TPLT(400) ,TOFPLl (AOO) ,TOTPLl (400) 

DIMENSION F PLQ n 400 ) ,FPLOTl( 400) .ITPLOT (400), ITPLTl (400), 
2TPLOT(400),rPL(JTl(400) 

DIMENSION FOIFF(400),IDIFF(400),TDIFF(400),IADIFF(400) 

DIMENSION NUMPLT(5) 

DATA G/32.1725/ 

IF(Y.LE.O.O) NT0=0 

IFINDUM.EQ. 1 ) GO TO 2 

NP=NP+1 

YSFT=0.0 

YB=Y 

IFIY.LE.0.0) M2=MDIS 
MDBAR=(M2+MDI S)/2. 

SUMMT=SUMMT+MDBAR*DELTAT 

PRES=( I .+BOT/2.*ME«ME)**(-GAM/BOT) 

ALT=HBO(T/TfiM‘*17./3. ) 

PATM= 14. 696/EXP (0.43 I 03E-04*ALT) 

IF(MOIS.LE.O.O.OR.PONOZ.LE.O.O) GO TO 45 

CF=CAPGAM*S0RT(2.*GAM/8OT#( l.-PRES**(BOT/GAM) )) +AE/ AT* ( PRES-PATM/P 
lONOZ) 

F=ZETAF*COS(THETA)*PONOZ*AT*( ( 1 . +COS ( ALFAN ) ) /2. *CF+ ( 1 .-COS( ALFAN ) ) 
1/2.*AE/AT*(PRES-PATM/P0N0Z) ) 
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IF(F.LE.O.O) F=0.0 
IF(V.LE.O.O) F2=F 
FBAR=(F+F2 )/2. 

ITOT=ITOT+FBAR«DELTAT 

M2=MDIS 

F2 = F 

IF(PHEAD.GT.PHMAX) PHMAX=PHEAD 
GO TO A7 
A5 F=0.0 

47 WRITE(6,1) T,YB,PONOZ,PHEAO*SUM4B,F* ITOT 

1 F0RMAT(//,13X,*TIME= • . F7. 3 . 12X , • Y= • . F6. 3 . / , 13X , 

2*P0N0Z= •,IPE1U4,* PHEAD= *,IPE11.4,» SUMAB= *«1PE1U4« 
3* F= '.IPEU.A,* ITOT= ',IPE11.4) 

IF(IPO.EQ.O) RETURN 
TPLOT(NP)=T 
FPLOTlNP)=F 
ITPLOT(NP)=ITOT 

IF(TPL0T(NP).LT.100.) GO TO 50 
NTO=NTO+l 
TOTPLT(NTO)=T 
TOFPLT(NTO)=F 
50 RETURN 

2 NP=NP+2 
NT0=NT0+2 
I0P = l 

IF(KPLT-l) 4000,4000,4001 

4000 NP2=NP-2 
NT02 = NTC)-2 
WRITE (I, 4002) NP2 

4002 FORMAT! 14) 

WRITE (1,4003) (FPLUT! I ) , ITPLOT! I ),TPLOT ( I ) , I=l,NP2) 

4003 F0RMAT(E12.5) 

WRITE! 1,4002) NT02 

WRITE (1,4003) (TOFPLT! I ) , TOTPLT (I ) , I = I , NT02 ) 

GO TO 1004 

4001 REWIND I 

IF!IP0.NE.3) WRITE!6,9998) 

9998 FORMAT!//, 20X, 'TABULATED IMBALANCE DATA',/, 

213X,' TIME *,10X,' FDIFF »,10X,» IDIFF *, 
210X, • I ADIFF • ) 

READ!!, 4002) NP21 

READ! 1,4003) (FPLOTK I ) , ITPLTU U.TPLOTK I), I = 1,NP2U 
READ! 1,4002) NTOl 

READ! 1,4003) (TOFPLl! I ) , TOTPL I !I ) , I = 1 , NTOl ) 

NPl=NP21+2 

.► IF! IPO. EG. 2) GO TO 8888 

► IFUCOUNT.E0.2 ) YSFT=1.5 

► IFINUMPLT ! I ).NE.O) GO TO 7001 
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► CALL PLOTITIFPLOTI ,TPL0T1*NPI ,FPLOT,TPLQT,NPt ‘FHRUST 

► 2‘TIME ( StCS) -11,0. 0,400000. 0,0. 0,10. 0*9,0,YSFTI 
►7001 XSFT=18.0 

► lF(r\IUMPLT(l)..NE.O) XSFT=9.0 

► NTl=fjroi + 2 

► IFCNUMPLT (2).NE.O) GO TO 7002 

► IF(NUMPLT ( 1 ) .tO.O) YSFT = 0.0 

► CALL PLOriT(TOFPLl,raTPLl,'4n,TOFPLT,TOTPLT,.MTO,«TH«UST (LBS)*«12, 

► 2'TIME (SECS) -11,0. 0,400000. 0,100. 0,2.0, XSFTjYSFT) . 

► 7002 XSFT=18.0 

► IF(.MUMPLT ( I ). 'JE.0.AMD.NUKPLT(2).NE.0) XSFT=9.0 

► 8888 COiMTLMUE 

IF(NPl-NP) 2000,2000,2001 
2000 NX = .NP-2 
NY=NPl-2 

CALL MTLRPt TPL0T,FPL0T,MX,TPL0T1,FPL0T1,NY,FDIFF,0 ) 

CALL INTERP(TPLOT,ITPLOT,NX,TPLOTl,ITPLri,NY,IDIFF, 1) 

TOIFIG=TPLOr( 1) 

FDIFIG=A9S(FDIFF( I ) ) 

DO 3000 J=2,NX 

IF(rPLaT(J).GT..02<'TB) GO TO 3001 
IFIABS{FDIFF(J)).LT.ABS(FDIFF(J-1)) ) GO TO 3000 
FDIFIG=ABS(FDIFF( J) ) 

TOIFIG=TPLOT( J) 


3000 CONTINUE 

3001 CONTINUE 

DO 2004 1=1, NX 
2004 TDIFFIi)=TPLOT(I) 


2003 


DUM1=0.0 

lADIFFd )=ABS(FDIFF(1 ) /2 . ) ♦TPLOT ( 1 ) 

DO 2003 I =2, NX 

FBARI= (FDIFF( n+F0lFF(I-l))/2. 

DUMl=ABS (FBARI )♦ ( TPLOT ( I )-TPLOT ( I-l ) ) 

IADIFF( I )=I ADIFF ( I- I ) +0UM1 

IF( IP0.NF.3) WRITE (6,9999) (TPLOK I ) ,FDIFF( I ) , I D I FF ( I ) , I ADI FF ( I ) , 


21=1, NX) 

9999 FORMAT ( 1 3X , 1 PE 1 1 . 4 , I OX , I PE 1 1 . 4 , 1 OX , 1 PE 1 1 . 4 , lOX , I PEI 1 . 4 ) 
TI=AMIN1 ( TWl, TW2 ) 

CALL INTRPl ( IDIFF,TPLOT,NX,TI ,DIT1,0) 

DIT=IDIFF(NX)-DIT1 

CALL INTRPl ( I AOI FF, TPLOT,NX,T I , AD IT 1,0) 

ADI T=IADIFF (NX )-ADI T1 

CALL INTRPl ( FDIFF ,TPLOT ,NX,TMAX0,DFMQ,0 ) 

CALL INTRPl ( FDIFF, T PLOT, NX, TWl, DFW 1,0) 

CALL INTRPl ( FDIFF , TPLOT, NX, TW2,DFW2,0) 

CALL INTRPl (TPLOT, F PLOT, NX, lOOOOO. ,T100K2,1 ) 

CALL INTRPl (TPLOT I, FPLOT I, NY, 1 OCOOO . , T 1 OOK 1 , 1 ) 

T100K = AM,AXl( TlOOKl, T100K2) 
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CALL INTRPl (FDIFF,TPL0T,NX,T100K,DF100K,0) 

► IF(IPO.E0.2) GO TO R887 

► CALL SCALE( IADIFF,8.0,MX, I ) 

► YSCAL1=-ABS{ 8.0»I ADIFF ( NX+2 ) ) 

► YSCAL2=ABS(2.0«IADlFF(NX+2) ) 

► NX=NX+2 

► IF(NUKPLT(3).NE.O) GO TO 7003 

► IF(NUMPLT( 1 ).E0.0.0R.NUMPLT(2).E0.0) YSFT=0.0 

► CALL PLOri (TPLOT.FDIFF, NX, 'THRUST IMBALANCE (L8S)*f22, 

► 2'TIME (SECS)', -ll,-A00000. ,100000. ,0.0, 26.0, A. 0,XSFT,YSFT1 

►7003 XSFT=9.0 

► IFINUMPLT (3).NE.0) XSFT=18.0 

► IF(NUMPLT(A).NE.O) GO TO 7004 

► I F( NUMPLT ( 1 ) .EO.O.OR.NUMPLT ( 2 ) . EQ. 0. OR . NUMPLT I 3 ) . EQ .0 ) YSFT = 0.0 

► CALL PLOTl (TPLOT, IDfFF,, NX, • IMPULSE IMBALANCE I LB-SECS )', 27, 

► 2'TIME (SECS) ' ,-ll,YSCALl, YSCAL2, 0.0, 26. 0,4.0, XSFT,YSFT) 

►7004 XSFT=9.0 

► IF(NUHPLT(3).NE.0.ANU.NUMPLT(4).NE.0) XSFT=18.0 

► IFINUMPLT (5).NE.O) GO TO 7005 

► IF (NUMPLT ( 1 ). EO.O.OR.NUMPLT (2 ). EO.O.OR.NUMPLT ( 3). EO.O.OR.NUMPLT (4) 

► 2.EQ.0) YSFT=O.G 

► CALL PLOTl (TPLOT, lADIFF, NX, 'ABS. IMPULSE IMBALANCE ( LB-SECS 32, 

► 2' TIME (SECS) ' ,-ll, lADIFFlNX-l ), IADIFF(NX),0.0,26.0,0.0,XSFT,YSFT) 

► 7005 CONTINUE 

► NX=NX-2 

► 8887 CONTINUE 

GO TO 1004 


2001 NX=NPl-2 
NY=NP-2 

CALL INTERP (TPLOT 1 , FPLOT L , NX , TPLOT, FPLOT , NY , EDI FF , 0 ) 
CALL INTERP (TPLOT I , I T PLT 1 , NX , TPLOT, ITPLOT,NY,IDIFF,l) 
TDIFIG=TPLOTU 1) 

FDIFIG = AHS( FOIFF( I ) ) 

DO 3002 J=2,NX ; 

IF(TPLOT( J).GT/.02«TB) GO TO 3003 
IF(ABS(FDIFF( J) ).LT.A8S(FDIFF(J-1))) GO TO 3002 
FDI FIG=ABS ( FOI FF ( J ) ) 

FDIFIG=Fl;IFF( J) 

TDIFIG=TPLari{ J) 

3002 CONTINUE 

3003 CONTINUE 

00 2005 1=1, NX 
2005 TDIFFI I )=TPL0T1( I ) 

DUM1=0.0 

lADIFFd )=A13S( FDIFFd )/2.)*TPL0Tl( 1) 

DO 2002 1=2, NX 

FBARI = (FniFF( I ) + FDIFF( I-l) )/2. 

DUMl = ABS (F BARI )<=( TPLOT 1( I ) -TPLOT I ( I-l ) ) 
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> IADIFF( I ) = IADIFF( I-l )+DU«l 

IF( IPO.Mt.3) WRITE (6,9999) (TPLGTK I) , FO I FF ( I) , 1 0 IF F ( I ) , I AO ! FF ( I) , 
2 1=1, MX) 

TI=AMI.a ( rWl,TW2 ) 

CALL I URPK IOIFF,TPLnTl,^)X,TI,OITl,0) 

DIT=IDIFF(NX)-DIT1 

CALL I.aRPl ( 1 ADIFF,TPLOTl, )X,TI ,ADITI,0) 

A0IT=IADIFF(MX)-ADITI 

CALL INTRPl(FDIFF,TPLOTl,NX,TMAXQ,DFM0,0) 

CALL I'JTRPl(Ff)IFF,T'’LOTl,.VX,TWl,OFWl,0) 

CALL INTRPKFDIFF.TPLOTl ,NX,TW2,DFW2,0) 

CALL INTRPK rPLnr,FPLOT,NX, lOOOOO. , T100K2, I ) 

CALL IVr<Pl(TPLOTI,FPL()Tl,Wy,10COOO.,T100Kl,l) 

TlOOK=AMAXl( TlOOKl , T100K2 ) 

CALL INTRPKFDIFFfTPLOn ,NX,T100K,DF100K,0) 

1F(IP0.EQ.2) GU TO 100^ 

CALL SCALE! IADIFF,R.O,NX, I ) 

YSCALI=-A8S(8.0*IADIFF(NX+2) ) 

YSCAL2 = A8S(2.0*IADIFF(i'JX+2) ) 

NX = fMX + 2 

IFCNUMPLT (3).;'4E.O) GO TO 7006 

IF(NUMPLT{ l).EQ.O.OR.NUMPLT(2).eQ.O) YSFT=0.0 

CALL PLOrUTPLOTl,FOIFF, NX, ‘THRUST IMBALANCE (LBS)', 22, 

2 'TIME (SECS) • ,-11,-400000., 100000. ,0.0, 26.0, A. 0,XSFT,YSFT) 
j XSFT=9.0 

IF(NUMPLT(3).NE.O) XSFT=18.0 
IF(NUMPLT(4).NE.O) GO TO 7007 

IF( NUMPLT ( I ) .EO.O.fJR.NUMPLT (2 ) .EG.O.OR.NUMPLT ( 3) .EO.O ) YSFT = 0.0 
CALL PLOTl ( TPLOTl , I 01 FF,-\X, • I MPULSE IMBAL ANCE ( L B-S EC S ) ' , 27 , 

2 'TIME ( SECS ) • ,-I 1, YSCAL1,YSCAL2,0.0,26.0,4.0,XSFT,YSFT) 

' XSFT=9.0 

IF(NUMPLT(3) .NE.O. AND.NUMPLT(4) .NE.O) XSFT= 18.0 
IF(NUMPLT(5).NE.O) GO TO 7008 

IF(NUMPLT( 1 ) .EG.O.OR.NUMPLT (2 ) .EG.O.OR.NUMPLT I 3) .EG.O.OR.NUMPLT (4) 
2.EQ.0) YSFT=0.0 

CALL PL0T1(TPL0T1,IADIFF,NX,‘ ABS. IMPULSE IMBALANCE ( LB-SECS )', 32 , 
2* TIME ( SECS ) • ,-l 1, I ADIFF (NX-1 ) , I AD I FF ( NX ) , 0 . 0 , 26. 0 , 0 . 0 , X S FT , YSF T ) 

! CONTINUE 
NX=NX-2 
h CONTINUE 


RETURN 

END 
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► SUBROUTINE PLOT I T ( Y I , XI , NPl , Y2 , X2 . NP2 ♦ YHDR , MY , XHDR t NX , SY 1 , S Y2 , 

► 2SXl,SX2,XSFT,YSFT) 

► DI MENS I 0 J XHDR ( 8 ) , YHDR {8),Xl(NPl),Yl(NPl)tX2( MP2 ) , Y2 ( NP2 ) 

► Nl=NPl-2 

► NSl=;NPl-l 

► N2=NP2-2 

► NS2=NP2-1 

► XI(NS1)=SX1 

► X1{NP1)=SX2 

► X2(NS2)=SX1 

► X2(NP2)=SX2 

► Yl{rgsiJ = SYl 

► Y1(NP1)=SY2 

► Y2(NS2)=SY1 

► Y2(NP2)=SY2 

► CALL PLQT( XSFT,YSFT,-3) 

► CALL AXI S (0.0, 0.0, YHDR.NY, 8. 0,90.0, SYl, SY2) 

► CALL AXIS(0. 0,0.0, XHDR, NX, lA. 0,0.0, SX1,SX2) 

► CALL LINE(Xl,Yl,Nl, 1,0, I ) 

► CALL LINE ( X2,Y2,N2, 1,0,2) 

► NPLOT^NPLOT+l 

► RETURN 

► END 
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SUBROUTINE OVAL 
REAL Ml.Nl 

C0MM0N/CUNST4/DEL0I , DOt D 1 , ZC, XT f ZO 
COMMON/VARI A7/Y 

commqn/ovalm/z.zq 

COMMON/OVALM2/KK I , I I 

COMMON/QVALA/CHIH,CHlN»SEN,SEH,CHlNH, AZ.BZ 
COMMON/OV ALB/CHINN, CM I NAV.SENN 

COMMON/OVALC/RCNDCN.KONOCH.RONOGNf RONDGH,EXN,EYN,EXH,EYH, 
2ALPHAN, ALPHAH 
DATA PI/3.1A159/ 

KKI=KKI+1 

IF(KKI.GT.l) GO TO 8 

AGN=( KUNDGN+SQKT ( RONOGN^^^Z + DI^^Z ) ) /2. 

BGN=AGN-KGNDGN 

AGH=(RONDGH+SQRT ( R0NDGH*»2+D I **2 ) )/2. 

BGH=AGH-RONDGH 

DTH=2.*PI/I I 

KKM=0 

KKL=0 

KKJ = 0 

KKXT=0 

KKP = 0 

AX = 0. 

AZ=0. 

BZ = 0. 

ACN=(RONI)CN+(RCNDCN<'<‘2+(DO-ZC )<'*2)«*.5)/2. 

8CN=ACN-R0NDCN 

ACH=(R0N0CH+{RCNDCH*<'2*(D0 + ZC)«'«=2)**.5)/2. 

BCH=ACH-RCNDCH 

A1N=(C0S ( ALPHAN) ) *♦2* ( ACN/6CN ) ( S I N { ALPHAN ) ) **2 
A1 H= I COS ( ALPHAH) ) **2 +( ACH/ BCH ) ««2« ( S I N ( ALPHAH ) )**2 
BIN=( ( ACN/8CN) **2-1 . ) ♦ S I N ( 2 . * ALPHAN ) 

BIH=( {ACH/BCHM=*2-1. )*SIN{2.*ALPHAH) 

CIN = 2.»( EXN*COS( ALPHAN)- (ACN/BCN H‘*2»EYN*S INC ALPHAN ) ) 
ClH=2.«{EXH*CnS( ALPHAH) -{ACH/BCH)*»2#EYH«S INC ALPHAH) ) 

D1N = 2.«C C ACN/BCN) ♦*2*EYN*COSC ALPHAN )-EXN«S IN CAL PHAN ) ) 
DIH=2.*C C ACH/BCH ) 2«E YH«=COS C. AL PHAH ) -E XH'f'S I N C ALPHAH ) ) 

ElN = ( SIMC ALPHA')) ) **2+ C ACN/liCN ) *<'■2* C COS C ALPHAN ) )*«=2 
EIH=C SIN C ALPHAH) ) C ACH/BCH ) C COS C ALPHAH) )**Z 
F1iN = ACN*^Z-£XN*4'2-I C ACN/BCN ) *EYN ) <=*2 
F1H=ACH**Z-EXH«<‘2-C C ACH/ RCH ) *EYH) ♦♦Z 
SENNO=PI* CDO-ZC) 

SENO=SENNO 
SEHO=PIO C DO + ZC ) 

8 KK=0 
YO=Y 

3 IFCKK.EQ.l) Y=YO+ZQ/2. 
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IF(KK.EO.l) GO TO 5 
2 IF(KK.EQ.2) Y=Y0-Z0/2. 

IF(KK.EQ.2) go to 6 

IF(KK.E0.0.AN0.XT.GT.0.) Y = Y0*-XT + 2Q/2. 

IF{ KK.EQ.O.ArJD.XT.GT.O. ) GO TO 7 

KK=l 

GO TO 3 

5 THETA=0.0 
SUM0=0. 

DO 12 1 = 1, II 
THETA = THETA + [)TH 

M1 = A1,M«{ cost THETA ) J **2 + fl IN* S I N ( THETA ) #005 I THETA ) ♦ 

2EliM«‘I SINI THETA) )«*2 
N1=CIN«'C0S ( THETA )+DlN»SlN( THETA ) 
RC=(-Nl+S0RT(iNi*«24^A.<'Ml«FlN) )/t2.*Ml ) 

IFIRC.lt. 0. ) RC=(-N1-SQRTINI*«'2+A.«M1*F1N) )/l2.*Ml) 

RG=SQRT (!./(( COS ( THET A ) / ( AG J+Y ) ) ♦♦2 + ( S INI THETA) /( BGN + Y ) l**2 I I 
IFIRG.GE. RC) KKM=1 
IF(RG.GE.RC) RG=0. 

SUMO=SUMO*RG«DTH 

12 CONTINUE 
IF(KKM.EO.l) SEN=SUHO 
IFl SUMQ.LE.O. ) SEN = 0. 

IF(KKM.EO.O) GO TO 9 
CHIN=SEN/SE NO 
CHINAV=( I .+CHINN) /2. 

IF! XT.LE. 0.0) CHINA V= 1.0 

IF(KKJ.EO.l) CHINAV=( l.-AX)*CHIN+AX*CHINN 
CHINH=( l.+CHlN)/2. 

9 KK = 2 

IF( Z.GE.O.O.AND.KKM.EQ.O) GO TO 62 
GO TO 2 

6 THETA=0.0 
SUM0=0.0 

DO 13 1=1,11 . 

THETA=THETA+DTH 

Ml=AlH*(CnS( THETA ) ) «*2*B IH* S I N ( THE T A ) ^COS ( THE TA ) ♦ 
2E1H*{SINITHETA) )^=!'2 
N1=C1H*C(JS( THETA ) +D 1 H*S I N ( THE T A ) 

RC=(-N1 + SQRT(N1*«2 + 4.*M1*FIH) )/(2.*MI ) 

IFIRC.LT.O. ) RC=(-Nl-SQRT (Nl*<=2 + A.*MI«FIH) )/{ 2.*Ml) 

RG=SQRT (!./(( COS I T HET A ) / ( AGH + Y ) ) **2 ♦ ( S I N ( THET A ) / ( BGH+Y ) )**2)) 
IF(RG.GE.RC) KKL=1 
IF(RG.GE.RC) RG=0. 

SUMO=SUMU+RG*DTH 

13 CONTINUE 
IFlKKL.EO.l) SEH=SUMO 
IFISUMO.LE.O. ) SEH=0. 
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CHIH=SEH/SEHO 
IF(KKL.EQ.O) CHIH=1.0 
CHINH={ l.+CHIH)/2. 

IFCKKM.EQ.l) CHINH=(SEM+SEH)/(SENO+SEHO) 

GO TO 62 
7 THETA=0.0 
SUM0=0. 

DO 1 1 I = 1 , 1 1 
THETA=THETA+OTH 

Ml=AlN*(COS( THETA) )*»2+8l I N( THETA )*COS( THETA )♦ 
aElN^'C SIM( THETA) ) ^«2 
Nl=ClN«COS(THETA )+Dl,M*S IN(THETA) 

RC=(-N1 + SQRT( Nl#*2 + A. «M1«=FIN) )/(2.*Ml) 

IF(RC.LT.O. ) «C=(-iN/l-SQRTl 'a**2+4.*Ml»FlN) )/(2.*Ml) 

RG=SQRT( 1 ./( ( COS ( THETA)/ (AGN+Y) ) «*2+ ( S IW ( THETA )/{ BGN+Y ) )*«2) ) 
IF(RG.GE.RC) KKJ=1 
IF(RG.GE.RC) RG=0. 

SUMO=SUMG+RG*OTH 
11 CONTINUE 

IF(KKJ.EO.l) SENN=SUMO 

IF( SUMO.LE.O. ) SENN=0.0 

IF(KKJ.EQ.O) GO TO 9 

CHINN=SENN/SENNO 

KKXT=KKXT+1 

IF(KKXT.EQ.l) YXIP=Y 

AX = ( Y-YXIP)/ (XT4^00/2.-DI/2.-YXIP) 

IFIAX.LE.O.) AX=0. 

IF( AX.GE.l.O) AX=l.O 
CHINAV=1 .-AX+AX*CHIMN 
KK=l 

IF( AX. LE. 0.5. AND. XT. GE. 0.02097*00) GO TO 9 
GO TO 3 

62 Y=YO 

IF(KKL.EQ.O.AND.KKM.EQ.O) GO TO 63 
KKP=KKP+1 

IF(KKP.EQ.l) YZIP=Y 

AZ=(Y-YZ IP)/( ABS(Z ) /2 . +D0/2 ,-0 I /2.- YZ I P ) 

IFIAZ.LE.O.) AZ=0. 

BZ=l.-AZ 

63 CONTINUE 
RETURN 
END 
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subroutiinie setup 

INTEGER TEMPCn,CnOE 
REAL TI200) 

REAL ANSI 60 I 

REAL TEMPAI 10» ,CnNST(60l 
INTEGER ORDER ( 601 tCNSTNH 
OEAL FXARFAI2, 1001 
PEAL XXI 105) . YY( 105) 

REAL PSEU00(105) 

REAL X(^0, 105) ,Y( 105) ,FX(40, 105 ) 

C « IF the dimension OF X AND FX APE CHANGED M AND N SHOULD * 

c * also be reset ♦ 

REAL MODEf MEAN.Mlf M2,K,INC 
r.ONMON/SEFD/lX 
inptnm=o 
r.NSTNM = 0 
N = 105 
NI=100 
NSI=10 
M = AO 
MM = 0 
MI I=NI+1 
MSIUMSI+l 
READ! 5, 100) IX 
^0 CONTINUE 

PEAD(5,106) NAMl,NAM2 

106 F0RMAT(2A4) 

READ! 5, l02)Cnoe,Xl , X2f X3,X4f X5f X6tX7 

WRITE (6, 107) NAMl ,NAM2,Cflt)E,Xl,X2,X3f X4,X5,X6,X7 

107 FORMAT! IX,2A^^,5X, I2»5X,7( 1PEU.4.3X) ) 

IF(C0DE.EQ.90) GU TO 399 
INPTNM=INPTNM+1 
TF(C0DF.E0.60)GC TO 356 

ORDER ( INPTNM)=MM 
TEMDCD=COOE/10 

GO TO (31,32»33,3tt35),TEMPCO 
31 CONTINUE 
N0I=X4 
NOIl = NOH-l 
X(M.M, I ) = X2 
DO 311 1=2, NOI 
X(M,M, I) = X(mm,i_ij4.x3 
311 CONTINUE 

DO 312 1=1, NO! 

Y(I)=0. 
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312 CONTINUE 
H=X3 

STAPTR=X2-X3/2. 

SUM=0. 

NQV=X1 

NOC=( Xl+9. ) /lO. 

00 313 JJ=1,N0C 

READ! 5, 104) (TEMPAI I) ,1=1,10) 

WRITE(6,109) ( TFMPAI IJ , 1=1, 10) 

00 314 J = l, 10 

IF( JJ*104-J.GT.NOV)GO to 317 
no 315 1=1, NOl 

IF (TEMPA(J) .LT.XINM, I )+X3/2. )G0 TO 316 

315 CONTINUE 
GO TO 314 

316 CONTINUE 

Y( 1 )=YII )«-l. 

SUM=SUM*-1. 

314 CONTINUE - 

313 CONTINUF 

317 CONTINUE 
IFICOOF.EQ.IDGO TO 99 
FXIMM, 1 )=0. 

00 318 I =2, NO II 

FXIMM, 1 )=FX(MM, I-l )+-Y ( I-D/SUM 

318 CONTINUE 
GO TO 30 

32 CONTINUE 
N0I=X1 
X(MM, l)=X2 
00 220 1=2, NOI 
X( MM, I ) = X( VM, I-D + X3 

220 CONTINUE 

READ(5,104) (Y( I ),I=l,NOI ) 

WRITE(6,109) ( Y( I ) , I=1,NCI ) 

H = X3 

STARTR=X2-X3/2 . 

IF(CnOE.Eg.21)GO TO 99 
SUM=0. 

00 222 1=1, NOI 
SUM=SUM+Y( I ) 

222 CONTINUE 
NOI1=NOI+1 
FXIMM, 1 )=0. 

DO 221 I=2,N0Il 

FXIMM, I ) =FX (mm, I-l )+YU ) /SUM 

221 CONTINUE 
GO TO 30 
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33 CONTINUF 
MEAN=X1 
S2 = X1 
U2=X2 

U3 = X3 
U4 = XA 
H = X5 

STARTR=X6 
SUMX=X7 
GO TO 331 

34 CONTINUE 
NOI = Xl 

Xf MM, 1 )=X2 
DO 341 1=2, NOI 
X|MM, I J=X(MM, I-n*X3 
341 CONTINUF 

REAO(5,104MFXIMM, II, I = l,NOn 
WRITE(6,109I ( FX( MM, I I , I = l,NOI I 
109 F0RMAT(5X, IPIOEII.41 
Go TO 30 

35 CONTINUF 
C0DE=Cn0E-50 

GO TO(35l, 352, 353, 354, 3551 , CODE 
351 CONTINUE 
MEAN=Xl 
SIGMA=X2 

IF(X6.F0.0. IX6 = mFAN“3.<‘SIGMA 
IFIXT.EO.O. )X7*MFAN+3.*SIGMA 
X0 = X6 
XN = X7 

1351 CONTINUE 

H=(XN-XOI /FLOAT'NI I 
D=H/FLnATIN5I I 
X(MM,l)=XO 

INC=CXN-XO)/FlOAT(NII 
00 201 1=2, Mil 
XIMM, I )=X(MM, I-ll+H 

201 CONTINUF 

DO 202 J=2,NI1 
TC U = X(MM, J-l) 

00 203 KK=2,NSI1 
T< KK I =Ti KK-l l+D 

203 CONTINUF 

DO 204 L=l,NSIl 

Y(LI=( l./( S0RT(6.2832)«SIGMAI l♦^EXPC-•5*^T(LI-MEANI/SIGMA)♦♦2) ) 

204 CONTINUE 

CALL CAPEAI Y,FX,M,N,MM,NSI ,J,D| 

202 CONTINUE 
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DO 205 I=2,NI1 

FX(MM, n=FX(«M, I)/FX(MM,NI II 
205 CONTINUF 
GO TO 30 

352 CONTINUE 

INC=( X2-Xl)/FLnAT(Nn 
XIMM, 1)=X1 
DO 3521 1=2, Nil 
X{MM, I ) = X(MM, I-U MNC 

3521 continue 

INC=1 ./FLOATINI ) 

FXIMM, 1)=0. 

DO 3522 1=2, NIL 
FX{MM,I)=EX(mm,I-1)+ImC 

3522 CONTINUE 
GO TO 30 

353 CONTINUE 
MEAN=X1 
SIGMA=X2 
XO=MEAN 

IFIXT.EO.O. ) X7=MEAN*3.*S IGMA 
XN = X7 

GO TO 1351 

354 CONTINUE 

355 continue 
GO TO 30 

356 continue 
CNSTNM=CNSTNM+1 

ORDER { INPTNM) = IOO+CNSTNM 
CONST(CNSTNM) =X1 
GO TO 30 
99 NEAN=0. 

SUMX=0. 

S1=0. 

52 = 0. 

53 = 0. 

54 = 0. 

S5=0. 

DO 200 L=1,N0I 

I=N0I-L+1 

SUMX=SUMX4-YIL) 

S1=S1+Y( I ) 

52=S2+S1 
S3=S3+S2 
S4=S4+S3 
S5=S5+S4 
200 CONTINUE 

MEAN=S2/SUMX 
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S2=S2/SUMX 
S3=S3/SUMX 
S^ = S4/S'JMX 
S5 = <?5/SUMX 
U2 = 2.*S3-S2*( 1. + S2 ) 

U3 = 6,*S<f-3,*a2*( 1.+S2 )-<;2«( U + S2) »( 2.+S2 ) 

U-V=24.*S5-2.*U3*(2.*( 1. + S2 Ul. I-U2<=(6.*( I . +S 2 ) * ( 2 . + S2 )- 1 . ) 

9 -S2*l l.*S2)*(2.+S2)*(3. + S2) 

331 CONTINUE 

B1=U3««2/U2«*3 

B2=U^/U2**2 

K=(Bl*(32 + 3. ) *<'21 /( 4.* (2.>» 02-3. *31-6. I *(^.*82-3. *01 1 I 
IF(K)l,98,94 

I P=(6.*(B2-0l-l.l )/(6.f3.*Ol-2.*02» 

COM=01<'(R+2.) **2 + l6.*(«+l. J 
AlA2=.5*S0PTt'12)*SQ^»T(C(JM) 

CON12=P«(R*2. ) *SORT(0l/CnM» 

IF(U3.LT,0, JC0Ml2=-CaMl2 
M2 = .5*(R-2.<-C0M12) 

MI = .5*(R-2.-C:JM12I 

Y0=(SUMX/A1/V2) l**Ml>!'vi2«<'M2l/CMl*M2l ♦* (M1 + M2)*GAMMA( MI+M2 + 2. ) / 

9(GAVMA(Ml + l . ) *GAMMA(M2«>1. n 
A2 = AlA2/(Ml/M2 + l. ) 

AI=A1A2-A2 

MO()F = MEAM-,5*J3/U2*(R*2. I /(R-2. » 

MOr)E = MOOt*H + STAPTP 
]NC=AlA2/FLnAT(MJ 
X( MN, I ) =MOOF*-( -A1 ) *H 
X I WM, NT I ) =Nnf)F +A2*H 
H= ( X{ MM»NIII-X (MM, U I /FLOAT (Nl I 
X( mm, 2 ) = STARTR 
DO 706 1=3, Nl 
XCMM, I ) = X(MM, l-n+H 
706 CONTINUE 

PSEUDOt I )=-Al 
PSEUDOINI I/) =A2 
H=A1A2/')JI ■ 

00 701 1=2, N I 
PSEUDO! I )=PSEUOn( I-l )*H 
701 CONTINUE 

D=H/FLOAT(NSI ) 

00 702 J=2,NI1 
T( l)=PSEUOO( J-1 ) 

00 703 KK=2,NSI I 
T(KK)=T(KK-IM-D 
703 CONTINUE 

00 704 L=l,NSI I 

Y( L) =Y0*( I . 4-T( L )/ Al I **Ml*( l.-TILJ /A2 »**M2 
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704 rONTlNUE 

C&LL CAqEA(Y,FXtMtN,MM,NSI,J,D) 

702 CONTINUE 

00 705 I =2, Nil 

FX(MM,1)=FX(MM,I )/FX(MM,NIl) 

705 CONTINUE 
GO TO 30 

94 IFIK-l )4,96,6 
4 CONTINUE 

R= (6.«(Q2-Bl-l. ) )/ (2.*B2-3.*Bl-6 , ) 

Ml =.5«(R«-2. ) 

COM=SORT( 16.*( P-l. )-Bl*(P-2. )**2J 
V= (-R« (R-2, 1«SQPT( Bl ) )/COM 
IF (U3.GE.0. )C,0 TO 44 
V=ABS( V> 

44 CONTINUE 

A1 = S3RT(U2/16. )«COM 

MOOE = MEAN-(U3=i=(P-2. ) )/( ( 2. +U2 1 « ( R + 2 . » ) 

THETA = ATAN( V/R » 

IF (R .LE. 10. )G0 TO 48 
Al^AK'H 

YO = SUMX/A 1«SQR T(P/6.2832 )’S'( EXPICOSI THFTA)**2/ 13.«R )-l. / 
9( U.’i'R ) -THETA* VI ) / ( COS ( THETA H *« ( R+ 1 1 
48 CONTINUE 

OR IGIN = MEAN+V*A1/R 
H=2. «OPIGlN/Fl.CAT(NI ) 

0=H/FL0AT(NSI ) 

X( mm, I )=-0R IGIN 

no 711 I = 2 , Nt 1 

X( MM, I )=X(MM, I-U+H 

711 CONTINUE 

on 712 J=2,NIl 
Tl 1)=X{MM, J-1 ) 

DO 713 KK=2,NSI1 
T( KK J =T( KK-l ) +0 

713 CONTINUE 

DO 714 L=1,NSI 1 

Y(L)=YO*( l.+T( L J’!'*2/A1**2 »**(-Ml)*EXP(-V*ATAN(TIL)/Al J ) 

714 CONTINUE 

CALL CARFAI Y, FX,m,N,MM,NSI ,J,0) 

712 CONTINUE 

DO 715 1=2, Nil 

FXIMM, I > =FX (MM, I) /FX (MM,NI 1 ) 

715 CONTINUE 

DO 716 1=1, Nil 

X{ MM, I ) =X( MM, I } +-0R IGIN 

716 CONT INUE 
GO TO 30 
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6 continue 

rMEAN=MeAN 
ME AN=HEAN-IMEAN 

R=(6.«(R2-B1-1.) )/(6.«-3.*fil-2.*82) 

COM=Hl*(^f2, ) ) 

AI=.5«S0RT(U2) «SQRT(C0M) 

IF (U3.lt. 0.)Al = -(ABS(Am 
COMl2=(P«(R+2 . ) )/2 .«SORT (Bl/COM) 

Ml=-( (R-2. ) /2.-COM12 I 
M2=(R-2. »/2.+CnMl2 

Y0 = ( Al««(Ml-M2-l . I -'GA(’m4( MI-M2-1. n«(GAMMA(Ml ) /GAMMA ( M2+ 1. ) l*SUMX 
OR IGIN = MEAN-(Al=c(vl-l.) )/(Ml-M2-2.) 

MOnE = MP&N-.5*U3/U2-(P*-2. )/(R-2. ) 

XN = AU-XN/H 
SAVEH=H 

H= (XN-Al I/FLOATINI ) 

D=H/FLOAT(NSI ) 

X(MM,U = A1 
DO 721 I=2tNIl 
X( MM, I )=X( MM, I-l ) 

721 CONTINUE 

on 722 J=2,NIl 
T( n=X(MM, J-l) 
no 723 kk=2,nsii 

TlKK)=TlKK-l)4-D 
723 CONTINUE 

DO 72't L = l ,NSI 1 

Y(L)=Y0*( T(L)-A1 )#«M2*T(L)»*(-Mn 
72A CONTINUE 

CALL CAPEA( Y,FX,m,N,MM,NSI ,J,0) 

722 CONTINUE 

DO 725 1=1, Nil 

FX(MM, n=FX(MM, I ) /FX(MM,NI n 

725 CONTINUE 

DO 726 1=1, Nil 

X(MM, I )=|X(MM, I )-Al)*SAVEH 

726 CONTINUE 
GO TO 30 

98 WRITE(6, 1031 
GO TO 399 
96 CONTINUE 

WR ITE(6, 105) 

399 CONTINUE 
RETURN 

<[!([♦#){( entry point 

ENTRY input 
REWIND A 

DO 500 J=l , INPTNM 
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IF (ORDER (J ) .GT.100)G0 TO 501 
RN0=04N0U( IX) 

DO 502 1=1, Nil 

IF (RN0.LT.FX(0RDER ( J) , I ) )G0 TO 503 

502 CONTINUE 

503 CONTINUE 

ANSI J)=X(OPOER(J), I) 

GO TO 500 
501 CONTINUE 

. ANS( JI=CONST(ORDEP (J)-IOOJ 
500 continue 

WR ITE(A, 101 ) (ANS( I ) , 1=1, INPTNMJ 
ENOFILE 4 
REWIND 4 

100 FORMAT ( I 10) 

101 format (E16. 9) 

102 FORMAT! 12, 2X, 7E10.0) 

103 FORMAT! • • , 'K=0* ) 

104 FQRMAT!10E8.0) 

105 FORMAT! • • , *K= 1. • ) 

RETURN 

END 
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SUBROUTI rslE CAREA ( Y , FX ,M, N , MM» NS I , J, D ) 
REAL FX(K,N),Y(N) 

NSI l=NSl+l 
NSIO=NSI-l 
FX(MM, 1 )=0. 

SUM=0. 

DO 201 I = 3,iMSI0,2 
SUM=SUM+4.#Y( I-l ) ♦2.«Y(I ) 

201 CONTIMUE 

AREA=D/3.*(Y( I )+SUM+Y(NSIl ) ) 

FX( MM, J) =FXC KM, J-1 ) +AREA 

RETURN 

END 

FUNCTION RANDU(IX) 

IX= 1X065541 
IFI IX)5,6,6 

5 IX=I X+2147483647+1 

6 RANDU=IX 

RANDU=RANDU0.4656613E-9 

RETURN 

END 
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TABLE IV-2 (CONT'D) 


SUBROUTINE I NTERPIXl ,Y1 »N1 ,X2,Y2,N2 ,YDIFF, ICHK ) 

DIMENSION XI (Nl ) , YL INI ) , X2 (M2 ) ,Y2(.N2 ) ,YDIFF ( NU 
DO 100 I=lfNl 
N3=N2-1 
CO 200 J=^l,.N3 

IF! I .GT.N2. AND. ICHK.EC.O ) YD I FF ( I ) =Y 1 ( I ) 

IF( I .GT.N2. AND. I CHK .EQ. I ) YD I F F ( I ) = Y I ( I) - Y2 (N2 ) 

IF( I.Gr.N2) GO TO 100 

IF! ABSI XI I I )-X2IJ) ) .GT. l.E-S) GO TO 1 
YDI FF ( I )=Y1 ( I )-Y2( J ) 

GO TO 100 

1 IF( XI ( I ) .LT. X2 ( J 1 .OR.Xl { I ) .GE.X2I J + 1 ) ) GO TO 2 

YDI FF { 1 ) = Y1 { I )-( ( Y2 (J+1 )-Y2 (J))/(X2(J + l)-X2(J)))«'(Xi(I )-X2( J ) ) 
2-Y2U) 

GO TO 100 

2 IF(X1( I )-GE.X2(J+1).AN0.J+1.LT.N2) GO TO 200 
IFIJ.EQ. 1 ) GO TO 3 

YDIFF ( I )=Y1 ( I )-( ( Y2 ( J )-Y2( J-1 ) ) / ( X2 (J ) - X2 I J- 1 ) ) ) * ( X 1 ( I ) -X2 (J- 1 ) ) 
2-Y2 ( J-1 ) 

GO TO 100 

1 YDIFF ( I )=Y1 { I )-( Y2 ( J ) / X2 ( J ) ) 4rXl ( I )-Y2( J ) 

200 CONTINUE 

100 IF( ABS(YOIFF( I ) ).LT.ABS(Y1( I )«l.E-5) ) YDIFF(I)=0.0 

IF(Nl.eQ.N2.AN0.ABS(Xl(Nl)-X2(N2n.LT.l.E-5) YD 1 F F ( N 1 ) = Y 1 ( N 1 ) 
2-Y2IN2) 

IF( AOSIYDIFF(Nl) ) . L T . ABS ( Y 1 I N I ) =» 1 . E-5 ) ) YDI FF I N 1 ) =0 .0 

RETURN 

END 


SUBROUTINE PLO T 1 ( X , Y , N , YHOR » NY , XHDR , NX , SYl, SY2,SXl, SX2,XY» 
2XSFT,YSF r ) 

DIMENSION XINI.YIN) 

DIMENSION XHDR(8)tYHDR{8) 

XIN-1)=SX1 

X(N)=SX2 

Y(N-1)=SY1 

Y(N)=SY2 

CALL PLOT! XSFT,YSFT,-3) 

CALL AXIS ( 0.0, 0.0, YHDR, NY, 8.0,90.0, SYl ,SY2) 

CALL AXIS(0.0,XY,XHUR,NX,5.0,0.0,SX1,SX2) 

Nl=N-2 

CALL LINE(X,Y,N1, 1,0,1) 

KPLOT=KPLOT+ I 

RETURN 

END 
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SURROUTINE PAIR 

REAL lOIFF, IMAXI, iMlNlf IMAX2, IMIN2t IMAX, IMIN 
COMMlN/PAIPl/TWl ,TW2iOTW,FWl,Fw2,OFWl#OFW2,DFWf TMAXO.DFMO. 
2F01FF,T01FF 

CnMMnN/PAIR2/FMAXl , TFMXI iFMINl, TFMNlf 
2 FMAX2f TFMX2fFMIN2,TFMN2 

C□MM0N/PAIR3/AF^•AXf TFMAXt AFMAXT.TFMAXT 
coMMnN/nuri /FniF iG»TniFiG,niT,AoiT 
DIMEMSI?N FDIFF(400),Tr)lFFI400) 

CQMMON/nuT2/nFlOOK,T100K 
CaMMON/TOFF/OFTni , DFT02f TOFTni, r0FT02 
FMAX=FniFF( 1) 

FM IN=FOIPF ( I ) 

FMAXl = Fn IFF ( 1 ) 

FMINl=FDIFF(l) 

TFMXl=TOlFFI n 
TFMMl = Ti)IFF ( 1 ) 

T = A^^IN1( TWltTW2) 

DG 6 I=2fN 
K=! 

IFITOIFF ( I )-T) 7,7,a 

7 FMAX=AMaX1 ( F0IFF( I ) .F'^AX) 

IFIFMAX.GT.FyAXn TFMXl = TDIFFm 
FM AX1=FMAX 

FMIN=AM,inK FDIFFI I I,FMIN) 

IFIFMIN.LT.FMIt^l) TFMNI*TDIFF( I ) 

FM^^J]=FMI^J 
6 rONTlNUE 

8 F'^AX=FniFF(K ) 

FMIN=FDIFF(K) 

FMAX2=FDIFF(K) 

FM1N2=FDIFF(K) 

TFMX2=TDIFF(K) 

TFMN2=TDIFF(K) 

00 9 I=K,N 

FMAX=AMAX1 I FDIFFI I ) ,FN1AX) 

IF (FMAX.GT.FMAX2) TFMX2= TO I FF ( I | 

FMAX2=FMAX 

FMIM=AMINl( FDIFF( I ),FVIN) 

IFIFMIN.LT.FMIN2) TFMN2= TO I FF ( 1) 
r'^IN2 = FMlN 

9 CONTINUE 

AFMAXl = A9S( F>^AX1 I 
AFMIN1 = ABS1 FNIINI ) 

IFIAFMAXI.GE.AFKINU TFMAX=TFMXl 
IFIAFMlNl.GT.AFMAXn TFMAX=TFMNl 
AFMAX=AMAX1 ( AFMAXI ,AFMINl) 

AFMAX2 = ABS( FMAX2 » 
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C 

C 

C 

C 

C 

c 

c 

r 

c 

c 

c 

c 

c 

c 

c 

c 

r 

c 

c 

c 

c 

c 

r. 

c 

c 

c 

c 

c 

c 

c 

c 

r. 

r. 

r. 

c 

c 


AFMIN2=ABS( FMIN2) 

IF (AFMAX2.GE.AFKIN2 > TFMAXT=TFMX2 
IF { AFMIM2.GT.AF.VAX2 ) TFMAXT = TFMN2 
AFMAXT=AMAX1 UFMAX2t AFMI N2 ) 

DTW=ABSIDTW) 

OFw=AtJS( DFW> 
nFWl = ABS< OFWl ) 

DFw2 = ABSIOrv.2 ) 

DFNia^ABSI OFVO) 

FOIFIG=ARS( FOIF IG) 

OF 100K=ABS(DF100K J 

4:4 « 4 * 44444 # 444444 >( 1 ' 44444444444444444*444 444444 4 444 44 44 44 44 


4 4 

♦ OUTPUT MOTOR PAIR DATA 4 

4 4 

* FMAXl ,FMIN1 ,TFMXI AMO TFMNl ARE THE MAXIMUM AND MINIMUM 4 

4 VALUES OF THRUST IMBALANCE DURING EWAT AND THE TIMES 4 

4 at which they OCCUR RESPECTIVELY 4 

4 FMAX2iFMIN2iTFvx2 AMO TFMN2 ARE THE MAXIMUM AND MINIMUM 4 

4 VALUES OF thrust IMBALANCE DURING TAIL-OFF AND THE TIMES 4 

4 AT WHICH THEY OCCUR RESPECTIVELY 4 

4 TDFTOlf TDFT02 AMO OTW ARE THE WEB TI.MES FOR THE FIRST AND 4 

4 SECOND MOrOPS TO BEGIN TAILCFF AND THE ABSOLUTE VALUE 4 

4 OF THF DIFFERENCE IN WEB TIMES RESPECTIVELY 4 

4 FWl,FW2 AND DEW ARE THE, THRUSTS AT WEB TIME EUR THE FIRST 4 

4 AND SECON!) MOTORS TO BEGIN TAILOFF AND THE ABSOLUTE 4 

4 VALUE OF THE DIFFERENCE IN THRUSTS AT wEB TIME 4 

4 RESPECTIVELY 4 

4 DFTOl and DFTP? ARE THE ABSOLUTE VALUES OF THE THRUST 4 

4 IMBA! ANCES WHICH EXIST WHEN THE FIRST AND SECOND MOTORS 

4 BEGIN TAILOFF RESPECTIVELY 

1000 CONTINUE 

4 DEMO AND TMAXQ ARE THE ABSOLUTE VALUE. OF THE THRUST * 

4 IMBALANCE V.HEN THE MAXIMUM DYNAMIC PRESSURE OCCURS ON 4 

4 THF VEHICLE AND THE TIME AT WHICH IT OCCURS RESPECTIVELY 4 

4 AFMAX AND TFMAX ARE THE ABSOLUTE VALUE UF THE MAXIMUM THRUST 4 

4 IMBALANCE OuRlN EWAT AND THE TIME AT WHICH IT OCCURS 4 

4 RESPECTIVELY 4 

4 AFMAXT AND TFMAXT ARE THF ABSOLUTE VALUE OF THE MAXIMUM 4 

4 THRUST IMBALANCE DURING TAIL-OFF AND THE TIME AT WHICH 4 

4 IT OCCURS RESPECTIVELY 4 

4 FOIFIG and TDIFIG ARE THE ABSOLUTE VALUE OF THE MAXIMUM 4 

4 THRUST imbalance DURING THE INITIAL PART OF OPERATION 4 

4 AND THE TIME AT WHICH IT OCCURS RESPECTIVELY * 

4 DIT and ADIT are ThE THE TOTAL IMPULSE IMBALANCE AND THE * 

4 ABSOLUTE VALUE OF THE TOTAL IMPULSE IMBALANCE DURING 4 

4 TAIL-OFF * 

4 DFIOOK AND TIOOK ARE THE ABSOLUTE VALUE OF THE THRUST 4 
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C ♦ IMBALAMCE when the last M0T0« reaches lOOK AND THE * 

C « Tiwp AT WHICH IT OCCURS RESPECTIVELY » 

IFITWI-TW2) 700,700,701 

700 DFT01=DFWl 
DFT02=0FW2 
CO TO 702 

701 DFT01=DPW2 

DFT02=0FW1 ! 

FWl=FW2 

FW?=FWl 

702 CONTINUE 
TDFT01=AMIN1ITWI,TW2) 

T0FT02 = AMAXH TWl, TW2) 

WRITE(6,n 

1 FORMAT!//, 20X, 'VOTOR PAIR DATA*) 

WRITE(6,2) FMAxI, TFMX l ,FMINl,TFMNl, 

2FMAX2, TFMX2,FM IN2, TFMN2,DFTni.DFT02, 

3TDFT01 ,T0FT02,0TW, FWI,FW2,OFW,OFMQ,TMAXO, 

3AFMAX, TFMAX, AFMAXT ,TFMAXT , FO I F IG, TO I F I G, DI T , A D I T , OF I OOK , TIOOK 

2 F3RMATU3X, *FMAXl= • , IPE 1 1 . 4 , I 3X , • TF MXl = •,IPEIU4,/, 

213X,*FMINI= ' , IPE 11.4, 13X, • TFMNl= *,IPF11.4,/, 

213X,»FMAX2= ' , IPE 1 1 . 4, I 3X , • TFMX2= *,IPE11.4,/, 

213X, •FMIN2= ', IPE 11.4, 13X, *TFMM2= •,IPE11.4,/, 

213X,*0FT01= • ,IPE11.4,13X,*0FT02= *,1PE11.4,/, 

213X, * TDFT01= • , IPF 1 1 .4, 1 3X , • TDFT02= • , IPE 1 1 . 4 , 1 3 X, • DTW= *,IPE11.4, 
2/,l3X,*Fwl= », IPE11.4,13X,'FW2= * , IPE 1 1 . 4, 1 3X , • OF W= ', 

21PE1I.4,/, 

213X,*DFMQ= • , IPE11.4,13X,«TMAXQ= •,1PEII.4,/, 

213X,'AFMAX= ' , IPF 1 I . 4, I 3X , * TFMAX= ',1PF11.4,/, 

213X, • AFMAXT= • , IPF I I . 4 , I 3X , * TFMAXT= *,IPE11.4,/, 

213X,*FUlFir.= • , lPEll.4,13X,»Tr)IFIG= *,1PF11.4,/, 

213X,*OIT= • ,1PE11.4,13X, *ADIT= »,1PE11.4,/, 

213X, • 0F100K= » , IPE 1 1.4, 13X, »T100K= *,IPF11.4) 

return 

END 
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SUBROUTINE I NTRP U Y , T , N « T T . DY , I CHK ) 

DIMENSION Y(N),T(N) 

Nl=N-l 

DY=0.0 

IF(ICHK) 2,2,3 

2 DO I 1 = 1, 

IF( TT.GE. T( I) .AND. TT.LT.K I + l ) ) DY=( ( Y ( I ♦ I ) -Y ( I ) ) / ( T ( I + 1 )-T ( I ) ) ) 
2»ITT-T( I ) )+Y(l ) 

IF(OY.NE.O.O) RETURN 
I CONTINUE 

3 DO A I = 1 , N 1 

IF( TT.LE.T (I ) .AND.TT.GT.TI I + l ) ) DY={ (Y(I + l)-Y(I))/mi + l)-T(I))) 
2*( TT-n I J )+Y( I ) 

IF(DY.NE.O.O) RETURN 
A CONTINUE 
RETURN 
END 


SUBROUTINE SIGBAR(X,XI,XI2,SIGX,BX, I COUNT , N , S I G 1 , S I G2 ) 
XN=FLOAT(N) 

IFI ICOUNT.GT.2) GO TO 1 

XI2=0.0 

XI = 0.0 

1 XI2=XI2+X*^2 
XI=Xl+X 
BX=XI/XN 
XI S=XI#*2 

SIGX = S0RT { (XI2/XN)-{XIS/XN=!'«2) ) 

SIG1=SQRT(XI2/XN) 

SIG2=SQRT(XI2/(2.*XN) ) 

RETURN 

END 
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V. SAMPLE STUDY 


In keeping with the present interest in very large SRMs the case 
selected for the sample study is a 146-inch diameter motor of the type 
being considered for use on the Space Shuttle. However, the study does 
not constitute a prediction of the imbalance characteristics of the 
Space Shuttle. No effort has been made to select the most recent design 
considered for the Shuttle or to minimize or maximize imbalance, Although 
the statistical characteristics of the input variables should be reason- 
able for the most part, in some cases their selection has been somewhat 
arbitrary since the purpose of this study is merely to demonstrate the 
setup procedures, format and computational capabilities of the computer 
program. 

The SRM for the sample case has three center segments consisting of 
circular perforated grains, an aft segment with a circular perforated 
grain and a forward segment with a truncated (slotted tube) nine point 
star grain. The fixed values for the sample case are given in Table V-1. 

It will be noted that both the head end and aft end domes are represented 
by hemispherical closures which is seemingly inconsistent with the fact 
that the head end of the foremost circular perforated grain is flat. 

This is an attempt to artificially correct for the curvature of the star 
grain segment located at the extreme head end. Experience with the simpli- 
fied computer program indicates that this is a satisfactory alternative 
procedure to specification of an effective length of star grain. Of course, 
the effects of the head end closure burning surface geometry could be 
represented more precisely by specifying tabular input values, but it 
should be kept in mind that only the burning surface defined by equations 
is subject to statistical variation. In the present case, because the 
entire head end segment is consumed far in advance of tailoff, only minor 
errors in the predicted behavior during the critical tailoff portion of 
the traces should be encountered as a result of the approximation used. 

Table V-2 gives the input data for the statistical variables of the 
sample case. Included in Table V-2 are brief comments on the sources of 
the data and the method in which it is applied to the present study. 

It will be noted that in a number of cases the convention is adopted of 
taking the drawing tolerances as representing ±3 standard deviations in 
a normally distributed population of a variable. Also, where more than 
one dimension controls a variable input dimension, the standard deviation 
of the variable is taken as the square root of the sum of the squares of 
the standard deviations of the controlled variables (assumed to be nor- 
mally distributed and uncorrelated.) An example of this is the standard 
deviation of the average outside diameter of the circular perforated 
grain which is calculated based on the standard deviations of the outside 
diameter of the case, and the thicknesses of the case wall , liner and 
insulation. 
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TABLE V-1. FIXED INPUT VALUES OF SAMPLE CASE 


Options 


Grain Configuration 


lEO = 1 


INPUT = 2 


IPO = 1 


GRAIN = 3 


NUMPLT(J) 00000 


STAR = 2 


Propellant Characteristics 
L = 1367.23 
TAU = 39.740 


NT = 0 
ORDER = 1 
COP = 2 


C. P. Grain Geometry 

Basic Performance Constants 

XTZO = 0.120 

DELTAY =0.040 

II = 26 S = 1 

XOUT = 1000.00 
DPOUT = 10,000.00 

ZETAF = 0.9600 Basic Star Geometry 

TMAXQ = 60.0 

TB = 122.2 NS = 1 

HB = 130,000 

PREF = 560.00 NP = 9 

DTREE = 57.285 

PIPK = 0.00150 NN = 0 

CSTART = 0.0000380 

CSTARP = 0.0057000 

PTRAN = 0.0 

GAMP = 0.00527000 
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table v-2. input data for statistical variables of sample case 


triable 

Code 

n 

n 

Source or Comment 

BHO 

51 

0.06350 

0.0000105 

Ref. 8, p. 4.5-26 (2/3 of total 
variation ascribed to propellant 
weight) . 

A1,A2 

21 

11 

0.03655 

X3 » 0.00001 



Data cards for A1 

and A2: 1,3,5 

,2,13,16,10,12,4,7,1 



Reference Values: 

a = 0.0366045 

C^ » 0.060% 



Source: Ref. 8, pp. 4.5-24, 4.5- 

analysis scaled. 

-25, MINUTEMAN data pair 

N1,N2 

60 

0.035 



ALFA 

60 

0.0 



BETA 

60 

0.0 



roal 

51 

4.350 

0.04 

Ref. 1, Weighting error distribu- 
tion, p. 44. 





From drawing tolerance: 

DE 

51 

145.67 

0.033333 

3a « 0.100 

DTI 

51 

54.430 

0.0100 

3a » 0.030 

THETA 

60 

0.0 



alfan 

60 

11.250 



LTAP 

60 

176.5 



XT 

51 

3.0400 

0.02357 

3a » 0.05 /2 controlling dimensions 

ZO 

51 

2.41 

0.02357 

3a » 0.05 /2 

ZC 

51 

0.0000 

0.02357 

3a - 0.05 /2 

RONDCN 

53 

0.0 

0.083333 

3a “ 0.250 

RONDCH 

53 

0.0 

0.083333 

3a = 0.250 

RONDGN 

53 

0.0 

0.033333 

3a =» 0.100 

RONDGH 

53 

0.0 

0.033333 

3a = 0.100 

EXN 

51 

0.0 

0.05 

3a - 0.150 

EYN 

51 

0.0 

0.05 

3a - 0.150 
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TABLE V-2. INPUT DATA FOR STATISTICAL VARIABLES OF SAMPLE CASE (CONT'D) 


irlable 

Code 

3Q 

3^ 

Source or Comment 

EXH 

51 

0.0 

0.05 

3a = 0.150 

EYH 

51 

0.0 

0.05 

3a = 0.150 

ALPHAN 

52 

0.0 

360 

Random orientation of mandrel 
and case. 

ALPHAH 

52 

0.0 

360 

Random orientation of mandrel 
case. 

ERSEF 

51 

0.00763 

0.00032 

POSEIDON data per Ref. 8, 
p. 4.5-26, Cv = 4.2%. 

TGR 

51 

60.0 

0.2333 

Ref. 8, p. 4.5-26, 6a = 1.4*F 

TIGR 

11 

40.0 

0.3740 

X3 = 0.0040, X4 = 15 


Data card for TIGR: 0.3777, 

0.3811, 

0.4030, 

0.3980, 

0.3744, 

0.3795 

0.4266, 

0.4300, 0.4334, 

0.4300, 

0.3980, 

0.3997, 

0.3862, 

0.3845 

0.3895, 

0.3963, 0.4013, 

0.3929, 

0.4097, 

0.4081, 

0.3980, 

0.4030 

0.3827, 

0.3963, 0.3980, 

0.3996, 

0.3895, 

0.3929, 

0.4081, 

0.4132 

0.4215, 

0.4148, 0.3946, 

0.3744, 

0.3845, 

0.4136, 

0.4148, 

0.4030 


0.4013. 

Reference values: “ 0.40, C^ = 3.83%. 


Source: Artificial data based on C^ of large nxsnber of motors. 


DO 

51 

143.080 

0.01462 

3a = /0.0322 + 2 X 0.022 + q.012 

DI 

51 

63.590 

0.033333 

3a = 0.10 

THETAG 

60 

10.1990 



LGCI 

51 

1135.58 

0.577 

3a = 1.0/3 segments 

LGNI 

51 

51.20 

0.33333 

3a = 1.0 

THETCN 

60 

0.0 



THETCH 

60 

90.0 



LGSI 

51 

189.15 

0.3333 

3a = 1.0 

RC 

51 

71.540 

0.00731 

3a = /0.032^ + 2x0.02^ + 0.01"^/2 

FILL 

51 

2.010 

0.011111 

3a = (1/9 points x 0.1) 

RP 

51 

12.000 

0.01667 

3a = 0.05 

RIS 

51 

63.540 

0.01667 

3a = 0.05 
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Not only must the procedures used In manufacture and quality control 
of the motor production be recognized when specifying the Input character- 
istics, but also the way a particular variable is used in the program. 

Thus, when a dimension (or other characteristic) of a variable is subject 
to random variation and the effect of the variation is averaged in the 
program, the standard deviation in the variable should be reduced. For 
exan^jle, the standard deviation in the fillet radii of the star points 
is reduced by the because the nine star points each have equal effects 
on the burning surface. Similarly, the propellant average burning rate 
variation between pairs may be reduced substantially if propellant from 
the same mixer batch is divided between the pair of motors as was assumed 
in this sample case. 

In treatment of characteristics such as burning rate , propellant 
density, and grain temperature, it is assumed that the primary concern is 
the variation between two motors of a pair and variation between pairs is 
not considered. Thus, for example, the standard deviation in these charac- 
teristics is only that within a pair and excludes any between pair varia- 
tion as experienced from change of a lot of propellant between pairs or from 
different ambient temperature histories. 

An actual printout of the statistical input data is shown in Table 
V-3 to demonstrate the computer program format. Table V-4 gives the com- 
plete set of variables (fixed and distributed) selected by the Monte Carlo 
program for the first motor of the first pair. Table V-5 illustrates the 
format for the printout of transient values for one SRM of the sample case 
and also shows the propellant weight and initial and final seed numbers for 
the configuration. The initial seed number may be used to repeat the 
calculations for the configuration. The tabular output of motor pair 
imbalance data shown at the bottom of Table V-5 may be omitted by proper 
specification of IPO. 

Figure V-1 through V-5 illustrate the graphical data which may be 
obtained using the CalComp plotter. Any or all of the plots may be 
omitted by proper choice of the NUMPLT(J) input array (See also IPO). 

The figures are derived from the same initial seed number given in Table 
V-5. 


Tables V-6 and V-7 illustrate the statistical analysis of the output 
data that may be obtained from the program by use of the appropriate value 
of IPO. Table V-6 is a facsimile of the computer output for 25 motor 
pairs of the sample study and Table V-7 is a compilation of selected 
statistical characteristics of 50 motor pairs. The latter was obtained 
by separate calculations from three groups of data. 

A histogram (Fig. V-6) of maximum thrust imbalances during 
tailoff for the 50 motor pairs demonstrates that the program results are 
generally consistent with the type of behavior which would be expected. 



TABLE V-3. PRINTOUT OF STATISTICAL INPUT DATA FROM SAMPLE CASE 
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TABLE V-4. PRINTOUT OF INPUT VARIABLES FOR ONE SRM FROM SAMPLE CASE 


CONFIGURATION NUMBER 1 


OPTIONS 

1E0» I 
1P0= I 

NUMPI.T(J)= 0 0 0 0 0 


PROPELLANT CHARACTERISTICS 
RHO= 0.063507 
Al» 0.03665 
Nl= 0.350 
A?= 0.03667 
N2= 0.350 
ALPHA= 0.0 
8ETA= 0.0 
ROAL= 6.2S52 
CSTARN= 5.I632E 03 
GAHN>= I.IAIAE 00 


BASIC MOTOR DIMENSIONS 
L= 1367.23 
TAU= 39.7^0 
DE= I.A575E 02 
OTIx 5.AA37E 01 
TMETA= 0.0 
ALFAN= l.^SOE 01 
LfA»»= 1.7650E 02 
Xr= 3.0A70E CO 
20= 2.A1A2E 00 

ZC= 5.6571E-03 
RUNDCN= 1.9500E-01 
RONDCHs 2.5000E-02 
RONDCN= 3.6000E-02 
RONOGH= 2.0000E-02 
EXN= A.2001E-02 
EYN= 2.1001E-02 
EXH= 3.9C01E-02 
EYH= -2.0999E-02 
ALPHAN=» 2.5200E 02 
ALPHAHs 3.3L20E 02 


BASIC PERFORMANCE CONSTANTS 
OELTAY* O.OAO 
11= 26 
XUUT» 1000.00 
DPOUT» ICCOO.OO 
2erAF= 0.0600 
TB* 122.2 
HD= 130000. 

ERREF= 0.00782 
PREF* 560.00 
DTKEFs 57.285 
TGR* 59.916 
P1PK= 0.00150 
CSTART= O.0C0O380 
PTRAN= 0.0 
CSTARP= 0.0C57000 
TIG* 0.3906 
GAKP= 0.0052700 
THAXQ* 60.000 


CRAIN CONFIGURATION 
INPUT* 2 
GRAIN* 3 
STAR* 2 
NT= 0. 

ORDER* I 
COP* 2 


C.P. GRAIN GEOMETRY 
DO* IA3.067 
01= 63.575 

XT2U* 0.120 
S* 1. 

THETAG* 10-19900 
LGCI* 1135.91 
LGNI* 51.15 
THETCN* 0.0 
THETCH* 90.00000 


BASIC STAR GEOMETRY 
NS* 1. 

LCSI* 189.21 
NP« 9. 

RC* 71.527 

fill* 2.006 

NN* 0. 


TRUNCATED STAR GEOMETRY 
RP* 12.005 
RIS* 63.555 
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TABLE V-5. PORTION OF PRINTOUT OF TRANSIENT VALUES FOR ONE SRM FROM SAMPLE CASE 
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Figure V-1. Thrust versus time for one pair of SRMs of 
sample case (CalComp plot-reduced size) . 
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Figure V-2. Thrust during tailoff versus time for one pair of 
SRMs of sample case (CalComp plot-reduced size) . 
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TABLE V-6. OUTPUT VALUES FOR FINAL MOTOR PAIR AND STATISTICAL CHARACTER- 



ISTICS OF 

TWENTY-FIVE MOTOR PAIRS 

FROM THE 

SAMPLE 

STUDY. 




MOTOR PAIR DATA 




FMAX1= 

3.5491E 

04 

TFMX1= 

1.0226E 

02 



FMIN1= 

-1.7347E 

04 

TFMN1= 

8.0128E 

00 



FMAX2= 

2.0555E 

05 

TFMX2= 

1.1119E 

02 



FMIN2= 

1.5032E 

03 

TFMN2= 

1.2134E 

02 



TDFT01= 

1.1066E 

02 

TDFT02= 

1.1099E 

02 

DTW= 

3.2643E-01 

FW1= 

2.0811E 

06 

F172= 

2.0809E 

06 

DFW= 

1.7500 02 

DFT01= 

2.2095E 

04 

DFT02= 

1.1096E 

05 



DFMQ= 

6.6945E 

03 

TMAXQ= 

6.0000E 

02 



AFMAX= 

3.5419E 

04 

TFMAX= 

1.0226E 

02 



AFMAXT= 

2.0555E 

05 

TFMAXT= 

1.1119E 

02 



FDIFIG= 

1.6448E 

04 

TDIFIG= 

2.2855E 

00 



DIT= 

8.6436E 

05 

ADIT= 

8.7547E 

05 



DF100K= 

2.8583E 

04 

T100K= 

1.1869E 

02 




STANDARD DEVIATIONS AND MEANS FOR MOTOR PAIR DATA 
STD. DEV. MEAN 


AFMAX 

0.1016E 

05 

0.1952E 

05 

TFMAX 

0.3787E 

02 

0.8279E 

02 

AFMAXT 

0.6072E 

05 

0.9811E 

05 

TFMAXT 

0.5590E 

00 

0.1114E 

03 

DTW 

0.1282E 

00 

0.176 IE 

00 

FWl 

0.4726E 

04 

0.2079E 

06 

FW2 

0.5726E 

04 

0.2076E 

06 

DFW 

0.2551E 

04 

0.6369E 

04 

DFTOl 

0.6512E 

04 

0.6369E 

04 

TDFTOl 

0.1291E 

00 

0.1107E 

03 

DTF02 

0.2757E 

05 

0.2961E 

05 

TDFT02 

0.1947E 

00 

0.1109E 

03 

DFMQ 

0.2461E 

04 

0.3679E 

04 

FDIFIG 

0.4574E 

04 

0.6234E 

04 

TDIFIG 

0.1553E 

00 

0.2233E 

01 

DIT 

0.4186E 

06 

0.5927E 

05 

ADIT 

0.2386E 

06 

0.3585E 

06 

DFIOOK 

0.8192E 

04 

0.1228E 

05 

TIOOK 

0.1531E 

00 

0.1186E 

03 


ALTERNATE DISPERSION VALUES FOR THRUST IMBALANCE DATA 

SIGMA 1 SIGMA 2 

AFMAX 0.2200E 05 0.1556E 05 

AFMAXT 0.1154E 06 0.8159E 05 
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TABLE V-7. SELECTED STATISTICAL CHARACTERISTICS OF FIFTY MOTOR 
PAIRS FROM THE SAMPLE STUDY 


STANDARD 


PARAMETER 

MEAN 

DEVIATION 

Absolute value of maximum thrust Imbalance 
during web action time (AFMAX) Ibf. 

19,620.40 

9250.22 

Time pf AFMAX(TFMAX) sec. 

83.89 

36.59 

Absolute value of maximum thrust 
Imbalance during talloff 
(AFMAXT) Ibf. 

110,346.00 

61,130.86 

Time of AFMAXT (TFMAXT) sec. 

111.60 

0.93 


Absolute value of the difference In time 0.20 0.14 

at which the two motors of a pair begin 
talloff (DTW) sec. 


Absolute value of the thrust Imbalance 2954.46 3965.88 

at input time of maximum dynamic 
pressure (DFMQ) Ibf. 


-51,059.09 461,769.56 


406,400.00 237,49794 


Absolute value of thrust Imbalance 8554.64 13,469.31 

when last motor of pair reaches 
100,000 lb. thrust during talloff 
(DFIOOK) Ibf-sec. « _ 

Time of DFIOOK (TIOOK) sec. 118.66 0.29 


Algebraic value of the Impulse imbalance 
during talloff (DIT) Ibf-sec. 


Absolute value of the area between the 
thrust-time traces of the pair during 
talloff (ADIT) Ibf-sec. 
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Figure V-6. Histogram of absolute thrust Imbalance values 
for fifty motor pairs from the sample study. 


VI. CONaUDING REMAKKS 


A technique has been established for statistically investigating the 
thrust Imbalance of pairs of SRMs firing in parallel. The computer program 
based upon the analysis permits the Imbalance characteristics of a large 
number of SRM pairs to be evaluated In a reasonably short time. 

It remains to demonstrate the accuracy of the program by comparisons 
of theoretical Imbalance results with those from real SBM populations. 
Preliminary investigations of this type were conducted during the program. 
The results, although encouraging, are too incomplete to warrant reporting. 
This Is largely due to the difficulty encountered In obtaining specific 
data to define confidently the statistical distributions of Input variables 
for past rocket motors. The necessary data Is often Incomplete or not 
readily accessible. 

Additional areas for extended effort Include Investigation of methods 
for accounting for effects of radial, axial and circumferential temperature 
gradient differences between motors of a pair and incorporation of between 
pair variations of propellant characteristics into the analysis. Ability 
to treat the between pair variations would Improve the accuracy of the 
program as a device for predicting the absolute performance characteristics 
of SBMs. 

Finally, It Is noted that a large number of improvements In the basic 
accuracy of the simplified computer program presented In References 3, 4 
and 5 have been Incorporated Into the program presented In this report. 
Those using the earlier program may wish to adopt these improvements for 
the purpose of design and performance analysis of single rocket motors. 
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APPENDIX 

PERFORMANCE SENSITIVITIES 
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Figure A-1. Thrust versus time for two SKMs with propellant density p 
difference of 3%. 




Figure A-2. Thrust versus time for two SSMs with burning rate coefficient 
a- and & 2 > difference of 3%. 



Figure A~3. Thrust versus time for two SRMs with burning rate exponent 
n, and n 2 » difference of 3 %. 






Figure A-4. Thrust versus time for two SBMs with characteristic velocity 
C* difference of 3%. 






Figure A-5. Thrust versus time for two SRMs with ratio of specific heats 
Y difference of 3%. 




Figure A~6, Thrust versus time for two SKMs with nozzle exit diameter 
D- difference of 3%. 





Figure A-7. Thrust versus time for two SRMs with nozzle initial throat 
diameter D* difference of 3%. 





Figure A-8. Thrust versus time for two SRMs with nozzle exit divergence angle 
a difference of 0.3® (2.7%). 



Figure A-9. Thrust versus time for two SRMs with aft end grain tapered 
length difference of 3%. 
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Figure A- 10. Thrust versus time for two SRJfe with aft end radial displace- 
ment of tapered grain x,^ difference of 0.3 inches (9.9%). 



Figure A-11. Thrust versus time for two SRMs with radial displacement of 

the main portion of the grain due to bore taper z difference 
of 0.3. inches (12.4%). ° 
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Figure A-14. Thrust versus time for two SEMs with out-of-round conditions of the 
grain exterior at the aft end AR _ of 0 and 0.5 inches (aan= 90®). 



Figure A-15. Thrust versus time for two SRMs with out-of-round conditions of 
the grain bore at the aft end AR of 0 and 0.5 inches. 



Figure A-16. Thrust versus tine for two SRMs with out-of-round conditions of the grain bore at the 
head end ARgh of 0 and 0.5 inches. 



Figure A-17. Thrust versus time for two SRMs with eccentricities of centers of 
grain exterior and interior at the aft end e of 0 and 0.5 inches 






Figure A-18. Thrust versus time for two SRMs with eccentricities of centers 

of grain exterior and Interior at the head end exh of 0 and -0.5 Inches 






Figure A-19. Thrust versus time for two SRMs with eccentricities of centers of 
grain exterior and interior at the head end exh of 0 and -0.5 
inches for Figs. A-18 and A-19 is 0 and + 0.5, respectively) 




Figure A- 20. Thrust versus time for two SRMs with radial nozzle throat 
erosion rate difference of 20%. 



IFlgure A-21. Thrust versus time for two SRMs with bulk propellant 
temperature T difference of 3°F. 
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Flijure A— 23i Thrxist versus time for two SRMs with average grain outside 
diameter Dg difference of 0.3 Inches (0.2%). 





Figure A-24. Thrust versus time for two SRMs with average grain inside 
diameter difference of 0.3 inches (0.5%). 





Figure A-25. Thrust versus time for two SHMs with angle of aft burning 
surface 6 difference of 3" (28.8%). 






Figure A-26. Thrust versus time for two SRMs with length of circular 
perforated grain 1^,,^ difference of 3.0 inches (0.3%). 





Figure A-27. Thrust versus time for two SRMs with length of grain 
associated with 0 difference of 3X. 



Figure A-28. Thrust versus time for two SRMs with length of star grain 

difference of 3.0 Inches (1.6%). 



Figure A-29. Thrust versus time for two SRMs with radius of star grain 
fillet difference f of 0.3 Inches (2.4%). 
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Figure A-31. Thrust versus time for two SRMs with average web thickness 
of truncated star grain t of 0.3 inches (3.7%). 


